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ABSTRACT 


In  accordance  with  specified  cont  cvietural  requirements,  the 
present  program  will  analyze  any  presently  foreseeable  turbuir.achlne- 
driven  cryogenic  aystem  in  general,  employing  helium  as  the  working 

fluid . and  will  readout  size,  weight  and  performance  capability  of 

each  major  component. 

Likewise,  with  the  planned  broad  input  spectrum,  the  effect  of 
numerical  changes  in  various  selected  variables  or  input  parameters 
may  be  studied,  either  individually  or  in  combination  with  others 
as  deaired. 

The  cryogenic  load,  also  an  elective  input  variable,  may  be 
either  wholely  within  the  vapor  phase  or  within  the  gas  phase  or 
partly  within  both  phases  as  required.  These  input  data  are  selected 
from  the  attached  T-S  diagram. 

It  will  be  noted,  however,  that  the  present  analytical  procedure 
specifically  conaiders  the  porous-plate  type  heat  exchanger,  as  recently 
developed,  wherein  very  high  numerical  values  of  "thermodynamic  effectiveness" 
(  can  be  achieved. 

Upon  activating  either  one  of  two  input  codea,  the  computer  will 
evaluate  heat  exchanger  performance  based  on  use  of  either  A1-1100-F  or 
A1-3003-F  as  the  plate  material. 
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I.  INTRODUCTION 


The  primary  purpose  of  the  present  work  is  that  of  generating  an 
analytical  procedure  for  both  study  and  design  of  turbomachine-driven 
cryogenic  systems  in  general,  utilizing  helium  as  the  working  fluid. 

In  consideration  of  the  comparatively  large  number  of  inter-related 
variables  necessary  to  broad  analytical  study  of  these  systems,  and  in 
order  to  find  either  the  individual  or  collective  effect  of  these 

various  parameters  on  overall  system  size,  weight  and  performance . 

a  rather  broad  input  spectrum  was  indicated. 

With  this  approach,  the  computer  readout  not  only  indicates  the 
comparative  effect  of  changes  in  the  various  design  parameters,  but 
will  guide  the  designer  in  optimizing  a  system  design  for  a  particular 
application. 

The  present  program  is  so  arranged  that  data  specifying  the  selected 
cryogenic  load  is  accepted  as  inputs. ... .That  is  to  say,  upon  referring 
to  the  T-S  diagram  Fig.  1,  one  selects  a  set  of  temperatures  and  pressures 
that  simultaneously  defines  ths  desired  cryogenic  load.  Thus  with  this 
election,  the  selected  cryogenic  load  may  be  either  entirely  within 
the  vapor  phase,  or  within  the  gas  phase,  or  partly  within  both  areas. 

It  will  be  noted,  howevsr,  that  only  those  systems  employing  a 
jT-valve  (Jpul e-Thompson  expanoion  valve)  are  considered  within  the 
present  program,  tharefore  selected  numerical  values  for  T14  (referring 
now  to  both  Figa.  1  and  2)  must  be  well  below  the  inversion  temperature 

for  helium . i.e.,  wall  below  40°K  or  72°R . Thus  for  comparatively 

larga  cryoganic  loads  primarily  within  the  vapor  phase,  one  would  select 
T14  at  or  below  10°R. 
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In  general  therefore,  referring  again  to  Fig.  1,  the  computer 
accepts  input  data  for  a  specified  cryogenic  load,  and  with  further 
input  data  defining  various  selected  design  parameters. ... .the  machine 
then  defines  the  number  and  design  of  the  heat  exchangers  and  turbines 
necessary  to  reach  540°R  within  the  external  low-pressure  compressor 
line . i.e.,  outside  of  the  vacuum  or  cryogenic  flask. 

It  will  also  be  noted  that,  depending  upon  the  selected  numerical 
values  of  such  inputs  as  pressure  ratio,  turbine  efficiency  7C»  effec¬ 
tiveness  £  ,  and  etc.,  some  systems  may  require  three  turbines  while 
others  will  require  only  two. 

With  the  exception  of  station  designations  and  actual  numerical 

values . the  equations  and  calculation  procedure  for  HX-3,  5  and  7 

is  exactly  similar  to  that  of  HX-1. 

Likewise,  with  the  same  exceptions,  the  equations  and  calculation 
procedure  for  HX-4  and  6  is  exactly  similar  to  that  of  HX-2. 
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II.  THE  HELIUM  T-S  DIAGRAM 


In  order  to  select  numerical  load  data  required  for  inputs  to 
the  mair  computer  program,  a  portion  of  the  helium  T-S  diagram  is 
shown  in  Fig.  1.  This  map  covers  the  temperature  range  of  6  to  16°R, 

and  pressures  of  10  to  300  psia . within  which  area  very  nearly  all 

presently  specified  low  temperature  cryogenic  loads  will  occur. 

As  an  example  in  the  use  of  Fig.  1,  let  it  be  agreed  that  our 
desired  load  is  to  be  wholely  within  the  vapor  phase,  at  say  6.894°R  and 
10  paia,  and  that  our  load  exit  is  to  be  on  the  saturated  vapor  line  at 
this  same  temperature  and  pressure.  One  then  selects  a  JT-valve  inlet 
temperature  and  pressure  comprising  T14  and  P14.  The  associated  constant 
H  line  is  then  followed  down  to  the  point  of  intersection  with  the  10 

psia  line,  which  establishes  T15,  P15  and  H15  entering  the  load . Note, 

here,  that  H15  equals  H14.  Then  with  known  values  of  A  Hfg  along  the 
various  constant  pressure  lines  (see  Table  1),  H16  is  readiiv  determined. 
The  load  AHL  is  then  equal  to 

AH  -  H16  -  H15  **  BTU/lb  (1) 

L 

and  with  assigned  mass  flow  W,  the  cryogenic  load  in  watts  is 

Watts  -  1^  -  1054.54  W  Al^  (2) 

The  tabulated  values  of  H  along  both  the  saturated  liquid  and 
saturated  vapor  lines,  and  tabulated  values  of  the  heatsof  vaporization 
A  Hfg;  have  been  reproduced  from  ref.(l). 

Then  as  required  by  the  program,  all  lead  input  data  comprising 
T14,  P14,  T15,  P15,  T16,  P16  and  W,  are  then  known. 
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1 1 1 .  THE  SPECIFIC  HEAT  EQUATION 

la  helium,  between  temperatures  of  about  40  and  540°R,  and  for 
pressures  between  10  and  300  psia,  the  specific  heat  is  essentially  a 
straight  line  function  ranging  from  approximately  1.38  or  so  at  40°R, 
down  to  about  1.250  3TU/lb°R  at  540°R. 

In  that  temperature  region  between  about  6  and  40°R,  however, 
and  for  pressures  between  10  and  300  pala.....the  specific  heat  exhibits 
very  wide  excursions  and  tends  to  plot-  as  a  rapidly  varying  carpet,  with 
some  disagreement  among  various  authorities. 

In  view  of  the  considerable  complications  Involved  in  programming 
derived  specific  heat  equations  for  this  region,  and  since  the  operating 
conditions  of  heat  exchanger  HX-1  will  invariably  lie  within  this  question¬ 
able  area,  as  may  part  of  HX-2 . a  special  matrix  type  computer  sub¬ 

routine  has  been  devised,  with  which  cp  can  be  determined.  These  cp 
values  will  closely  match  the  tabulated  U  values  of  ref.  (1) . l.e.  where 


dH  -  cp  dT 


(3) 


IV  THE  PORUS  PLATE  TYPE  HEAT  EXCHANGER 


^  . 


V. 


The  porous  plate  type  heat  exchanger  considered  here  is,  as  its 
name  implies,  constructed  of  a  comparatively  large  number  of  thin 
perforated  plates  (usually  of  aluminum),  with  thin  plastic  separators 
inserted  between  each  plate . see  ref.  (5). 

These  plastic  separators,  here  assumed  as  Teflon,  are  bonded  to 
the  aluminum  plates  and,  due  to  their  comparatively  low  thermal  con¬ 
ductivity,  also  serve  to  reduce  the  endwise  heat  leakage.  Should  the 
plastic  or  other  low-thermal-conductivity  separators  not  be  employed, 
the  endwiBe  heat  loBS,by  conduction  through  the  aluminum  plates* could 
reach  such  values  that  a  high  thermal  effectiveness  a  could  not  be 
attained. 

The  present  state-of-the-art  is  such  that  practical  numerical  values 
of  thermal  effectiveness  £  can  approach  99%.  At  thiB  point,  however, 
any  further  increase  in  f  is  achieved  only  with  considerably  increased 

heat  exchanger  surface  area  and  physical  size . Thus  the  economics 

of  a  particular  application  enters  at  this  point. 

It  therefore  la  suggested  that  numerical  values  of  c  ,  whether  inputs 
or  computed,  ba  limited  to  about  .985  or  .987  throughout  comparative 
studies.  Thereafter,  should  the  computer  readout  values  for  e  ^.987, 
ona  merely  alters  the  associated  inputs  and  continues. 

Referring  now  to  Fig.  3 . Here  a  typical  heat  exchanger  plate  is 

shown  in  order  to  define  some  of  the  associated  nomenclature.  It  ie 
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to  be  noted,  however,  that  the  present  program  will  accept  input  data  for 
any  of  several  different  mechanical  configurations,  which  variations  will 
be  discussed  later. 

As  shown  in  Fig.  3,  flow  through  the  plates  occurs  through  a 
large  number  of  closely  spaced,  small  diameter  holes. 

The  generalized  mathematical  relations  between  hole  diameter  dh, 
porosity  f,  hole  spacing  s,  and  heat  t'ansfer  surface  ares  Axl  per 

plate  per  unit  of  plate  face  area. ..  .have  been  derived  as  follov.s . 

Upon  inspection  of  the  magnified  or  blown  up  area  of  Fig.  3-B,  one 
can  write 


a  *  b 

(*£) 

■  —  ■ 

a  •  b 


Now  from  the  figure. 


a  ■  2  (s  sin  tC  ) 


and 


b  -  2  (a  cos  ) 


(A) 

(5) 

(6) 

(?) 
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imM|>piTiTiirtnnr — *  »>« 


As  employed  within  the  present  program,  o C  is  assumed  and  held 
constant  at  30  degrees. 


*  a  «  b  -  *  s2  sin  30°  •  cos  30° 


-  1.732060  b> 


.  .  <r  m 


1.73206  s2 


.906894  dh2 


.906894 


Here  it  is  suggested,  for  mechanical  and  other  reasons,  that 
porosity  (The  restricted  to  the  range  of 

.40  4  <r  4  .60 

Then  the  number  of  holes  m  per  inch  2  of  plate  face  area  on  the 
hot  or  driving  side,  becomes 


if  dh2 
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and  the  theoretical  heat  transfer  surface  area  per  inch  of  plate  face 
area  is  then 

Ax  i  *  (m  Tf  dh  tp)  +2  (1-ff*)  (15) 

thereafter,  the  total  hole  area  required  for  a  given  flow  rate 

at  given  velocity,  Reynold's  Number,  etc . la  readily  found  by 

solving  the  continuity  equation  wherein 

W  V 

—  •  >°v  -  (16) 

A  & 

2 

and  total  hole-flow  area  in  inches 

144  W  nr 

ah  [  »  —  (17) 

.  V 

Therefore  total  required  face  area  (on  one  side  of  the  plate  and 

2 

for  the  hot  or  driving  side  only)  in  inches* 

ah  i 

Afl  »  -  (18) 

cr 

Particular  attention  is  directed  to  the  fact  that  area  requirements 
are  first  generated  cr  found  on  the  hot  or  "driving  side"  of  the  heat 
exchangers. 

In  addition  to  the  major  advantage  of  very  large  heat  transfer  surface 
area  per  unit  of  total  volume . 2600  ft^/f.^  or  more  being  readily 


attainable . the  comparatively  thin  plates  employed  (.004  to  .010 

inches  thick)  involve  such  hole  length/diameter  ratios  thatj  at 
Reynold's  Numbers  of  2000  and  less.. ...a  laminar  boundary  layer  will 
not  be  established. 

Under  these  conditions,  i.e.,  boundary  layer  thickness  approaching 
zero,  it  is  known  that  relative  Nusselt  Numbers  (and  heat  transfer  film 
coefficients)  can  easily  reach  numerical  values  approximately  ten  times 
that  found  in  say  a  long  tube  after  the  boundary  layer  has  been  established. 

For  Instance,  in  accordance  with  ref.  (4) 


at  constant  wall  temperature.  Reference  (3)  also  reports  similar 
findings  within  this  area. 

In  recent  exploratory  runs  with  the  present  program,  utilising  a  plate 
thickness  tp  of  .0063  Inches  end  hole  diameter  dh  of  .0080  inches . Nusselt 
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Numbers  approaching  30  have  been  achieved  at  reasonable  Reynold's 
Numbers  NRe  and  Prandtl  Numbers  NPr. 

Nov  since  the  film  coefficient 

NNu  •  K 

h -  (20) 

dh 


or 


12  NNu  •  k 

h  -  -  (2D 

dh 

(where  K  is  the  thermal  conductivity  of  the  gas)  then  the  capabilities 
of  the  plate  type  heat  exchanger  are  readily  visualized. 

Simultaneously  referring  now  to  the  input  section  for  say  HX-1,  and 
to  Fig.  3-A. ... .While  the  computer  evaluates  almost  all  of  the  physical 
plate  dimensions,  working  with  certain  "key"  ratios  and  border  dimensions 
as  inputs,  some  discussion  of  the  related  concept  is  in  order. 

In  genersl,  any  desired  numerical  value  can  be  employed  for  ths  various 
input  dimensions  and  rstios,  however,  some  degree  of  relative  compatibility 
must  also  be  observed  in  selection  of  the  various  inputs. 

For  instance,  in  apecifying  the  C  factor  and  in  order  to  maintain 
approximately  equal  heat  flcw-path-lengths  in  both  the  hot-to-cold 

direct Iona . it  is  obvious  thst  a  hot  atrip  must  alwaya  lla  batwaen 

two  cold  strips.  It  will  also  be  noted  that  an  external  or  outside  cold 
flew  strip  must  ba  one  half  the  height  cf  an  internal  cold  flow  strip . 
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fT 


i.e. ,  Y2  as  compared  with  Y3  ,  and  that  both  Y2  and  Y3  will  be 

greater  than  the  associated  hot  strip  height  Y i * . Thus  the  selected 

C  factor  should  always  be  an  odd  number,  auch  as  3,  5,  7,  or  9  eu . 

C  -  Nh  +  Nc  (22) 

Then  with  specified  inputs  C,  Nh,  Nc,  Fs,  Ra,  Bx,  By  and  B*t  the 
computer  determines  all  other  plate  dimensions,  including  the  various 
heat  flow-path  lengths. 

Upon  specifying  the  shape  factor  Fs,  the  computer  determines  all 
X  and  Y  dimensions  required  to  handle  the  specified  flow  at  the  specified 
temperature,  Pressure  and  Reynold's  Number  conditions. 

Then  with  Aft  known  from  prior  machine  calculation,  the  computer 
selects 


X  "  assume 

(23) 

X*  -  X  -  (2  Bx) 

(24) 

Y '  •  JL  -  <2  By) 

Fs 

(25) 

,  -  f(C-l)  B’l  -  Nh 

Y 1  » - - - - - 

(26) 

2  (C-l) 

calculated  »  Nh  (X  Yj') 

(27) 

The  machine  then  iterates  on  X  until  Afi  calculated  equals  Aft 
required. 
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Thereafter, 


Y2*  -  Yl  ' 

i  i 

Y3  -  Yi 


Yl' 

1 1  “  - 

24 


12 


Ra 

2 

•  Ra 


h 


Y3 


f 


12  *  Ra 


(28) 

(29) 

(30) 

(31) 

(32) 


The  computer  then  evaluates  required  flow  areas  and  various  pressure 


drops 


V.  ANALYTICAL  TECHNIQUE  FOR  THE  COUNTER-FLOW  HEAT  EXCHANGER 

Within  the  present  program  the  idea  of  the  f-Ntu  or  effectiveness- 
number  of  heat  transfer  units  approach  is  employed  in  counter-flow  heat 
exchanger  analysis,  wherein  thermal  effectiveness  is  defined  as 

Q 

(.  «  -  (33) 

Qraax 

Co  (Thl  -  Th2) 

-  - = -  (34) 

Croin  (Tht  -  Tcf) 

Cc  (Tc2  -  Tci) 

-  -  (35) 

Cmin  (Thi  -  Tcj ) 

where 

Ch  *  (W  cp)  hot  side 
Cc  «  (W  cp)  cold  aide 

and  Cmin  is  the  smaller  numerical  value  of  the  two. 

In  each  case  cj>  la  either  the  hot  or  cold  aide  specific  heat 
integrated  over  the  particular  temperature  range  of  lntareat. 

Ntu,  usually  referred  to  as  the  number  of  heat  transfer  units, 
essentially  expresses  the  "thermodynamic  size"  of  the  heat  exchanger, 
and  la  defined  by 

A  •  Uavg 

Ntu  *  -  1  ■  ■  (38) 

Cmin 


(36) 

(37) 
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1 


A 


Referring  now  to  Fig.  4  and  deriving  further  working  equations  from 

those  of  ref.  (3) . It  will  be  noted  that  on  the  temperature-temperature 

plane,  and  with  given  endpoint  temperatures  but  with  no  flow  present,  the 
thermal  equilibrium  condition  of  a  counter-flow  heat  exchanger  may  be 
representated  by  a  straight  line  having  the  slope  dY/dX  equal  to  1.0. 

It  will  also  be  noted  that  with  no  flow  present,  the  endpolnc  temperature 
gradients  and  are  equal  to  zero. 

Upon  establishing  a  finite  flow,  however,  a  steady  state  operating 
line  may  ba  defined,  which  operating  line  will  have  a  slope  dY/dX  exactly 
equal  to  (W/c^. 

With  this  Idea  in  mind,  and  considering  the  6 -Ntu  approach,  many 
useful  relations  can  be  derived,  which  relations  permit  very  rapid  and 
direct  analysis  as  compared  with  the  older  log  mean  -At  approach. 

Hare  attention  is  directed  to  the  fact,  however,  that  the  relative 
disposition  or  location  of  the  endpoint  temperature  gradients  Vi  and 
Vg  must  be  carefully  observed  at  all  times. 

Undar  these  conditions  it  can  be  readily  shown  that 


(40) 
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and  that 


v£ 

v< 


1  -  6 

1-6  (Cmin/Cmax) 


(41) 


Likewise,  referring  simultaneously  to  Fig.  5, 


and 


V/  -  V2 


(42) 


(43) 


(44) 


(45) 


Therefore,  in  Fig.  4, 

Y2  ■  Yi  +  Ay 


(46) 
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and 


then 


X2  -  Xl  +  AX 


Vj  «  Y2  -  X2 


(47) 


(48) 


and 


Vg  «  Yi  -  Xi 


Further,  returning  to  equations  (33)  through  (35), 

Q 

6  .  - - 

Qmax 

Ch  .  AY 

m  —  ■ 

Cmin  (Thl-  Tci) 


(49) 


(50) 


(51) 


Cc  .AX 

Coin  (Thl-  Tci) 


(52) 


Upon  combining  equations  (40)  and  (41) ,  one  can  arrive  at  the 
basic  relation 


-Ntu  [l-(Ctain/Cmax)j 


(53) 
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where  Ch  ^  Cc. 

In  Chose  rare  cases  where  Ch  s  Cc,  that  is  to  say,  where 


or  vice-versa,  then 


=  1.0 


1  +  Ntu 


The  foragoing,  aid  various  additional  expressions  derived  there¬ 
from,  are  employed  throughout  the  main  analytical  program. 


Attention  is  also  directed  to  the  distinction  between  thermal 
effectiveness  and  €  ,  as  employed  here. 

Since  a  very  considerable  amount  of  iteration  on  temperature  is 
avoided  by  analytically  finding  actual  or  desired  temperature  based  on 

a  thermal  effectiveness  ft . and  finding  required  surface  area  based 

on  calculation  with  a  higher  thermal  effectiveness  £  (reflecting  the 
effect  of  heat  loss).,.,,  g/  and  «  are  not  to  be  considered  as  "ideal” 
and  "actual"  a«  usually  employed. 

Here, 


.  6r - 

l  1-  * 
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where  A  reflects  the  effect  of  loss  by  endwise  heat  leakage. 
This  approach  Is  believed  to  be  valid  since  it  can  be  shown 

that 


ill.  s  4  .  > 

Q  * 

where 

*  .  [(K/M 

Ctn  in 


(58) 


(59) 


Therefore,  with  the  inverse  problem,  equation  (57)  in  effect 
demand*  a  greater  heat  transfer  surface  area,  to  reach  given  temperatures, 
in  the  presence  of  heat  loss. 
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VI.  HEAT  LEAKAGE  PATHS  NOT  CONSIDERED 


In  addition  to  the  effect  of  heat  loss  by  endwise  conduction, 
as  discussed  within  Section  IV, there  remains  four  possible  heat  leakage 
paths  which  can  not  be  rationally  considered  within  tbe  present  generalized 

analytical  program . since  an  actual  mechanical  design  or  configuration 

would  first  be  necessary,  along  with  planned  grouping  or  arrangement  of 
the  various  components . 

These  additional  heat  leakage  paths  are, 

1.  Radiation  between  the  various  system  components  and  the 
enclosing  vacuum-flask  walls. 

2.  Heat  conduction  through  all  bracing  and  support  members. 

3.  Heat  conduction  through  all  thermocouple  wire  lesds. 

4.  And  heat  conduction  through  the  turbogenerator  lead  wires 
(assuming  the  turbogenerator  load  to  be  external  of  the 
cryogenic  section). 

Regarding  radiation  effects . Once  an  actual  design  and  grouping 

of  components  has  been  established,  the  effect  of  radiated  heat  loss 
from  the  various  components  can  be  readily  evaluated  with  known  areas, 
eaissivity  and  temperature  difference  by 


(60) 


whers  t  is  the  equivalent  emissivity. 
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Regarding  heat  loss  through  bracing  and  support  members 
possible  solution,  for  minimum  heat  leakage,  ia  that  of  employing  a  strong, 
low-therraal-conductivity  material  auch  as  Teflon  in  theae  areas,  however, 
an  actual  established  design  ia  required  before  radiation  and  conduction 
loas  can  be  evaluated. 

Regarding  heat  loss  through  wires .... .Here ,  again,  wire  length,  aize 
and  material  must  necessarily  be  kn^'-m  before  heat  lose  can  be  determined. 
However,  it  can  be  noted,  at  this  point,  that  one  or  two  "cold  sinka"  may 
be  provided  at  judicious  points  along  all  "bundles"  of  wire  leads,  by 
tapping  off  coclant  from  the  main  -tring  of  heat  exchangers.  The  required 
amount  of  refrigeration  could  not  be  large. 

The  item  1  and  2  heat  loss  can  and  should  be  charged  to  affected 
heat  exchangers  after  a  particular  aystem  design  la  established,  along 
with  planned  physical  grouping.  This  can  be  accompliahed  by  small 
modifications  within  the  main  program. 

In  the  case  of  radiation  effects,  and  considering  the  equations  of 

HX-1  as  nn  example . Qr  can  be  computed  and  added  to  the  numerator  of 

such  equations  as  (1Q2) ,  (IO5)  and  (108) . which  results  in  a  larger  ^ 

and,  thereby,  HX  design  for  greater  cooling  capacity. 

In  the  case  of  heat  lost  by  conduction  through  various  bracing  and 
support  members,  another  Q  loss  can  be  determined  by  a  group  of  equations 
similar  to  Ob)  through  U01) . .  „  .This  additional  Q  loss  would  also  be 
added  Into  the  numerator  of  equations  (iQ2),  (105)  *nd  (lQ8)»  again  using 
HX-1  as  an  example. 
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VII.  PROGRAM  OPERATIONAL  PROCEDURE 


Referring  now  to  Fig.  2 . The  computer  is  instructed  to 

"stack"  or  employ  any  number  of  turbogenerators  up  to  a  maximum  of 

three,  and  any  number  of  heat  exchangers  up  to  a  maximum  of  seven . 

as  required  to  reach  the  selected  temperature  of  540°R  in  the  external 
low  pressure  compressor  line,  shown  here  as  station  23. 

In  general,  however,  two  turbogenerators  and  five  heat  exchangers 
will  usually  be  sufficient  to  reach  any  desired  cryogenic  load  temper¬ 
ature,  in  helium,  from  540°R,  On  the  other  hand,  since  it  may  be 
desirable  to  study  very  low  pressure  ratio-low  efficiency  systems  for 
some  particular  application,  the  analytical  capability  to  handle  three 
turbines  and  sever,  heat  exchangers  has  been  incorporated  within  the 
program. 

In  good  rational  system  design,  this  ambient  condition  should  be 
reached  either  at  station  21  for  two-turbine  designs,  or  at  station  23 
for  three-turbine  systems. ... .This  means  that  a  final  heat  exchanger 
should  always  be  employed  above  the  final  turbine,  since  considerable 
"cooling  capability"  is.  always  present  within  the  turbine  exit  flow 
stream. 

Therefore,  in  those  cases  where  the  desired  final  exit  temperature 
(540*R  as  employed  here)  is  either  exceeded,  or  not  reached,  with  a  given 

set  of  inputs . one  merely  alters  various  inputs  until  the  desired  exit 

conditions  ara  achieved.  It  may  ba  noted  that  in  those  cases  where  it 
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is  desirable  to  retain  certain  input  values,  a  change  in  overall  system 
pressure  ratio  will  usually  serve  the  purpose. 

In  actual  practice  it  has  been  iound  that  little  is  to  be  gained 
in  driving  or  forcing  actual  heat  exchanger  effectiveness  £  to  values 
greater  than  say  .985  to  .987  since,  at  this  point,  heat  exchanger  size 
and  weight  tends  to  increase  rapidly. 

In  "driving"  or  loading  heat  exchangers  up  to  their  economic  limit, 
however,  it  should  be  observed  that  in  those  cases  where  the  "comer" 
temperatures  are  controlled  by  turbine  thermodynamics,  i.e.  turbine 
inlet  and  exit  temperatures,  nothing  other  than  added  weight  will  be 

gained  by  adding  plates  in  an  effort  to  Increase  £ . That  is  to  say, 

where  the  turbine  establishes  such  comer  temperatures  that  reflect  a 
heat  exchanger  effectiveness  of  say  only  ,820»nothing  is  to  be  gained 
by  designing  that  particular  heat  exchanger  for  an  4  of  say  .960..... 

This  situation  can  be  encountered  st  HX-2,  HX-4  and  HX-6. 

As  a  fine  point  in  further  refinement  of  best  exchanger  design, 
attention  is  directed  to  the  "fin  effectiveness  "  which  values  are 

readout  as  a  result  of  “quations  (126) (for  HX-1,3,5,7)  and  (125) 

(for  HX-2, 4, 6) .... .Experience  to  date  has  shown  thst  best  optimization 
will  usually  occur  when  numerical  values  of  fff  are  restricted  to 

.40  ^  4  -60  (61) 

In  those  casea  where  tends  to  fall  beyond  these  limits,  one 
can  usually  force  Iff  into  this  desired  range  by  altering  the  numerical 
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value  of  the  "shape  factor"  Fs,  which  factor  is  listed  as  an  input 
for  all  HX  components.  The  net  effect  of  the  foregoing  is  that  of 
optimising  (to  a  limited  extent)  the  various  heat  flow-path  lengths 
which,  in  effect,  tends  to  utilize  plate  material  more  economically. 

It  will  also  be  obaerved  that  turbine  efficiency  per  se  will  have 
little  effect  on  overall  system  performance  at  very  low  temperature 
levela.  However,  at  higher  temperature  levels,  i.e.  for  turbine  2  (and 
turbine  3  where  employed),  higher  turbine  efficiency  will  permit  greater 
work  extraction  and  will  result  in  the  more  effective  overall  system. 
This  effect  can  be  demonstrated  as  follows . 
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Within  all  analytical  work,  referring  again  to  Fig.  2,  particular 
attention  is  directed  to  the  fact  that  T16  through  T23  can  never  be 

equal  to  or  greater  (respectively)  than  T14  through  T7 . as  a  violation 

of  the  Second  Law  would  be  involved. 

In  general,  within  the  temperature-iterative  procedures,  the  program 
ia  so  arranged  that  the  machine  will  not  select  a  numerical  value  in 
violation  of  the  Second  Law. 
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VIII.  INPUT  PREPARATION 


Referring  now  to  Appendix  I. ... .Section  I  of  the  input'  data  is 
concerned  with  the  selected  cryogenic  load,  of  which  some  numerical 
values  ar^  employed  in  calculation  while  others  are  entered  and  readout 
for  identification  and  record  purposes  only. 

For  the  JT-valve  inlet,  one  selects  desired  numerical  values  for 
T14,  P14,  H14  and  S14,  either  from  the  attached  T-S  Diagram  Fig.  1,  or 
from  the  Tables  of  ref.  (1). 

For  the  JT-valve  exit  conditions  (load  inlet),  the  a3sociated  constant 
H  line  (Fig.  1)  is  followed  down  to  the  desired  load  inlet  pressure  and 
T15,  P15,  H15  and  S15  values  are  recorded.  (Note  here  that  H15  will  equal 
H14.) 

For  the  load  exit,  one  then  selects  desired  values  for  T16,  P16,  H16 
and  S16.  Then  -*ith  theBd  values  the  load  AHL  is  determined  by 

AHL  -  H16  -  H15  -  BTU/lb  (65) 

and  actual  watts 

Lw  «  1054.54  W|  *  AHL  (66) 


Input  data  for  Section  II . Section  II  is  concerned  with  the 

input*  for  the  firat  heat  exchanger  HX-1.  In  general,  a  major  portion 
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of  the  input  data  is  related  La  plate  and  hoLe  sizing,  as  discussed 
within  Sec ti  IV  of  the  text,  and  is  self  explanatory.  It  should  be 
noted,  hovevc . ,  particularly  In  relation  to  Reynold's  Number  NRei,  that 
the  computer  is  instructed  to  handle  laminar  flow  only,  therefore  numerical 
values  initially  selected  for  NRei  should  be  restricted  to  2,000  or  less . 
In  those  cases  where  li.niting  environmental  conditions  are  encountered, 
the  computer  will  readout  a  message  requesting  a  smaller  NRei. 

Values  for  the  plate  porosity  o',  as  discussed  elsewhere  within  the 
text,  should  be  limited  to  the  range  of  ,40  to  .60. 

While  the  number  of  heat  transfer  unita  Ntui  can  theoretically  reach 
any  numerical  value,  the  computer  is  also  instructed  to  readout  a  mesa age 
calling  for  less  Ntui  when  limiting  environmental  conditions  have  been 
encountered.  While  any  value  that  the  machine  will  accept  can  be  employed, 
numerical  values  in  the  order  of  10  or  so  are  suggested  for  initial  runs. 
Thereafter,  in  order  to  drive  the  heat  exchanger  to  its  economical  limit, 
or  in  optimizing  a  heat  exchanger,  one  can  increase  these  valuea  until  thia 
or  some  other  limit  is  encountered. 

The  remaining  dimensional  and/or  shape  factors  or  ratios  C,  Nc,  Nh, 

Fs  and  Ra  are  discussed  in  the  text,  within  Section  IV.  For  Instance,  C 
is  the  sum  of  all  the  flow  strips,  i.e., 

C  -  Nc  +  Nh  (67) 

Nc  and  Nh  are  the  selected  number  of  cold  and  hot  flow  strips  or 
face  areas  (see  Fig.  3),  and  will  always  be  in  the  following  relation 
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Nh=  Nc  -  1 


(68) 


and 

N c  =  Mh  t  1  (£,9) 

The  factor  Ra  reflects  the  ratio  by  which  Y3'  is  greater  than  Yl ' . 

For  instance  Fig,  3-A,  as  drawn,  reflects  an  Ra  factor  of  4,0. 

The  shape  factor  Fa  is  numerically  selected  st  some  ^alue  greater 
than  1,0,  in  thoae  cases  where  it  is  desirable  to  have  a  rectangular 

rather  than  square  crosa-aection  heat  exchanger . without  specifying 

either  X  or  Y . In  those  cases  where  a  square  croas-section  is  desired, 

the  numerical  value  of  1.0  ia  inserted. ... .With  a  specified  value  for  Fs, 
the  computer  will  determine  the  required  X-Y  values  in  consideration  of 
the  local  environment  of  temperature,  pressure  and  flow. 

The  flow  factor  Rf  ia  not  direct ly -employed  for  calculation  purposes 
within  the  program,  but  ia  entered  as  a  record  of  the  mass  flow  ratio 
within  the  heat  exchanger.  Therefore 

W  cold  side 

Rf  - -  (70) 

W  hot  aide 

It  will  be  noted  that  for  HX-1,  3,  5  and  7,  Rf  will  always  equal 
1.0  since  the  cold  aide  flow  must  equal  that  of  the  hot  side.  Within 
HX-2,  4  and  6,  however,  the  cold  side  flow  will  always  be  jreater  than 

the  hot  side  by  an  amount  equal  to  the  associated  turbine  l low . Hence 

Rf  reflects  the  choice  of  turbine  flow  relative  to  the  hot  side  flow.  Thus 
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The  numerical  value  for  (IP-31)  "CODE,"  will  in  all  cases  be 
aelected  as  either  1  or  2.  Upon  inaerting  the  value  of  1,  the  computer 
will  assume  the  plate  material  A1-1100-F,  and  will  determine  HX-1  aize 
based  on  this  material.  Upon  inserting  the  CODE  value  2,  the  computer 
will  evaluate  the  heat  exchanger  (HX-1)  based  on  use  of  A1-3003-F 
material. 


Input  data  for  Section  III . Section  III  la  concerned  with  the 

input  data  for  heat  exchanger  HX-2  and  turbina  1. 

While  the  input  valuea  for  HX-2  and  turbine  1  may  be  numerically 
different  from  those  of  HX-1,  they  are  aimilar  in  character,  with  certain 
deletions  and  additional  data. 

Since  the  preaence  of  turbine  1  now  controls  two  of  the  HX-2  end¬ 
point  or  "comer"  temperatures,  Ntui  becomea  a  calculated  result  and, 
therefore,  can  not  be  used  aa  an  input. 

Also,  since  it  la  desirable  to  aee  the  effect  of  various  turbine 
efficiency  levels  on  overall  ayatem  performance,  turbine  efficiency 

is  added  as  an  input . In  running  the  program  it  is  suggested  that  numerical 

values  be  restricted  to  the  range  of  .30  to  ,90  or  ao. 


Since  a  turbine  "bridges"  HX-2,  it  will  be  noted  (see  Fig.  2) 
that  the  cold  side  flow  will  be  greater  than  that  of  the  hot  side, 
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therefore  three  flows  are  now  specified . Wi  for  the  hot  aide, 

Wtt  for  turbine,  and  W2  for  the  cold  side. 

No  restriction  is  placed  on  numerical  values  for  theae  flows 

other  thsn  the  following . Upon  inspection  of  Fig.  2,  it  will  be 

v  noted  that  W|  must  obviously  equal  the  flow  through  HX-1 . Wtt  can 

be  any  desired  value,  limited  primarily  by  acceptable  Reynold's  Number, 
pressure  drop,  size  of  the  equipment,  and  similar  limits  imposed  by 

all  upatream  heat  exchangera . W2  must  obviously  equal  Wi  plus  Wtl , 

In  regards  to  the  hest  exchanger  material  CODE,  either  1  or  2 
is  inserted  here  in  accordance  with  the  desired  material. 


Input  data  for  all  other  sections . Here,  it  will  be  seen  that 

all  input  data  for  HX-3,  HX-5,  and  HX-7  is  exactly  similar,  in  character, 

to  that  of  HX-1 . and  sll  input  data  for  HX-4  and  HX-6  (including  the 

respective  turbines  2  and  3)  is  exactly  similar,  in  character,  to  that 
of  HX-2  and  turbine  1. 
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IX.  CONCLUSIONS  AND  RKt'.l I >AT IONS 

From  various  exp loratory  runs  in  checking  out  the  complete 
program,  it  is  felt  that  the  porous-plate  type  heat  exchanger  is 
probably  one  of  the  few  configurations  capable  of  reaching  very  low 
temperatures,  with  reasonable  component  size  and  weight. 

It  has  been  noted,  however  that  while  the  stack  of  plates  with 
plastic  separators  are  bonded  into  a  solid  assembly  as  it  were,  these 
heat  exchangers  are  mechanically  very  weak  in  the  endwise  direction  and 

would  tend  to  separate  a  appreciable  internal  pressure  is  applied . 

Therefore,  with  the  present  state-of-the-art,  these  unita  must  be  tightly 
clamped  or  held  In  bolted  fixtures  such  that  endwlae  "blowout"  pressure 
can  be  resisted. 

Their  radial  resistance  to  Internal  pressure  Is  apparently  much 
higher  as  seen  by  inspection  of  Fig.  3,  and  they  have  been  known  to 
withstand  differential  preaaurea  of  aome  2  or  3  atmcapherea.  It  ia 
conceivable,  however,  that  the  plastic  aeparatora  would  or  could  blowout, 
at  pressure  differentials  exceeding  say  4  or  5  atmoapherea . 

Since  it  may  be  desirable  to  design  future  cryogenic  systems  with 
pressure  ratios  in  the  order  of  8  or  10/1  or  more,  it  is  recommended 
that  tests  be  made  with  a  heat  exchanger  first  wrapped  with  a  plaatlc 
layer,  e.g.  Teflon,  then  tightly  wrapped  with  c  single  layer  of  amall 

diameter,  high  strength  wire . in  addition  te  endwise  bolting  or 

clamping. 
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X.  NOMENCLATURE 


Symbol  • 

Entity: 

Units: 

A 

Area  in  general 

ins2  or  ft2 

Af  / 

Plate  face  area,  hot  side 

ins  2 

Akl 

Face  area  of  heat  loss  path 

ft2 

Axl 

Heat  transfer  surface  area  per  in2  face 
area,  hot  side 

ins2 

Axp 

Total  heat  transfer  surface  area  per  plate, 
hot  side 

ins2 

Ax  tot.hs 

Total  heat  transfer  surface  area,  hot  side 

ft2 

Av 

Total  hot  side  heat  transfer  surface  area 
per  unit  total  volume 

ft2/ft3 

ahi 

Total  hole  flow  area,  hot  side 

ins2 

ah2 

Total  hole  flow  area,  cold  side 

ins2 

Bx 

External  border  dimension 

ins 

By 

External  border  dimension 

ins 

B' 

Internal  border  dimension 

ins 

C 

Configuration  factor 

-  Nh  +  Nc 

Cc2 

Heat  flow  potential,  cold  side 

=*  W  cpcm 

Chi 

Heat  flow  potential,  hot  side 

*  W  cphm 

cp 

Specific  heat  at  constant  pressure 

BTU/lb  °R 

cp 

Mean  specific  heat 

Bill /lb  °R 

cpcra 

Mean  specific  heat,  cold  side 

BTU/lb  °R 

cphm 

Mean  specific  heat,  hot  side 

BTL'/lb  °R 

dh 

Hole  diameter 

ins 

e 

Base  of  natural  logarithm 

2.7182818  etc- 
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Symbol : 

Fs 

H 

Hf 

AH 

AHfg 

AHL 

AHt 

h 

K 

K[ 

Km 

Kp 

L 

LG 

Lw 

h 
1 2 
23 

l  c 
Nc 
Nh 


Entity: 

Shape  factor 
Total  heat 

Total,  heat  In  fluid  or  liquid  atate 

Total  heat  in  gas  state 

Total  heat  change 

Latent  heat  of  vaporization 

Total  heat  load 

Total  turbine  work 

Film  coefficient 

Thermal  conductivity 

Thermal  conductivity,  loss  path 

Mean  thermal  conductivity 

Plate  thermal  conductivity 

A  length 

Logic  gate  in  computer  program 
Total  cryogenic  load 
A  heat  flow  path  length 
A  heat  flow  path  length 
A  heat  flow  path  length 
A  heat  flow  path  length,  loaa 
Number  of  cold  flew  strips,  see  Fig.  3 
Number  of  hot  flow  strip*,  see  Fig.  3 
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Units : 

X/Y 

BTU/lb 

BTU/lb 

BTU/lb 

BTU/lb 

BTU/lb 

BTU/lb 

BTU/lb 

BTU-ft/hr- ft2  °R 
BTU-ft/hr-ft2°R 
BTU-ft/hr-ft2oR 
BTU-ft/hr- ft2 °R 
BTU-ft/hr-ft2oR 
ina  or  ft 


Watts 

ft 

ft 

ft 

ft 
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Symbol : 

rr  v.  1. 1  v : 

Uni ts: 

NNu 

Meuse It  Number 

•-h*  rh/K 

NPr 

1'randtl  Number 

cp  •  M/K- 

NRe 

Reynold's  Number 

dh  V 17C"  si 

Ntu 

Number  of  heat  transfer  units,  see  text 

A'.'/  Cm  in 

ni 

Number  of  holes  per  in^  face  area,  hot  side 

- 

np 

Number  of  pistes  in  heat  exchanger 

- 

ns 

Number  of  spacers  in  hest  exchanger 

np  +  1 

P 

Pressure 

psia 

Ap 

Pressure  difference,  actual 

psia 

Ap' 

Pressure  difference,  loss 

psia 

Pcm 

Mean  Pressure,  cold  side 

psia 

Phm 

Mean  Pressure,  hot  side 

psia 

Q 

A  quantity  of  hest 

BTU 

Qt 

A  quantity  of  heat,  loss 

BTU 

Qmax 

Maximum  available  quantity  of  heat 

BTU 

Rs 

Pace  area  ratio  per  plate 

cold  side/hot  side 

Rf 

Flow  ratio 

cold/hoc 

S 

Fn trophy 

BTU/lb  °R 

a 

Hole  spacing,  center  to  center 

ins 

T 

temperature 

°R 

Tas 

Mean  temperature,  cold  side 

°R 

Thm 

Mean  temperature,  hot  side 

°R 

Til 

Mean  temperature,  cold  end 

°R 
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Symb  u  1  : 

Entity: 

Un  its: 

m 

Mean  temperature,  liot  end 

UR 

Tv 

Average  temperature,  hot  to  cold  side 

°R 

IP 

Elate  thickness 

ins 

is 

Sparer  thickness 

ins 

U 

Overall  heat  transfer  coefficient 

BTU-ft/hr-ft 

V 

Velocity 

f t/sec 

AT 

Specific  volume 

ftYlb 

W 

Flow  in  general 

Ibs/sec 

Wi 

Flow  at  various  points  in  system, 

see 

Fig.  2 

lbs/sec 

W2 

Flow  st  various  points  in  system, 

see 

Fig.  2 

lbs/sec 

W3 

Flow  st  various  points  in  system. 

aee 

Fig.  2 

lbs/ sec 

W4 

Flow  at  vsrlous  points  in  system. 

see 

Fig.  2 

lbs/sec 

Wtf 

Turbine  1  flow 

lbs/sec 

Ut2 

Turbine  2  flow 

lbs/sec 

Wt3 

Turbine  3  flow 

lbs/ sec 

X 

External  dimension,  see  Fig.  3 

ins 

X* 

Internal  dimension,  see  Fig.  3 

ins 

Ax 

Temperature  difference,  cold  side 

°R 

Y 

External  dimension,  see  Fig.  3 

ins 

Y* 

Internal  dimension,  see  Fig.  3 

ins 

Yi’ 

Internal  dimension,  see  Fig.  3 

ins 

Y2  * 

Internal  dimension,  see  Fig.  3 

ins 

Y3* 

Internal  dimension,  see  Fig.  3 

ina 

Ay 

Temperature  difference,  hot  side 

°R 
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vtnbol : 


v? 


i 


72 

Z 

co 


OC 

7 

? 

e 

g; 

V 

A 

tT 

cr 


Entity;  Unj.  t  : 

Turperature  gradient,  hot  end,  see  Figs.  A  6  5  CR 
I'erp^  future  gradient,  cold  end,  see  Figs.  4  &  5  °R 
Maximum  temperature  difference  °R 

Infinity 


Creek: 


An  angle  degrees 

Specific  heat  ratio  cpA'v 

Mean  specific  heat  ratio  cp/cv 

Actual  heat  exchanger  effectiveness,  with  loss  Q/Qmax 

Ideal  exchanger  effectiveness,  no  loss  Q/Qmax 

Fir.  effectiveness,  see  equation  (124)  in  HX-1 

_2 

Turbine  efficiency  Z  x  ]0 

A  heat  loss  function 

Viscosity  lbs/ft  sec 

Circumference  to  diameter  ratio  of  the  circle  3.14154 


Porosity  hole  area/ir.“face  area 


Subscripts: 

1  Location  identification,  aee  various  Figures 

2  Location  identification,  aee  various  Figures 

3  etc  Location  identification,  aee  various  Figures 

avg  Average 

c  Cold 

h  Hot 
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I  '•< 1  H 

•  tkl  • '  *  a 

'li.:  n:.o dynamic  properties  ct  helium  from 

.  c  54UDR  between  10  and  1500  psia. 

:■ jb  le..>  o  i  <-al  Note  154-A,  Jan,  1962, 

*  .  W.,  and  bu’xe,  W.  M. 

,  -yogenie  engineering. 

Wiley,  196.!. 

3.  Kays,  W.  M. ,  and  London,  A.  L, 

Compact  heat  exchangers. 

McGraw-Hill.  2nd  ed.,  1964. 

4.  Knudsen,  I.  G.,  and  Katz,  D,  L. 

Fluid  dynamics  snd  heat  transfer. 

McGraw-Hill,  1958, 

5.  Fleming,  R.  B . in 

Advances  in  cryogenic  engineering. 

Plenum  Press,  Vol.  14,  1969. 


36 


XII.  FIGURES 

AND  TABLES 

Figure 

Title 

Page 

1 

Temperature-Entropy  Diagram  for  Helium 

39 

2 

Cryogenic  System  on  the  T-S  Plane,  External 
Compressors  Not  Shown 

41 

3 

Generalized  Plate  Configuration  and  Dimensions 

42 

4 

Counter-Flow  Heat  Exchanger  on  the 
Temperature-Temperature  Plane 

44 

5 

Counter-Flow  Heat  Exchanger  on  the 

Temperature-Area  Plane 

45 

Table  I 

Heats  of  Vaporization 

46 
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warn* 


J.4 


39 


2.0 


GMS/CNT 


7 


F/$.  Z  C&Y0MMC  jvjrsp?  TH£  T-S  XlftNE. 
£K7tKNfiL  C0MXXXJ9OX5  HOT  MOM 


hott  I  rum 


7c, 


APE  A 


N/G.  S  COUNTER -F10W  HEAT  EXCHANGE#  ON  WE 
TEMPERATURE -AREA  WANE 


TAliLt  1  HKATS  OF  VAPORIZATION: 


T°R: 

l* 

psia: 

SI: 

Hg: 

Sg: 

AHfg: 

6.0 

5.58 

2.70 

.630 

12.81 

2.315 

10.11 

6.2 

6.44 

2.81 

.653 

12.91 

2.271 

10.10 

6.4 

7.35 

3.02 

.675 

12.99 

2.231 

9.97 

6.6 

8.35 

3.20 

.700 

13.06 

2.195 

9.86 

6.8 

9.44 

3.38 

.725 

13.13 

2.160 

9.75 

7.0 

10.64 

3.58 

.750 

13.17 

2.125 

9.59 

7.2 

11.93 

3.79 

.778 

13.21 

2.095 

9.42 

7.4 

13.31 

4.01 

.805 

13.23 

2.050 

9.22 

7.6 

14.80 

4.28 

.831 

13.23 

2.010 

8.95 

7.8 

16.43 

4.56 

.861 

13.21 

1.970 

8.65 

8.0 

18.13 

4.84 

.890 

13.15 

1.928 

8.31 

8.2 

19.90 

5.16 

.902 

13.05 

1.886 

7.89 

8.4 

21.85 

5.50 

.957 

12.90 

1.841 

7.40 

8.6 

23.91 

5.88 

.996 

12.72 

1.794 

6  -  84 

8.8 

26.28 

6.34 

1.042 

12.45 

1.738 

6.11 

9,0 

28.53 

6 ,38 

1.098 

12.05 

1.672 

5.17 

NOTES. 

.....All  tempeiature  » 

°R 

All  pressure  * 

psia 

All  H 

at 

BTU/lb 

All  S 

. 

BTU/lb  -  ®R 

From  ref.  (1) 

Use  with  Fig.  1  in  preparing  load  Inputs. 
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APPENDIX  I 


INPUT  DATA 

Note .....  Readout  all  inputs  for  record  purposes. 


I.  LOAD 

DATA: 

At  JT 

-valve 

inlet: 

Input 

No : 

Entity: 

Numerical 

Value : 

Units : 

(IP-1) 

T14 

Temperature 

°R 

(IP-2) 

P14 

Pressure 

psia 

(IP-3) 

H14 

Enthalpy 

BTU/lb 

(IP-4) 

S14 

Entropy 

BTU/lb°R 

At  JT- 

-valve 

exit: 

(IP-5) 

T15 

Temperature 

°R 

(IP -6) 

P15 

Pressure 

psia 

(IP-7) 

H15 

Enthalpy 

BTU/lb 

(IP-8) 

S15 

Entropy 

BTU/lb*R 

At  load  exit: 

(IP-9) 

T16 

Temperature 

°R 

(IP-10) 

P16 

Pressure 

psia 

(IP-11) 

H16 

Enthalpy 

BTU/lb 

(IP-12) 

S16 

Entropy 

BTU/lb°R 

(IP-13) 

Ahl 

Load 

BTU/lb 

(IP-14) 

Lw 

Load 

Watts 
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J  H_HX-U 

)ATA: 

Inpiu 

Nck_ 

Entity: 

Numerical 

Value: 

Units : 

< IP-15) 

dh 

Hole  diameter 

inches 

(IP-L6) 

tp 

Plate  thickness 

inches 

(IP-17) 

ts 

Spacer  thickness 

inches 

(IP-18) 

cr 

Porosity 

hole  ares/face  sres 

(IP-19) 

NRe  l 

Reynold's  Number,  hot  side 

dh  VfrtTM. 

(IP-20) 

Ntui 

Number  of  heat  transfer  units 
hot  side 

» 

• 

(IP-21) 

C 

Configurstion  factor 

(IP-22) 

Nh 

Number  of  hot  flow  strips 

(IP-23) 

Nc 

Number  of  cold  flow  strips 

see 

text 

(IP-24) 

Fs 

Shape  factor 

(IP-25) 

Ra 

Face  ares  rstio 

(IP-26) 

Rf 

Flow  ratio 

1 

t 

(IP-27) 

Bx 

External  border  dimension 

Inches 

(IP-28) 

By 

External  border  dimension 

inches 

(IP-29) 

i 

B 

Internal  border  dimension 

inches 

UP-30) 

W| 

Flow 

lbs/ 

sec 

(IP-31) 

CODE 

Plate  material 

see 

note 

HI.  HX-2  AND  TURBINE  1  DATA: 


(IP-32) 

dh 

Hole  diameter 

inches 

(IP-33) 

tp 

Plate  thickness 

Inches 
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Input 

No: 

Entity: 

Numerical 

Value: 

Units : 

(IP-34) 

ta 

Spacer  thickness 

inches 

(IP-35) 

<r 

Poroaity 

hole  area/face 

(iP-36) 

NRei 

Reynold’a  Number,  hot  side 

dhVAT-* 

(IP-37) 

C 

Configuration  factor 

(IP-38) 

Nh 

Number  of  hot  flow  strips 

(IP-39) 

Nc 

Number  of  cold  flow  strips 

see 

text 

(IP -40) 

Fa 

Shape  factor 

(IP -41) 

Ra 

Face  area  ratio 

(IP-42) 

Rf 

Flow  ratio 

1 

* 

(IP-43) 

Bx 

External  border  dimenaion 

inches 

(IP-44) 

By 

External  oorder  dimension 

inches 

(IP-45) 

b’ 

Internal  border  dimension 

inches 

(IP-46) 

Wl 

Hot  side  flow 

lbs/sec 

(IP -47) 

Wtl 

Turbine  flow 

lbs/sec 

(IP-48) 

W2 

Cold  side  flow 

lbs /sec 

(IP-49) 

Turbine  efficiency 

%  X 

CM 

1 

O 

rH 

(IP-50) 

CODE 

Plate  material 

•  - 

see 

note 

IV.  HX-3 

DATA: 

(IP-51) 

dh 

hole  diameter 

inches 

(IP-52) 

tp 

Plate  thickneaa 

inchea 

(IP-53) 

ts 

Spacer  thickneaa 

inches 

area 
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Input 

No: 

hut  it 

Numerical 

Value: 

Un its : 

(IP-54) 

<r 

I’m  os  I  ty 

hole  area/face  area 

( IP-55) 

NRei 

Reynold’s  Number,  hot  side 

dhV/OTM 

(11-56) 

Ncui 

Number  of  heat  transfer  units, 
hot  side 

(IP-57) 

C 

Configuration  factor 

(IP-58) 

Nh 

Number  of  hot  flow  strips 

(IP-59) 

Nc 

Number  of  cold  flow  strips 

( IP-60) 

Fs 

Shape  factor 

see 

text 

(IP-61) 

Ra 

Face  area  ratio 

(IP-62) 

Rf 

Flow  ratio 

’ 

f 

(IP-63) 

Bx 

External  border  dimensions 

inches 

(IP-64) 

By 

External  border  dimension 

Inches 

(IP-65) 

1 

B 

Internal  border  dimension 

inches 

(IP-66) 

Vi  2 

Flow 

lbs /sec 

(IP-67) 

CODE 

Plate  material 

see 

note 

V.  HX-4  AND  1URBINE  2  DATA: 


(IP-68) 

dh 

Hole  diameter 

inches 

(IP-69) 

tp 

Plate  thickness 

inches 

(IP-70) 

ts 

Spacer  thickness 

inches 

(IP-71) 

<r 

Porosity 

hole  area/face  area 

(IP-72) 

NRl*  l 

Reynold’s  Number,  hot  side 

dh  V  (ATM 
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Input 

No: 

Entity: 

Numerical 

Value : 

Units : 

(IP-73) 

C 

Configuration  factor 

| 

l 

(IP-74) 

Nh 

Number  of  hot  flow  strips 

1 

(IP-75) 

M 

Number  of  cold  flow  strips 

see  text 

(IP-76) 

Fs 

Shape  factor 

(IP-77) 

Ra 

Face  area  ratio 

(IP-78) 

Rf 

Flow  ratio 

f 

(IP-79) 

Bx 

External  border  dimension 

inches 

(IP-80) 

By 

External  border  dimension 

inches 

(IP-81)  . 

b' 

Internal  border  dimension 

inches 

(IP-82) 

W2 

Hot  side  flow 

lbs/sec 

(IP-83) 

Wt2 

Turbine  flow 

lbs/sec 

(IP-84) 

W3 

Cold  side  flow 

lbs/sec 

(IP-85) 

Turbine  efficiency 

%  X 

;  10~2 

(IP-86) 

♦CODE 

Plate  material 

see 

note 

VI.  HX-5  DATA: 


(IP -8 7) 

dh 

Hole  diameter 

inches 

(IP-88) 

tp 

Plate  thickness 

inches 

(IE-89) 

ts 

Spacer  thickness 

inches 

(IP-90) 

<r 

Porosity 

hole  area/face  area 

(IP-91) 

NRe  i 

Reynold's  Number,  hot  side 

dh  VfTirM 
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Input 

No: 

Un  t  : 

Numerical 

y:  Value: 

(IP-92) 

Ntui 

Number  of  heat  transfer  units, 
hot  side 

(IP-9  3) 

C 

Configuration  factor 

(IP-94) 

Nh 

Number  of  hot  flow  strips 

(IP-95) 

Nc 

Number  of  cold  flow  strips 

(IP-96) 

Fs 

Shape  fsctor 

(IP-97) 

Ra 

Face  srea  ratio 

(IP -98) 

Rf 

Flow  ratio 

(IP-99) 

Sx 

External  border  dimension 

(IP-100) 

By 

External  border  dimension 

(IP-101) 

» 

B 

Internal  border  dimension 

(IP-102) 

W3 

Flow 

(IP-103) 

CODE 

Plate  material 

VII.  HX-6 

AND  TURBINE  3  DATA: 

(IP-104) 

dh 

Hole  diameter 

(IP-105) 

tP 

Plate  thickness 

( IP-106) 

fs 

Spicer  thickness 

(IP-107) 

a ' 

Porosity 

(IP- 108) 

NRei 

Reynold's  Number,  hot  side 

(IP -109) 

C 

Configuration  factor 

(IP-110) 

Nh 

Number  of  hot  flow  strips 

Units: 


see  text 

> 

inches 
inches 
inches 
lbs/sec 
see  note 


inches 

inches 

inches 

hole  sres/face  sres 
dh  Vf4TM 
see  text 


$2 


see  text 


Input 

No: 

Enti  t; 

r. 

Numerical 

Value:  Units: 

(IP-111) 

Nc 

Number  of  cold  flow  strips 

(IP-112) 

Fs 

Shape  factor 

(IP-113) 

Ra 

Face  area  ratio 

see  text 

(IP-114) 

Rf 

Flow  ratio 

J 

(IP-115) 

Bx 

External  border  dimension  - 

inches 

(IP-116) 

By 

External  border  dimension 

inches 

(IP-117) 

B’ 

Internal  border  dimension 

inches 

(IP-118) 

W3 

Hot  side  flow 

lbs/ sec 

(IP-119) 

Wt3 

Turbine  flow 

lbs/sec 

(IP-120) 

W4 

Cold  side  flow 

lbs/eec 

(IP-121) 

'it 

Turbine  efficiency 

X  x  10-2 

(IP-122) 

CODE 

Plate  material 

see  note 

Vlll .  HX-7 

DATA: 

(IP-123) 

dh 

Hole  diameter 

inches 

(IP-124) 

tp 

Plate  thickness 

inches 

(IP-125) 

ts 

Spacer  thickness 

inches 

(IP-126) 

o' 

Porosity 

hole  area/face  area 

(IP-127) 

NRe  J 

Reynold's  Number,  hot  side 

dhV/^-A 

(IP-128) 

Ntui 

Number  of  heat  transfer  units, 
hot  side 

1 

(IP-129) 

C 

Configuration  factor 

see  text 

(IP-130) 

Nh 

Number  of  hot  flow  strips 

i 
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Input 

Ntmerical 

No: 

Entity: 

Value : . 

Units: 

(IP-131) 

Nc 

Number  of  cold  flow  strips 

(IP-132) 

Fs 

Shape  factor 

(IP-133) 

Ks 

Face  ares  rstlo 

1 

sec 

i  text 

UP0134) 

Rf 

Flow  ratio 

r 

(IP-135) 

Bx 

External  border  dimensiv.”. 

inches 

(IP-136) 

By 

External  border  dimension 

inches 

(IP-137) 

B' 

Internal  border  dimension 

inches 

(IP-138) 

W4 

Flew 

lbs/sec 

(IP-139) 

CODE 

Piste  material 

see  note 

Note . For  (IP-31-50-67-86-103-122  and  139), 

Insert  ’'l"  for  AL-UOO-F,  or 
Insert  "2"  for  AL-3003-F 
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APPENDIX  1! 


HX-1  PROGRAM . CALL  HX(J) 

j  *  1.0 


Inputs; 

Call  numerical  values  from  APPENDIX  I,  Sections  I  and  II. 


Initial  numerical  assumptions: 


Equa. 

Initial 

Na. 

value: 

(4) 

150.0 

(5) 

T14+10 

(6) 

1.02  P14 

03) 

2.0 

(35) 

T16+10 

(36) 

.98  P16 

Notes: 


Readout  last  result  af  all  equations  marked  with  a  star 
Do  not  stop  machine  at  TEST  9. 


HOT  SIDE: 


S  -  A/*  906894  Jhl 

v  or 


n 


4  cr 

~7Tdh* 


(1)* 

(2)* 


Ax  I  =  n  7Tdh  tpj  +  ''(l-tf) 


(3)* 


np  -  assume 

(Will  be  altered  only  by  TEST  8) 


(4)* 
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(15)* 


Y'  =  21  (2  By) 

Ft  ~ 


v.t.  Y‘  -  [  (C-i )  5'1  -  Nh 

rfc-lf  — 


06)* 


Afi  calculated  *  Nh  (X'  Y'i) 


07) 


TEST  1 

(17)  mutt  *  (12)  ±.001 

If  (17)  »  (12),  reduce  (13)  and  Iterate  from  (13). 
*  "  ,  Increase "  "  "  "  » 


Y1 2  *  Y'i  Kq 
2 


(18)* 


Y's  -  Y'i  •  Ra 

(19) 

;  Y'' 

6  "“2T 

(20) 

h  -  — 

* 2  12 

(21) 

,  Y'3 
<3  mTnra — 

(22) 

A  P',  .  370  x  Ip-6  V132  ^/Ttp/dh| 

^13  V—Nt.| 1 

* 

(23) 
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A  P|  *  np  •  AP'i 


(24) 


Pi3  calculated  =*  P 14  +  A P,  (25) 

TEST  2 

(25)  must  53  (6)  ±.001 

If  (25)  »  (6),  increase  (6)  and  iterate  from  (6). 

"  "  «  ",  reduce  "  M  . . 


Thm  “ 


Phm  -  .W*+™L  (29) 

(30) 


Chi 

*  cphm  •  Wi 

(31)’ 

NPn 

at  cphm*  ml 

(32)’ 

Kmi 
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04 


h  =  NNui  •  Km  i 
1-2  3T5  ’ 


COLD  SIDE 

T 1 7  =  assume 

(Will  be  altered  only  by  TEST  5) 

P17  *»  assume 


(Will  be  altered  only  by  TEST  3) 


Hr  17  =  2.863092 


Mi7  =  2.37008  x  10'7  (T17)  ,643 


NReZ 


V17  -  dh 
A?  17  *  jx  17 


370  x  10-6  V172 


(tp/dh) 

"Sfcr 


A  P,  -  ;‘.p  •  A  P', 


P17  calculated  s  P16  -  A  P2 

TEST  3 

(45)  must  »  (36)  ±.001 

If  (45)  >  (36),  increase  (36)  and  iterate  from  (36). 

"  "  <  ",  reduce  "  "  "  " 


TEST  4 

If  (45)  <  10,  stop  &  reodout  message. 
"  "  2  *  ,  continue. 


"REDUCE  NRei ,  or 
.  INCREASE  Ra  " 


_  8.55497  x  10 
^  (tl/-T16) 


-«  J  [(TIT)  1  “3-  (TW1'643] 


Km  2  = 


57.79  x  10"3  [  (.00355  T17)1,642-  (.  00355  T16)1  ‘ 642  ] 

}0355  (T17  -lfi6)  1  1.642 
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If  Chi  <  Cc2  If  Chi  >  Ccz  If  Chi  *  Cc2 

go  to  (55)  go  to  (61)  go  to  (67) 

!  *  * 


7,  =  (T14-T16)  e  Ntui  [l-(Chi  /Cc2 )  ] 


T17  calculated  *  T16  +  AX 

TEST  5-A 

(58)  must  =  (35)  ±.001 

If  (58)  >  (35),  Increase  (35)  and  Iterate  from  (35). 


T13  calculated  =  T1 4  +  AY 


TEST  6-A 

(59)  must  ■  (5)  ±.001 

If  (59)  >  (5),  increase  (5)  and  iterate  from  (5). 
"  "  <  ",  reduce  "  "  "  "  "  . 


-Ntu 


■'  E’-(Ch'/ c.0] 

--^tp-7a,/Ccl 


Then  ga  ta  (71) 


If  Chi  >  Cc2 


Vi  =  (T14-T16)  e 


N,“'['-(c‘*/ch,)] 


AX  *  V  1 

-  (T14  -  T16) 

Chi 

Cc2 

-  1 

AX 


T17  calculated  =  T16  +  AX 


TEST  5-8 

(64)  must  *  (35)  ±.001 

If  (64)  >  (35),  increase  (35)  and  iterate  from  (35). 
"  ••  <  "  ,  reduce  "  "  "  "  " 


(59) 


(60)* 


(61) 


(62) 


(63) 


(64) 


T13  calculated  -  T 1 4  +  AV 


(65) 


TEST  6-B 

(65)  mutt  =  >.i)  ±.001 

If  (65)  >  (5),  increase  (5)  and  iterate  from  (5). 
"  "  <  ",  reduce  "  "  "  "  "  . 


c  i  = 


1  -e 


1- 


Cc2 

‘'"CKT 


Ntu!  [  l-(Cca/ch|)] 

^03 


-  Ntu 


T\ 


(66)* 


Then  ga  ta  (71) 


If  Chi  S  Cc2 


±  10 


Z  = 


(T14-T16) 


(67) 


T17  calculated  *  (T16  +  Z)  -  (T14  -  T16)  (68) 


TEST  5-C 

(68)  mutt  =  (35)  ±.001 

If  (68)  >  (35),  increase  (35)  and  Iterate  from  (35). 
"  "  <  ",  reduce  "  "  "  "  " 


T 1 3  calculated  =  T16  +  Z 


(65) 


TEST  6~C 

(69)  mutt  a  (5)  ±.001 

If  (69)  >  (5),  increote  (5)  and  iterate  from  (5). 
"  "  <  "  ,  reduce  "  "  "  "  "  . 


If  0  <  Thm  <  113 
go  to  (73) 


If  113  <  Thm  <540 
go  to  (74) 


tf  540  <  Thm  <  oo 
go  to  (75) 


Kp  = 


[.1  (Thm 


“>]{t  [<•’  TW2-(.l  T'.  )2]  -  _^_[(.lTh4'4-(-'T''1'47]j 


(73) 


Th*n  go  to  (79) 


Kp  "[.1  (Thm -T'i  )]' 


f 

'  r  2.708  2.7081' 

2  2-i"| 

[(.1  Thm)  »  (.1  T'i  )  J 

+ 

9.551  [(.1  Thm)  -  ( .  1  T'i)  ] 

2.708 

l 

Thon  go  to  (79) 


(74) 


Kp  =  111. 74  =  constant 


(75) 


r 

If  0  <  Thm  ^  210 
go  to  (76) 


Thon  go  to  (79) 


go  to  (77)  go  to  (78) 


Than  go  to  (79) 


♦ 


If  (IP -31)  CODE  =1.0 


* 


If  (IP  -  31)  CODE  =2.0 


go  to  (87) 

* 


1 

£■1  (Tfe-  Tcm) 


go  to  (88) 

go  to  (89) 

(t 

r  2  2 

f  3.47  3.47* 

(.1  T’z)-  (.1  Tcm) 

-  1 

(.1  T4)  -(.1  Tcm) 

1 2 

l  -1 

TZ7 

L 

(87) 


Then  go  to  (93) 


2 


t/,.  snim 


AM  = 


"  r  * 

(  X2/F$  )  -  (X1  Y’)]  +  L(C-1)X‘  B1. 
~U4 


^  _  7.27x10-3 

‘"(TJ.-Tl,) 


1.585  1.585 

(TIz)  -  (T?i) 

1.585 


If  Ch  |  <  Cc  2 
go  to  (102) 

* 


Q?' 


K 1  •  Ak? 


If  Ch  |  >  Cc  2 
go  to  (105) 


If  Ch  |  3  Cc2  ilO' 

go  to  (108) 


I 


I 


(98) 


(99) 


(100) 


(101)* 


1 

1  -  X 


(103)* 


«  -  «  i 


TEST  7-A 

If  (103)  5  1.0,  stop  &  readout  message . "REDUCE  Nrui" 

"  "  <  ",  continue. 


Cht  /  \ 

1-  «  i  Cc2  j 
1 


1  - 

Then  go  to  LG-14 


lo9e 


(104)* 


If  Chi  >  Cc2 


Qt  (105)* 

3600  Cc2 


«  - 


* 

c  i 


J _ 

1  -  X 


TEST  7-B 

If  (106)  >  1.0,  stop  &  reodout  meisoge . "REDUCE  Ntui" 

"  "  <  ",  continue. 


(106)* 


Ntu 


•log  e 


|Cc2/ 

1  - « 

l  /chi 

L  i  ■ 

•  < 

1. 


(107)* 


Then  go  to  LG  “14 


np  calculated  = 


144  Ax  tot,  hot  side 
Axp 


013) 


Note . If  a  fraction  results,  go  to  next 

higher  whole  number. 


TEST  8-A 
(113)  must  =  (4) 

!f  (113)  >  (4),  increase  (4)  and  iterate  from  (4). 

"  "  <  ",  reduce  "  . . . 


Then  go  to  (120) 


If  Ch  t  >  Ccz 
- - - 


Ax  tot.  hot  side 


3600  Ntu  Ccz 

U 


(1)4)* 


Axp  3  Ax  i  •  Afi 


(115)* 


i  i  .  j  _  144  Ax  tot.  hot  side 

np  colcuioted  3  —  ■  ■  -  ■■ 

A  vn 


Note. ...  If  o  fraction  results,  go  to  next 
higher  whole  number. 

TEST  8-B 

(116)  must  3  (4)  +J 

If  (116)  >  (4),  increase  (4)  ond  iterate  from  (4). 

"  "  <  ",  reduce  "  "  "  "  " 


(116) 


Then  go  to  (120) 


74 


75 


core 


weight 


header  weight  -  .196  [(XY)  -  Afl  (Ro  +  1)  +  **  = 


tbs 


total  weight  =  (123)  +  (124)  =  lbs 


(123)* 


(124) * 

(125) * 


1 


r  h  /Ax?/  \  -i 

1  + 

i-2  V  / Nh  •  Y  i  )  (Y‘i)Z 

3  np2.,{x  'tp  .93041  [7-  (dtv^Zf]  }J 

TEST  9 

if  (126)  <  .40,  readout  message  . "INCREASE  (IP-24)." 

"  "  >  ,60,  readout  message  ...  ....  “REDUCE  "  ." 

Do  not  stop  mochlne  on  TEST  9. 


(126)* 


Av 


Ax  tot.  hot  side 
(  X  •  y  .  L  V 

V1725 — ; 


ft2 /ft 


FINAL  TEST 

If  (35)  ^  540,  stop  machine. 

"  "  ^  ",  coll  HX  (J)  J  =  2.0  and  continue. 


(127)* 


<  u 


FLOW  DIAGRAM  FOR  HX-1 


-  assume 
=  compute 
LG  =  logic  gate 
T  =  test 

C 

A 

A 

A 

LG-1 

C 


(71)-V(72) 


F 

ki. 


Call  HX  (J)  and  continue 

i  =2.0 


80 


APPENDIX  ill 


HX-2  AND  TURBINE  1  PROGRAM . CALL  HX  (J)  .  *2.0* 

J 


Inputs: 


Call  numerical  values  from  APPENDIX  I,  Section  III. 
Also  col!  last  result  for  TI3,  P13,  T17  ond  P 17  from 
output  af  HX-1 . 


Initlol  numerical  assumtions: 


Eqa. 

Initio! 

No: 

volue: 

(4) 

30.0 

(5) 

T13  +  5 

(6) 

1.02  P13 

(13) 

2.0 

(41) 

T17  +  5 

(42) 

.98  P17 

Notes: 

Readout  lost  result  of  all  equotlans  marked  with  o  star  11  £  11 
Do  not  stop  machine  far  TEST  9. ' 


81 


HOT  SIDE: 


s  =  "A/ 906894 

v  <r 


_  4cr 


Axi  =  (n  Jfd h  tp)  +  2  (J-cr) 


r>p  *  assume 

(Will  be  oltered  only  by  TEST  8) 


T12  ■  assume 

(Will  be  altered  only  by  TEST  3) 


P12  *»  ossume 

(Will  be  oltered  only  by  TEST  2) 


(11) 


ahi 


144  W1  /y-12 
VJ2 


(12)* 


r  Af  i  calculated  =  Nh  (X' Yt1)  (17) 

TEST  ! 

(17)  must  ■  (12)  *.001 

l  If  (17)  >  (12),  reduce  (13)  and  iterate  from  (13). 

"  “  <  "  ,  increase . . 


83 


Ra 

2 


08)* 


V'2  =  Y'l 


■  Y’l  •  Ra 


(19)* 


u  ■ 


III 

24 


(20) 


(21) 


v 


h  =•  Y'S 

12  •  Ra 


(22) 


AP'i 


370 x  10*6 

V  12 


(23) 


AP!  *  np  •  AP't 


(24)* 


P12  calculated  =  P13  +  APi 


(25) 


TEST  2 

(25)  must  -  (6)  ±.001 

If  (25)  >  (6),  increase  (6)  and  iterate  from  (6). 
"  "  <  "  ,  reduce  . .  "  "• 


Call  cp  subroutine  and,  with 

T12  and  PI 2  . get 

- - 

cP12  =  V 


(26) 


84 


TEST  3 

(31)  mint  ■  (5)  ±.001 

If  (3i)  >  *5)/  InertoM  (5)  and  It«rat«  from  (5). 
"  "  <  '*  ,  raduct  M  "  . . 


85 


AP2. « °~6  V ,  /_mm 

&  io  y  NRe2 


(49) 


A P?  *  np  •  AP*2 


(50)* 


P18  calculated  *  P 17  -  AP2  (51) 


TEST  4 

(51)  must  -  '  t.001 

If  (51)  >  (42/,  .screase  (42)  and  iterate  from  (42). 
"  "  <  "  ,  reduce . .  "  "  . 


TESTS 

If  (51)  10,  stop  &  readout  message . "INCREASE  Ra" 

"  "  >  "  ,  continue. 


.  _  8.55497  x  10-4 
^  (T18  -  T17) 


(52) 


Km2 


57.79  x  10“3  1  [{.00355  T18)1 ,642-  (.00355  T17)1 ,642]( 

[.  00355  (T 1  8  -  T17JU  ]  V.642 


(53) 


Tern 


(■ 


Km  2 


57.79  x  10 


-) 


1/.642 


.00355 


(54) 


Pcm  » 


(P17  +  P18) 
2 


(55) 


88 


Then  go  to  (70) 


TEST  6-6 

(64)  mutt  *  (41)  *.001 

If  (64)  >  (41),  increoie  (41)  and  Iterate  from  (41). 
"  "  <  "  ,  reduce  H  . , 


(62)* 


(63)* 


(64) 


(65)* 
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T18  calculated  =  T12  -  (T13  -  T17)  (67) 

TEST  6-C 

(67)  must  81  (41)  ±.001 

If  (67)  >  (41),  increase  (41)  and  iterate  from  (41). 

"  "  <  "  ,  reduce  . . 


Then  go  to  (70) 


91 


Thort  go  to  (78) 


%  & 


Kp  = 


07  [(J  T*)3'47-(.l  Tcm)3-4?] 

(86) 


Then  go  to  (92) 


Kp  = 


^'2  -  Terri)* 


(.1  T*) 


2.708 


t  L-' rcm)  2,70e^ 

2.708 


j  + 9.551  [(,  I  T*2)2  -  (j  Tcm)2] 


(87) 


Then  go  to  (92) 


Kp  “  111.74  =  conttonf 


Then  go  to  (92) 


Then  go  to  (92) 


(88) 


(89) 


(90) 
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Kp  *  92.25  *  constant 


(91) 


Then  go  to  (92) 


K  =  Kp  .93061 

4-5 


s  -  dh 

s  -  <4b/2) 
— 


(92)* 


U  „  — -1—  1-  (93)* 

x+k  +  k  k  + 1 r 

1-2  2-3  3-4  4-5  3-6 


ns  *  np  +  1 

(94)' 

ns  ‘  ts 

12 

(95) 

m-  (T,32+T'7) 

(96) 

n2mrnim 

(97) 

r(X2As)  -  (X*  V*71  +  Be-  1)  X'  B7] 

(98) 

144 

kI  r  7  27  *  10~3 
*(TZz  -TZl  ) 


1.585  1.585" 

g?z )  -(tzi) 

1.585 


(99) 


Q1 


K?  •  Ak2 


(100)* 
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(102)* 


TEST  7- A 


'"REDUCE  BORDER  DIMENSIONS 

If  (102)  >  1.0,  stop  &  reodout  message . OR  INCREASE  fs" 

If  (102)  <  1.0,  continue. 


Then  go  to  LG  -  14 


98 


<■  6'(-tV) 


(105)* 


TEST  7-B 


"REDUCE  BORDER  DIMENSIONS 

If  (105)  ^1.0,  stop  &  readout  message . OR  INCREASE  ts" 


If  (105)0.0,  continue. 


Ntu 


ioge 


1- 


6  (Cc2/Chi). 


T-_  «c2/£hl  ) 


(106)* 


Then  go  ta  LG  -  14 


-6 

If  Chi  s  Cc2  ±  10 


.  Ql 
*  36K)  Chi 


(107)* 


(108)* 


TEST  7-C 

If  (108)  ^1.0,  stop  &  readout  message 
If  (108)  <1.0,  continue. 


"REDUCE  BORDER  DIMENSIONS 
OR  INCREASE  ts" 


Ntu  = 


^  e 


(io?)* 


99 


Then  to  LG  -14 


If  Chi  <  Ccz 
go  to  (110) 


If  Chi  >Cc2 
go  to  (1 13) 


If  Chi  s  Cc2  ±  10 
go  to  (116) 


If  Chi  <  Cc2 


Ax  tot.  hot  side  = 


3600  Ntu  Chi 


Axp  =  Axi  ♦  Afi 


n/5  calculoted  = 


144  Ax  tot.  hot  side 


Note 


If  o  fraction  results,  go  to  next 
higher  whale  number. 


TEST  8-A 


(112)  must  =  (4) 


I  ( 1 » 2) >  (4),  increase  (4)  and  iterate  from  (4). 
If  (112)<  (4),  reduce  (4)  and  iterate  from  (4), 


Then  go  to  (119) 


If  Chi  >Cc2 


Ax  tot.  hoi  side 


3600  Ntu  Cc2 
U 


(1 13)* 


Axp  =  Axi  •  Afi 


(114)* 


np  colculated  - 


144  A.t  tot.  hoi  side 
Axp 


(MS) 


Note . If  o  fraction  results,  go  to  next 

higher  whole  numbei . 


TEST  8-B 

(115)  must  *  (4)  ^ 

If  (115)  >  (4),  increase  (4)  ond  iterote  from  (4). 
If  (115)  <  (4),  reduce  (4)  ond  iterate  from  (4). 

Then  go  to  (119) 


If  Chi  *Cc2±  10“6  | 

(116)* 
(117)* 


Ax  tot.  hot  side 


3600  Ntu  Chi 
— 0 - 


Axp  *  Axi 


Afi 


101 


144  Ax  tot.  hot  side 
up  calculated  -  7\xi  ” 


(118) 


Note  .....  If  a  froction  results,  go  tc  next 
higher  whale  number. 


TEST  8-C 
(118)  must  =  (4) 


If  (118)  >  (4),  increase  (4)  ond  iterate  from  (4). 

If  (1 18)  c  (4),  reduce  (4)  and  iterate  from  (4). 

Then  go  to  (119) 

width  X  =  equation  (13)  _  jnc[ies 
after  closure 


height  Y  =  pp  =  inches 

length  L  =  j^(np  •  tp)  +  (ns  •  ts)]  =  inches 


care 


ight  =  .098  np  fp<j[(XY)  -  (X'Y')j  +  [x'B'  (C-l)]  +  [Afi  (R»+l  )  (1-<T 
.078  ns  t*  j[(XY)  -  (X’Y')]  +  [x'B1  (C-l)] |  =  lbs 


header  weight  =  .196 


XY 


(XY)  -  Afi  (Ra+1)  +  -g- 


=  !bs 


total  weight  -  (122)  +(123)  =  lbs 


(119) * 

(120) * 
(121)* 


(122)* 

(123) * 

(124) * 
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If  (!25)  <  .40,  read  out  message . "INCREASE  ( 1  P—40) 

!f  {125)  >  .60,  readout  message . "REDUCE  (IP-40)." 


Da  not  stop  machine  an  TEST  9 


Ax  tat.  hat  side 
/X  • Y  •  L  \ 


/  X  •  Y  •  L  \ 

\  Tt5§  J 

A Hti  =  cp  T12  ??n-(P17/P12){  *~])/  7  ='  BTU/lb 


Turbine  output  =  1054.54  Wti  >AHtl  ~ 


FINAL  TEST 

If  (41)  2S  540,  stop  machine. 

If  (41)  <  540,  call  HX{J)  js3  and  continue. 


APPENDIX  IV 


HX-3  PROGRA 


CALL  HX(J)j*  3.0: 


Inputs: 


Call  numerical  values  from  APPENDIX  I,  Section  IV. 
Also  call  last  resut  far  T12,  PI 2,  T18  and  PI 8  from 
output  af  HX-2. 


Initial  numerical  assumptions: 


Equc. 

Initial 

Na: 

Value: 

(4) 

200.0 

(5) 

T12  +75.0 

(6) 

1.02  PI  2 

03) 

2.0 

(35) 

T18  +  75 

(36) 

.98  P18 

Nates: 


Readout  last  result  of  all  equations  marked  with  a  star 
Da  nat  stop  machine  at  TEST  9. 


HOT  SIDE: 


n  - 


4  ff 

TTdh7 


Ax  I  =  (n  Tf  dh  tp  )  +  2  (1  -  C  ) 


np  =  assume 

(Will  be  altered  only  by  TEST  8) 


0)* 

(2)* 

(3)* 

(4>* 


108 


Tt  1  =  assume 


(5)* 


'Will  be  altered  only  by  TEST  6) 
P1 1  =  assume 


(Will  be  altered  only  by  TEST  2) 


_  144  NRei  Vn  Uw 
VT1  ^ - 

1*4  W2  Yu 
ahr  «■ - yn - - 

Af,  -  ±- 

cr 

X  *  assume 

(Will  be  altered  only  be  TEST  1) 
X'  =  X  -  (2  Bx) 

v/,  * 

Y  =  "P*  “  (2  By) 


(6)* 


(8) 

(?) 

(10)* 

(11)* 

(12)* 

(13)* 


04)* 

(15)* 


109 


(16)* 


Y'i 


V  -  L(C-1)  0'!  -  Nh 

Fc^iT 


Af.  calculated  Mil  (X1  Y'i  ' 


(17) 


TEST  1 

(17)  must  -  (12)  i. 001 

If  (17)  >(12),  reduce  (13)  and  iterate  from  (13). 

If  (17)<(12),  incraase  (13)  and  iterate  from  (13), 


Y3  =  Y)  •  Ra 
Y‘l 

h  =  ~TT 


AP'i 


h=  ly 

,  Y^ 

T5”rRa' 

370  x  10"6  Vi  i  2  n  /  (tp/dh) 
/i>  i  j  V  NRei 


aPi  =  np  ♦  AP't 


P1 1  calculated  =  PI  2  +  APl 


(18)* 

(19) * 

(20) 

(21) 

(22) 

(23) 

(24) * 

(25) 


110 


TEST  2 


(25)  must  =  (6)  ±  .001 

If  (25)  >  (6),  increase  (6)  and  iterate  from  (6). 
If  (25) <  (6),  reduce  (6)  and  iterate  from  (6). 


imi 


=  8,55497  x  10*^  \  JT11)1,643  .(T12)1'643 


7TT1"  -  Tl2) 


1.643 


if*,  =  57.79  x  JO*3 

■puss  (ill . 


(.00355 Til)1  ’542,  (.0,3355  T|2)  1.642 
1.642 


J&L. 


1/.  642 


Thm  = 


,57.79  x  IQ"3 
.00355 


(26) 

(27) 


(28) 


Phm 


(Pll  »  P-l 2) 


(29) 


Call  cp  subroutine  and,  with 

Thm  and  Phm . get  j 

- 1 

cphm  =V 

(30) 

Chi  =  cphm  •  W2 

(3D* 

NPri - gi??  '  /■"" 

Kmi 

(32)* 

111 


Nnui  “ 


COLD  SIDE: 


3.  66  -*•  < 


,104 _ 

Cfp/dV»1 
.n$*»  “nET- 


1+ 


.016 

"  (/p MY  _ 

NRei  •  NPri 


(33)* 


h  - 

I  -  2 


12  NNui  •  Kmi- 
dh 


(34V 


T19  =  assume 

(Will  be  altered  only  by  TEST  5) 


(35V 


PP  =  assume 

(W5,l  be  altered  only  by  TEST  3) 


(36V 


*'9  =  2.863092 


(38) 


V,9=Ui^ 

(40) 

x<i9  =  2.37888  x  I0"7  (T19)’643 

(41) 

VI 9  •  dh 

12  Vis 

(42) 

.  r.  370  x  10  4  V^9  /(tp/dh) 

=  Vi.  V  nS.2 

(43) 

A  P2  =  np  *  A  P*2 

(44) 

P19  calculated  =  P18  -  AP2 

(45) 

TEST  3 

(45)  mutt  =  (36)  ±.001 

If  (45)  >(35),  Irwreaw  (36)  and  Iterate  from  (36) 

If  (45)  <(36),  reduce  (36)  and  Iterate  from  (36). 

TEST  4 

"REDUCE  NRei  OR 

If  (45)  <10,  stop  &  readout  montage . INCREASE  Ra" 

If  (45)- 10,  continue. 


113 


Km  2 


8.55497  x  i_-  ' 

“Tfrr-nrsr  " 


JiH/i 1  • ' i3  -  (Tieo1  -443. 1 

j  r  :m 

l 


(46) 


57.79  x  10 


-3 


[700155  (Tl9  -  T18)‘ 


T  ^ 


00355 


l.c',2  ,  Ttn.!  .642 

-  (.00355  Tl 8) 

rui - 


(47) 


Km? 


vl/.  642 


Tcm  = 


,57.79  x  10 
.00355 


-3 


(48' 


„  (PI  8  +  PI  9) 

Pcm "  2 

(49) 

Call  cp  subroutine  and,  with 
Tcm  and  Pern  ....  get 

cpcm=  ^ ^ 

(50) 

Cc2  =  epem  •  W2 

(S’) 

NPr,  -  epem  •  AM  mz 

Kmz 

(52) 

NNu2  =  3.66  W 


1 

l 


.104 

WM 


[NRez 


NPrzJ/ 


.016 

w 


|NRe2  •  NPrj 


.8 


) 


(53)* 
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f 


If  Chi  <  Cc2 
go  to  (55) 


I 


12  NNu2  •  Km2 
& - * - 


-6 


go  to  (61) 

I 


go  to  (67) 


I 


(54)’ 


If  Chi  <  Cc2 


V  1  =  (T12  -  T18)e 


Ntui  (Chi/Cc2)J 


(55) 


a  w  Vi  -  (T12  -  Tl  8) 

Ax - r^rr — 

ten  1 


(56) 


AY  =(^r) ' A  x 


(57) 


T19  ealculoted  =  T18  +  AX 


(58) 


TEST  5 -A 

(58)  must  =  (35)  ±  .001 


If  (58)  >(35) ,  increu:*  (35)  ond  Iterate  from  (35). 

If  (58) ^SS),  reduce  (35)  and  iterote  from  (35). 

Til  calculated  =  Tl 2  +AY  (59) 
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TEST  6- A 


(59)  mu st  *  (5)  A. 00 1 

If  (59)  >  (5),  increase  (5)  and  iterate  from  (5). 
If  (59)<  (5),  reduce  (5)  and  iterate  from  (5). 


€/= 


-Ntui  £l-  (Chi  /fcczjj 

,  <vcb.  \  -N,ui  0  -  tch.  /c«)]l 


\Za. 


) 


I 


Then  go  to  (71) 


If  Chi  >Cc2 


VI  -  (TI2-  T18)e 


Ntui[l-(Cc?/Chi)] 


AX  - 


Vi-  (T12  -  T18) 


AY 


-(&■)• 


AX 


T19  calculated  =  TT 8  +  AX 


(60)* 


(61) 

(62) 

(63) 

(64) 


TEST  5-B 

(64)  must  =  (35)  ±.001 

If  (64)  >(35),  increase  (35)  and  iterate  from  (35). 
If  (64)<  (35),  reduce  (35)  and  iterate  from  (35). 
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Til  calculated  =  T12  +  AY 


(65) 


TEST  6-B 

(65)  must  -  (5)  ±.001 

If  (65)  >(5),  increase  (5)  an-*  iterate  fram  (5). 
If  (65)<(5),  reduce  (5)  and  iterate  fram  (5). 


-Ntui  [l-(Cc2/Cht  )] 

1-6 


Then  go  te  (71) 


If  Chi  =  Cc2  ±  10"6 


z  =  (T12  -  T18) 


T19  calculated  =  (T18  +  Z)  -  (T12  -  T18) 


TEST  5-C 

(68)  must  =  (35)  ±.001 

If  (68)  >(35),  increase  (35)  and  iterate  from  (35). 
If  (68) <(35),  reduce  (35)  and  iterate  fram  (35). 


Tl  1  calculated  =  TT8  +  Z 


(67) 

(68) 


m 


♦  x 
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TEST  6-C 


(69)  must  (S')  00 i 


(f  (69)  >  5'»  , 
If  (69 )<  (5), 


u,.jsa  o)  or  I  iiorote  from  (5). 
reduc-o  '5')  o-'H  .  urate  fro..'  (5) . 


Mtui 

ti-  T+"  Ntui 


(70)* 


Then  to  to  (71) 


(71) 


T’g  =  Thm  - 


(Thm  -  Tcm) 

\Y'i  +  Y  2  +  Y'3  / 


(72) 


r 


If  (IP-67)  CODE  =1.0 
ao  to  LG-5 

t 


r 


If  (IP-67)  CODE  =  2.0 
go  to  LG-6 


LG-5 


If  0  <  Trim  <  1 1 3 
go  to  (73) 


|f  nxTlvn^540 
go  to  (74) 


I 


If  540<Thm<oo 

go  to  (75) 
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Then  go  to  (79) 


Kp  =  111.74  =  constant'  (75) 

Then  go  to  (79) 


Then  go  to  (79) 
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K.p  '  92.25  ~  constant 


(78) 


Then  90  to  (79) 


If  0<  T'i  <  1 13  If  1 13  <  T'i  <  540  If  540  <  T‘i  <  00 

go  to  (80)  go  to  (81)  go  to  (82) 


*  *  I 


(81) 


Then  go  to  (86) 
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Kp  =  1 11 . 74  =  constant 


If  0  <  T'l  <210 
go  to  (83) 

t 


Kp= 


1 

(T’l  -Tz 7 


Then  go  to  (86) 


(T'l)J-(T'2)J 


(T’l)2’'4  -  (Ts)2'16 


2.16 


Then  go  to  (86) 


Kp  -  86.0  + 


6.25 


r’i+T‘2)/2  -210 

550 


Then  go  to  (86) 


Kp  =  92.25  =  constant 


Then  go  to  (86) 


K  -  Kp  .93061 

4-  f 


s  -  dh 
s  -  (dh7Sy 


U  = 


!  2 


T  7  ^  .  ]z  +  __L  \ 
~r  T-  +  k  nr ) 

2-3  3-4  4-5  5-6  / 


ns  -  np 


4-  1 


ns  •  ts 
le  =  “T r~ 


(T12  +T18) 
Tli  =  - - 5 - - 


(Til  +T19) 
- 1  " 


Akl  = 


[(x2/f.)  -  W  V)1  *  l(C-t)  X'BL 

— - — nr 


,-3 


K1 


7.27  x  10 

7T72  -  T  i  r*T 


i  535  1.585 

(TH  )  -  (Till _ 


1  .585 


(93)  ^ 


(94)  - 


(95) 1 


(96) 


(97) 


(98) 


(99) 


(100) 


124 


f 


Q\ 


Kl  •  All 

le 


(101)* 


-6 

10 
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Ql 

i-.  DO  .  , 


n  02)* 


..  *.  I..LJ 

\  1  -  \  I 


Tt ST  7  -  A 

If  ^  1 03)  >  1  .0,  stop  &  readout  message . "REDUCE  Ni'ui" 

"  <  "  ,  continue. 


(106)* 


TEST  7-8 

If  M06)  >  1.0,  stop  &  readout  message . "REDULfc  Ntui" 

"  "  <  ",  continue . 
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TEST  7  -  C 
If  (109)  >  1 

II  ll  y 


r~ 

If  Chl  <  Cc2 
go  to  (1 1 1) 

I 


Then  go  to  LG  -  14 


If  Cht  =  Ccz  ±  10“6 


\  =  Q£ 

“3600  Chl  (lu8)* 


t  =  6i 


(109)* 


.0,  stop  &  readout  message . "REDUCE  Ntui" 

11  ,  continue. 


N?  j  =  € 

~\  -  €  010)* 


Then  go  to  LG  -  14 


\ 


I 
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If  Chi  <  C  2 


As  Uit.  hot  side  -  3600  Ntu  _ ChJ _ 

U  (Ml)* 


Axp  _  Ax  I  •  Afl  (1 12)* 


np  calculated  =  144  Ax  tot,  hot  side  _ 

Ax  p  (113) 


Note . If  a  fraction  results,  go  to  next 

higher  whole  number. 


TEST  8 -A 

(113)  must  -  (4)  +_  q 

If  ( 1  l  3 )  >  (4),  increases' (4)  3nd  iterate  from  (4), 
"  "  <  reduce  "  "  "  ", 


Then  go  to  (120) 


If  Chi  >  Cc2 


As  !-  t.  hot  side  3600  Ntu  _ Cc2 

U 


(114)* 


Axp  =  Ax  I  •  Afl 


(115)* 


np  calculated  i44  Ax  tot,  h.ot  side 

Axp 


(116) 
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Note 


If  a  fraction  results,  go  to  next 
higher  whole  number. 


TEST  8  -  B 

+  1 

(116)  must  =  (4)  -  0 

If  (116)>  (4)f  increase  (^)  and  iterate  from  (4). 
"  "  <  ",  reduce  "  »'  " 


Then  go  to  ( !  20) 


If  Cht  =  Co2  +  10  "6 


Ax  tot.  hot  aide  =  -16.9P  Ntu... . Chi 

U 


(H7) 


Axp  =  Ax  I  «  Af  I 


(118) 


np  calculated  =  1  14  Ax  tot,  hot  side 

Axp  (119) 


Note . If  a  fraction  results,  go  to  next 

higher  whol  i  number. 


TEST  8  -  C 

+  1 

(119)  must  =  (4)  -  0 

If  ( 1 19)  >  (4),  increase  (4)  and  iterate  from  (4). 
"  "  <  ",  reduce  "  "  ''  "  . 


Then  go  to  ( 1 20)  ■ 
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width  X  cquotioi'i  si  ->)  =  inche1. 
ufte.  closure 


(120)* 


height  . 


.  = inches 


021)* 


Fs 


core  length  L  =  jjnp  •  tp)  4  (ns  •  ts)j  =  inches 


(122)* 


?  weight  -  .098  np  tp  jjXY)  -  (X1  Y1)]  +  [X*  B'  (C  -  lj]  4  [Af  i  (Ra  4  1)(1  -  aj  'j 

.078  ns  ts|[(XY)  -  (X'  Y'f]  +  [X'  B'  (c  -  l]]J  =  lbs 


xy| 

(XY)-Afi  (Ra  +  D+  TJ  * 


header  weight  =  .  196  |(XY)  -  Af  i  (Ra 

total  weight  =  (123)  +  (124)  =  lbs 


lbs 


(123) * 

(124) * 

(125) * 


nf  - 


TEST  9 

If  (126)  <  .40,  readout  message . "INCREASE  (IP-60)." 

"  "  >  .60,  "  "  . "REDUCE  (IP-60)." 

Da  nat  stop  mochine  on  TEST  9. 


(126)* 


Av  -  _Ax  tot,  hat  side  rft2/f*3 
'  Y  *  L  \ 

\  1728  J 


(127)* 
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FINAL  TEST 

If  (35)  >  540,  stop  machine. 

"  "  <  "  ,  call  HX  (J)  and  runtime. 

j  =  4 


The  FLOW  DIAGRAM  for  HX-3  is  exactly  similar 
to  that  of  HX-1. 
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APPENDIX  V 


HX-4  AND  TURBINE  2  PROGRAM . CALL  HX(J)  .  ^  : 

1 :  4. 0 


Inputs: 


Call  numerical  valves  from  APPENDIX  1,  Section  V. 
Also  c<  it  last  result  far  T1 1,  PIT,  T19  and  P19  fram 
output  af  HX-3- 


Initial  numerical  assumptions: 


Equa. 

Initial 

Na. 

Value: 

(4) 

30.0 

(5) 

TIT  +10 

(6) 

1.02  Pll 

(13) 

2.0 

(41) 

T19  +10 

(42) 

.98  PI 9 

Notes : 


Readout  last  result  af  oil  equations  morked  with  o  star 
Do  :at  stop  machine  at  TEST  9. 
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HOT  SIDE: 


Axi  *  (nTfdh  tp)  +  2(1-  o') 

np  *  assume 

(Will  b«  altered  only  by  TEST  8) 
110  x  assume 

(Will  be  altered  only  by  TEST  3) 

P10  *  assume 

QWHI  be  altered  only  by  TEST  2!) 


»,  12  NRe!  yVio  *aio 

v,°  =  -  dh 

(10)* 

144  W2  v  10 
<,hl”-  ViO 

<nr 

Afi  -  ahi 

(12)* 

or 

X  -  assume 

03)* 

(Will  be  altered  only  by  TEST  1) 

X1  =  X  -  (2  Bx) 

(14)* 

v>  =  2L  -  (2  By) 

(15)* 

Y*  -  f  (C-l)  B'l  -  Nh 

Yl - 2(5-1) 

(16)* 

Afl  calculated  -  Nh  (X*  Vl  ) 

(17) 

TEST  1 

(17)  must  =  (12)  ±  .001 

If  {17)  >(12),  reduce  (13)  and  iterate  from  (13). 

»  »  <  “  ,  increase  "  "  "  ''  "• 


Yfc 

•  Vi  4 

(18)* 

Y*3 

-  Vi  •  Fa 

(19)* 
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Y‘3 
12*  Ra 


A  P'l  * 


370  x  10'6  V2to 


TEST  2 


API  =  np  •  AP'l 


P10  calculated  =  PI  1  +  APi 


(25)  must  -  (6)  *  .001 

If  (25)  >  (6)  increase  (6)  and  iterate  from  (6). 


h  h  <  »#  roduce 


II  M  II  M  It 


Call  cp  subroutine  and,  with 
T10  and  P10 . get 


cpio  = 


Call  cp  subroutine  and,  with 
T19  and  P19 . aet 


COLD  SIDE: 


T20  c  :sume 

(<D* 

(Will  be  altered  only  by  TEST 6) 


P20  -  assume 


(42)* 

(Will  be  altered  only  by  TEST  4) 
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nu  iwwi  jMuyTy 


JIM 


AP'2  - 


370  x  10~6  V220 
AT  20 


ff^/dh) 

NRez 


AP2  =  np  •  AP'2 
P20  calculated  =  PI 9  -  AP2 

TEST  4 

(51)  must  =  (42)  ±  .001 

If  (51)  >(42),  increase  (42)  and  iterate  from  (42). 

"  "  <  "  ,  reduce  "  "  "  "  “  . 


TEST  5 

if  (51  )<  10,  stop  and  readout  message 


11  II 


^  ",  continue. 


JUL  m2  = 


8.55497  x  10 
(T20  -  T19) 


-4 


"INCREASE  Ra" 


(T20)1,643  -  (T19)1,643 
1.643 


Km2  = 


57.79  xlO-3  J(.00355T20)1-642-<-00355Tl9> 


1.642 


[.00355  (T20  -  T19)] 


Tcm  = 


/  Kmt  \ 
\57.79x  10“3/ 


1.642 

1/.642 


705333“ 


(49) 

(50) * 

(51) 


(52) 


(53) 


(54) 


(55) 


NNu2  = 


Call  cp  subroutine  and,  with 
Tcm  and  Pcm  ........ .get 


cpcm  = 


Cc2  =  cpcm  •  W3 


NP,2  -  -CE”  •  .-fJ-? 
Km  2 


/ 


.104 


l.66+< 


1  + 


NRe2  •  Npr2j 

/ 

\ 

7  .616 

lr  dp/dh)  i 

i 

\[Nft#2  •  NPr2 

h 

5-6  = 


12  NNu2  •  Km2 


^gggBT-LLWIkJl  1 U _ 'UjipU  i-1  j  llUyULW. .  ■  J  ■  »-JlJLLH  Hij 


If  Chi  <  Cc2 


T20  calculated  =  T19  + 


(T10  -  Tl  1) 
(Cc2/Chi  ) 


(61) 


TEST  6  -  A 

(61)  must1  =  (41)  ±  .001 

If  (61)  >  (41),  increase  (41)  and  iterate  from  (41). 
"  "  <  ",  reduce  .  "  ", 


Then  go  to  (70) 


If  Chi  >  Cc2 


T20  collated  T19  +  ~  T11) 

l  (Ch'/Cc2)  J 

TEST  6  -  B 

(64)  must  =  (41)  -fc  .001 

If  (64)  >  (41),  increase  (41)  and  iterate  from  (41). 

"  "  <  "  ,  reduce  "  "  "  "  " , 


Then  go  to  (70) 
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(67) 


TEST  6  -  C 

(67)  must  =  (41)  ±  .001 

lf  (67)  >  (40,  increase  (41)  and  iterate  from  (41). 
"  "  <  ",  reduce  .  « 


»  143  V 


Then  go  to  (78) 


Then  go  to  (85) 


r 

If0<  T*i 
go  to  (82) 


Kp  =  111.74  **  constant 
Than  go  to  (85) 


go  to  (83) 


1 
If  540  <  r»  <o° 
ga  to  (84) 


If  0<  T'2<  113 
go  to  (86) 


If  113  <  T*2  <  540 
go  to  (87) 


If  540  <  T'?<  oo 
go  to  (88) 


„  1  J  49 

Kp  *  r.l  (T'2  -  Tcm)1  \T 


Lit’*)2  -  (.lTen>)Jj-  J^7  |(.  1  T'a)3'47  -  (.Item)3-' 


Then  go  to  (92) 


HI  (T'2  -  Tcm) 


,  1Ti  ,2.708  .  . T  ,2.708  r  o 

illlil — - +9.551  f(.lT'2)2-(.lTcm) 


Then  go  to  (92) 


Kp  =  111.74  =  constant 


Then  go  to  (92) 


If  0  <  Tie  <  210  If  210  <  Tfe  <  540  If  540  <  Tie  <00 

go  to  (89)  go  to  (90)  go  to  (91) 


Then  go  to  (92) 


Kp  =  92.25  =  constant 
Then  go  to  (92) 


WJUPPL* 


If  Chi  <  Cc2j 


v  -  Qj 
*  3600  Chi 


(101)* 


1  (102)* 

r~nr 


TEST  7  -  A 

If  (102)  >  1 .0/  stop  and  readout  message 


"REDUCE  BORDER  DIMS.  OR 
INCREASE  ts" 


"  "  <  "  ,  continue. 


Ntu  - 


Then  go  to  LG  -  14 


If  Chi  >  Cc2 


(103)* 


\  _  u  t 

3666  Cc2 


(104)* 


(105)* 
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TEST  7  -  B 


If  (105)  >  1.0,  stop  and  readout  message 
"  "  <  "  ,  continue 


"REDUCE  BORDER  DIMS.  OR 
INCREASE  ts" 


log 

Ntu  =  6 

1  -  *  (c7cm) 

.  1  -  « 

Cc2  / 

1  -  /Chi  ' 

(106)* 


Then  go  to  LG  -  14 


If  Chi  25  Cc2  ±  10"6 


X  "  3600  Chi  (107)* 


« 


=  «1 


1 

1  -  X 

i 


(108)* 


TEST  7  -  C 

If  (108)  >  1 .0,  stop  ond  readout  message 
"  "  <  "  ,  continue 


••REDUCE  BORDER  DIMS.  OR 
INCREASE  ts" 


Ntu  -  -r— - 

I  -  « 


(109)* 


Then  go  to  LG  -  14 
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If  Chi  <  Cc2 
go  to  (110) 


* 


If  Chi  >  Cc2 
go  to  (113) 


* 


If  Chi  s  Cc2 
go  to  (116) 


± 


10 


If  Chi  <  Cc2 


Ax  tot.  hot  lido  “ 


3600  Ntu 
U 


Chi 

( 


Axp  =  Axl 


Afl 


(' 


np  colculoted  = 


144  Ax  tot,  hot  tide 
Axp 


< 


Note . ..  If  o  froction  results,  go  to  next  higher  whole  number 


TEST  8  -  A 
(112)  mull  =  (4)  ^ 

if  (112)  >  (4),  increase  (4)  ond  iterote  from  (4). 
"  "  <  "  ,  reduce  "  “  "  " 


Then  go  to  (1 19) 

if  Chi  >  Cc2 


n 


Ax  tot.  hot  side 


3600  Ntu  *  Cc2 


Axp  ~  Axi  *  Af l 

(114) 

np  calculated  =  144  Ax  tot.  hot  side 

(115) 

Axp 

Note  ....  If  g  fraction  results,  go  to  next  higher  whole  number. 


TEST  8  -  B 
(115)  must  =  (4)  *q 

If  (115)  >  (4),  increase  (4)  and  iterate  from  (4). 
"  "  <  "  ,  reduce  "  "  . . 


Then  go  to  (119) 


If  Oil  ~  Ccz  ±  10 


-6 


Ax  tot.  hot  side  » '^22±!iii5!LL 
Axp  =  Axi  •  Afi 


(116)* 

(117)* 


np  calculated 


144  Ax  tot,  hot  side 
Axp 


(118) 


Note  ....  If  a  fraction  results,  go  to  next  higher  whole  number. 


TEST  8  -  C 
(118)im»t  =  <4>  tj 

If  (118)  >  (4),  increase  (4)  and  iterate  from  (4). 
"  11  <  " ,  reduce  "  "  ”  "  "  . 


Then  go  to  (119) 
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| 


s 


.1 

i 


I 

j 

i 

I 

% 

i 


/ 


width  X  -  equotion  (13)  =:  inches 
of  ter  closure 


(119)’ 


height  Y  =  -  inches 


(120)* 


core  length  L  =  £(np  *  tp)  +  (ns  •  ts)j  =  inches 


(121)* 


core 


weight  =  .098  np  tp  j[(XY)  -  (X'Y')]  +  [x'B1  (C-l)]+[Afi  (Ro  +  1)  (1  -rf)]j+ 
.078  ns  ts  |[(XY)  -  (X'Y')]  +  [x'B*  (C -1 )  ] |  =  lbs 


(122)* 


heoder  weight  =  .196  £(XY)  -  Afl  (Ra  +  1)  +  ]  =  lbs 


(123)’ 


total  weight  =  (122)  +  (123)  =  lbs 


(124)’ 


1  + 

&  ( AxP/Nh  •  r, 

(V»  )2 

■ 

3  nP  2-3  |X*  tf  3061 

'  *  -  dh 

L7^)] 

i 

(125)* 


TEST  9 

lf(125)<.40,  readout  message . "INCREASE  (IP-76)." 

If  "  >  »  "  "  . "REDUCE  "  ." 


Do  not  stop  rnochine  on  TEST  9. 
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Ax  tot.  hot  side 
~  X  *  Y  •  L  1 
- ^ - 


ft2/  ft3 


A  Hta 


cp  T10  ijt 


»  BTU/lb 


Turbine  output  =  1054.54  Wt2  *  AHt*  -  Watts 


(126)* 


(127) * 

(128) * 


FINAL  TEST 

If  (41)  >  540,  stop  machine 

“  "  <  "  ,  coll  HX  (J)  ond  continue. 

j*  s.o 

The  FLOW  DIAGRAM  far  HX  •  4  ts  exactly  similar  to  that  of  HX  -  2, 
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APPENDIX  VI 


HX  -  5  PROGRAM  .  CALL  HX(J)Jl5  0  : 


Inputs: 

Call  numerical  values  from  APPENDIX  I,  Section  VI. 

Also  coll  last  result  far  T10,  P10,  T20  and  P20  from  output  af  HX-4 

Initial  numerical  assumptions: 


Equa.  Initial 

No  Vglue 

(4)  250.0 

(5)  T10  +  100 

(6)  1.02  P10 

(13)  2.0 

(35)  T20  +  100 

(36)  .98  P20 


Notes: 

Readout  last  result  of  all  equotions  marked  with  a  star  "  *  " 
Do  not  stop  machine  at  TEST  9. 
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HOT  SIDE: 


4  cr 

TTdh 


Axl  -  (n  IT  dh  tp)  +  2  (1  -  cr  ) 


np  -  assume 

{Will  be  altered  only  by  TEST  8) 


T9  -  assume 

{Will  be  altered  only  by  TEST  6) 


P9  =  assume 


(Will  be  altered  only  by  TEST  2) 


(9) 


TEST  1 

(17)  muit  =  (12)  ±.001 

If  (17)  >  (12) ,  reduce  (13)  and  Iterate  from  (13). 
"  "  <  ",  increase  "  "  . . 
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i 


ah2  ~  ohl  •  Ra 


(39)* 


V21  - 


144  WS  f\f  21 
aK£ 


(40)* 


-7  643 

AjL2\  *  2.37888  x  10  (T21)' 


NRe2  = 


V21  dh 


12^?2l  *^2l 


AP'2  = 


370  x  10 

fir 


P±4 


(tp/dh) 

NRe2 


AP2  =  np  •  A  F*2 

P21  calculated  =  P20  -  A  P2 


TEST  3 

(45)  must  =  (36)  ±  .001 

If  (45)  >  (36),  increose  (36)  and  iterate  from  (36). 
"  "  <  "  reduce  . . 


(41) 


(42)* 


(43) 

(44) * 

(45) 


TEST  4 

If  (45)  <  10,  stop  ond  readout  message 
"  "  .>  "  ,  continue. 


'.'REDUCE  NRei  OR 
INCREASE  Ra" 
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T2I  calculated  T20  +  A  ^ 


TEST  5  -  A 

(58)  must  (35}  )  .001 

If  (58)  >  (351,  Increase  (35)  and  iterate  fnr.m  (35). 

II  It  <  It  H  M  II  U  II 


,  reduce 


T9  calculated  T10  +  A  V 


TEST  6  -  A 

(59)  must  (5)  -fc  .001 

If  (59)  >  (5),  Increase  (5)  and  iterate  from  (5) 


"  "  <  ",  reduce  "  " 


II  tl  It 


6-i 


1  - 


1  -  ~ 


CTT 

Cc2 


f 1  ■  (ehl  /  P 

Fv3I[ 


e 


-  Ntui 


1  [ 


1  - 


Then  go  ta  (71) 


If  Chi  >  Cc2 


Vi"  (T10  -  T20)  e 


N.ui  [l  -  (Cc2  /  ai)] 


AX 


VI 


-  (TlO  -  T20) 

“TCR 


CZT 


-  1 


AV 


Chi 


Cc2 


AX 


p8) 


(59) 


(60)* 


(61) 


(62) 


(63) 
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i*#F,nt«p 


T21  calculated  -  T20  +  /VX 


TEST  5  -  B 

(64)  must  =  (35)  ±  .001 

If  (64)  >  (35),  Increase  (35)  and  iterate  from  (35). 
"  ii  ^  ii  t  recjUCe  "  "  "  "  " 


T9  calculated  =  T10  +  Ay 


TEST  6  -  B 

(65)  must  =  (5)  ±  .001 

If  (65)  >(5),  increase  (5)  and  iterate  from  (5). 

C  reduce  "  11  "  "  " 

11  m  ^  ii#  reauce 


6i 


-  Ntui 

I  -  e 


[l  -  (Cc  2 /Ch i )  ] 


Then  ga  ta  (71 ) 


If  Cht  m  Cc2  ±  10 


(TIP  -  T20) 


Ntui  1 

1+  Ntui 

T21  calculated  (T20  +  Z)  -  (T10  -  T20)  (68) 


TEST  5  -  C 

(68)  must  *  (35)  *  .001 

If  (68)  >(35),  increase  (35)  and  iterate  fram  (35) 

"  "  <  "  ,  reduce  . 

T9  calculated  ~  T20  +  Z 


TEST  6  -  C 

(69)  must  «  (5)  ±  .001 

If  (69)  >  (5)  ,  increase  (5)  and  iterate  from  (5) 

"  "  <  "  ,  reduce  . 

>•  _  Ntui 

1  +  Ntui 


Then  go  to  (71 ) 


T'l  *  Thm 


/  Vi  I 

(Y)  +  Y'2  +  Yl5 

’(  Y’l  +  Y'2 
\Y\  +  Yfe  +  Y's 


(69) 


(70)  * 


T*  =  Thm 


(72) 


If  0  <  Thm<  113 
go  to  (73) 


* 


If  113  <  Thm  <  540 
go  to  (74) 


* 


If  540  <  Thm  <  oo 
go  to  (75) 


* 


Then  go  to  (79) 


Then  go  to  (79) 


Kp  ~  tU. 74  '  constant 


(: 


♦ 

If  0  <  Thm  < 
go  to  (76) 


Kp  (Thm  - 


Kp  *=  86.0  + 


Then  go  to  (79) 


210 


If  540  <  Thm  <  oo 
go  to  (78) 


I 


* 


TTT 


2.765  ,2  .2] 

— = —  [(Thm)  -  (Tl  )  j  - 


(Thmr  -  (Tl  ) 

Then  go  to  (79) 


/tl  .2.16  ,2.16 

(Thm)  -  (Tl  ) 


2.16 


6.25 


^[(Thm  +  T1 


i  )/2]’  ^  210 

330 


>} 


Then  go  to  (79) 

KP  =  92.25  =  constant 

Then  go  to  (79) 


>•  • 


M  =  “P  -»304l  f-~y 


Then  go  to  (86) 


Then  go  to  (86) 

Kp  =  111.74  =  constant 


Then  go  to  (86) 


r* 

l 


If  0  <  T  t  <  210 
go  to  (83) 

* 


go  to  (84) 

* 


go  to  (85) 


* 


Kp 


(T'i  -  Tz) 


r  " 

1 

2.765 

2  2 
(T’l  )  -  (T'z) 

(T'i  )2,16  -  (T'a)2,16 

2 

. 

2. 16 

s 

(83) 


Then  go  to  (86) 


Kp  =  86.0  +  j  6.25 


(T’l  +  T‘z)  / 2]  -  210  I 
330  / 


(84) 


'hen  go  to  (86) 


Kp  =  92.25  =  constont 


(85) 


Then  go  to  (86) 


(86)* 
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If  0  <  T*2  <113  If  1 13  <  T'2  <  540  IF  540  <  T*2  <  oo 

go  to  (87)  go  to  (88)  go  to  (89) 


Then  go  to  (93) 


(90) 


tttW  fa  >2-M  - 


/TU  ,2.16  ,  2.16 

(^2  )  -  (Tcni) 


2.16 


Then  go  to  (93) 


Kp  =  86.0  * < 

/ 

6.25  ^  &T‘2 

\ 

+  Tcm)/^]  -210  \ 

Then  go  to  (93) 

330  j  j 

Kp  =  92.25  a  constant- 

% 

Then  go  to  (93) 

K  =  Kp  .93061  f  s  -  dh 
♦•5  -  (dh/z  )  " 


(91) 


(92) 


(93)* 


U  = 


1-2  2-3 


1 _ 

+  UTiJTX 

k  k  h 

3-4  4-3  5-6 


(94)* 
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If  Chi  <  Cc2 

go  to  (102) 


If  Chi  >  Cc 2 
go  to  (105) 


If  Chi  =  Cc2 
go  to  (108) 


*  3600  Chi 


(102)* 


£  - 


£ 


^V) 


(103)* 


TEST  7  -  A 


If  (103)  >  1.0,  stop  ond  readout  message  . "REDUCE  Ntui" 


II  II 


<  "  ,  continue. 


Ntu 


(104)* 


Then  go  to  LG  -  14 


If  Cht  >  Cc2 


*  _  Q? 

A  "  i460  Cc2 


(105)* 


£ 


(106)* 


TEST  7  -  B 

If  (106)  y  1.0,  stop  ond  reodout  message . 'DEDUCE  Ntui" 

"  "  <  "  ,  continue 


Ntu 


Then  go  lo  LG  -  14 


If  Chi  s=  Cc2  ±  10"6 


01 

n~  5600  Chi 


£  = 


TEST  7  -  C 

If  (109)  >1.0,  stop  and  readout  message  . "REDUCE  Ntui" 

11  "  ^  /  continue 


(107)1 


(103)* 


(109)* 


Ntu  = 


1 


-  £ 


a  io)* 


Then  go  to  LG  -  14 


177 


If  Chi  <  Cc2 
go  to  (111) 

I 


If  Chi  »  CC2 

goto  (114) 

* 

If  Chi  «  Cc2 


If  Chi  =  Cc2*  10 
go  to  (117) 

I 


Ax  tot.  hot  tide 


3600  Ntu.  Chi 


Axp  *  Ax  t  •  Af  i 


op  calculated 


144  Ax  tot.  hot  lid* 


Note  .....  If  a  fraction  reiults,  go  to  next 
higher  whole  number. 


TEST  8-A 


(1 13)  muit  «  (4) 


If  (113)  *  (4),  lncrea«e  (4)  and  iterate  from  (4). 


II  „  II  II  II  II  II 

f  rMUCt 


Then  go  to  (120) 
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If  Chi  >  Cc2 


Ax  tot.  hot  side  =  3600  Ntu  .  Cc2 _ 

U  (114)* 


Axp  =  Ax  I  •  Afl 


(115)* 


np  calculated  =  114  Ax  tot,  hot  side 

Axp 


(116) 


Note.  ....  If  a  f raction  results,  go  to  next 
higher  whole  number. 


Test  8  -  B 
(116)  must  =  (4) 

If  (116)  >  (4),  increase  (4)  and  inerate  from  (4). 
"  <M  ,  reduce  "  "  "  "  "  . 


Then  go  to  ( 1  20) 


If  Chi  =  Cc2  +  10  '6 


Ax  tot.  hot  side  =  3600  Ntu  ■  Chi 

U 


(117)* 


Axp  =  Ax  l  «  Al  I 


(118)* 
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l  .tlculaled  -  144  Ax  tot,  hot  side 

Axp 


(119) 


Note . If  a  fraction  results,  go  to  next 

higher  whole  number. 


Test  8  -  C 
(119)  must  =  (4) 

If  (119)  >  (4),  increases  (4)  and  iterate  from  (4). 
"  "  <  ”,  reduce  "  "  "  .  "  ", 


Then  go  to  (120) 


width  X  =  equation  (13)  =  inches 
after  closure 


(120)* 


heigth  Y 


Fs 


=  inches 


(121)* 


core  length  L  =  £(np  •  tp)  +  (ns  .  ts)J  =  inches 


(122)* 


core 


weight  -  .098  np  tp  |[(XY)  -  (X '  Y  ’ )]  +  £x  B '  (C  -  1)J  +  [a£ l  (Ra+1  )(1  -O')] 
.078  ns.  ts  [(XY)  -  (X*  Y')]  +[x*  B’  (C-l)]  =  lbs  (123)* 


header  weight  =  .196  £(XY)  -  Afi  (Ra  +  1)  +  -X-Y  _  J  = 


ibs 


024)- 


total  weight  -  (123)  +  (124)  Ibs 


(125) 
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(126) 


Test  9 

If  ( 1 26) < .  40,  readout  message .  "INCREASE  (IP  -  96). 

"  "  >.60,  "  "  "  .  "REDUCE  (IP  -  96). 


APPENDIX  VII 


HX  -  6  AND  TURBINE  3 PROGRAM . CALL  HX  (J) 

j 1 6,0 


Imputs: 


Cal!  out  numerical  values  from  APPENDIX  I,  section  VII. 
Also  call  last  result  for  T9,  P9,  T21  and  P21  from  output 
of  HX  -  5. 
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InH  i  a  hmn_H'  rn  a  1  assumptions: 


K.q'iu,  Initial 

Xu;  Value: 


(4) 

Kb  0 

(•’) 

T9  -  10 

((’) 

1 . 02  P9 

(H) 

2.  0 

(41  l 

T  2  1  ^  lo 

(42) 

.r»8  P21 

Notes : 

Readout  last  result  of  all  equations  marked  with  a  star  " 
Do  not  slop  machine  at  TEST  9. 


s  -"W.  906894  £|h^ 

v  (T 


(D‘ 


n 


4  cr 
rr  dh2 


(2)* 


Ax  I  -  (n  7Fdh  tp)  4  Z(  I  -O') 


(3>* 


np  =  assume 

(Will  be  altered  only  by  TEST  8)  (4)* 


To  -  assume 

(Will  be  altered  only  by  TEST  3)  (e»)« 


P8  =  assume 

fVVill  be  altered  only  by  TEST  2)  (6) 
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v  , 


v*  -  [{C-  1)  b'  ]  - 


Nh 


2  yC~  l) 


(16)> 


Afi  calculated  =  Nh(X'Y'i)  '  (17) 


TEST  1 

(17)  must  =  (12)  .  001 

If  (17)  >(12),  reduce  (13)  and  iterate  from  (13). 


=  Y'l  Ra 

2  (18)* 

Y^  =  Y*l  .  Ra  (19)* 

li-  -  yl' 

24  (20) 


(21) 


Y*3 

12  •  Ra 


(22) 


ip;  - .  ^r<tp/di,) 

'V*8  *  NRe  1  (23) 

A  Pi  =  np  •  A  Pi  (24)* 


P8  calculated  =  P9  +  APl  (25) 
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TEST  2 
(2  5)  must  = 
If  (2  5)  >  (6), 

II  H  <  IT 


(fa)  .  001 

increase  (fa)  and  iterate  from  (fa), 
reduce  "  " 


Call  cp  subroutine  and,  with 

T8  and  P8 . get 

cp  8  - -- 

Call  cp  subroutine  i 
T2 1  and  P21  .  .  .  ge 

ind,  with 

t 

— 

cp  21  - 


(26) 


(27) 


cp  =  1cE_8_  .+  cp  21) 

2  (28) 


(29) 


(30) 


T8  calculated  = 


_ T2  ]_ 

(T21/T8) 


(31) 
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TPST  3 

(31)  must  -  (5)  t  .00) 

If  (  31)  >  (5).  ru  r  tMsi'  (6)  and  iterate  (rum  (5). 
"  '■  <  "  .  ruduru  '*  "  '•  "  . 


m  mi 


S ,  5  ‘i-?  •)  7  x  10 


~i 


gs,1-643  - 

1.643 


-1 


(32) 


Km  i  -  _ ^T-v  ?q  x  ^ 


-3 


[ 


00355  ( T8  -  T9) 


■■’3  l 


1.642  1.642 

(.00355  T8)  -  ..00355  T9) 

1.642 


(33) 


Kml 

1  /.  642 

57.79  x  10-3  | 

.00355  ( 34) 


Phm  =  _iPL±P9j _ 

2  (35) 


j  Call  cp  subroutine  and,  with 
j  Thm  and  Phm . get 

cphm  =  ^ — 1 

(36) 

'.hi  =  cphm  • 

W3 

(37)* 

NPrl  =  cP^m  ‘ 

M  ml 

Kml 


( 38)* 
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NNu  I 


-  3.  66  +  ^ 

VLNKel  ■  N rr  1  J /  l. 

1  + 

.016  \ 

(ip/ dl.) 

LNRei  •  NPrij  / 

h  = 

1-2 


1Z  NNul  ■  Km  I 
dh 


COLD  SIDE: 


TZ2  -  assume 

(Will  be  altered  only  by  TEST  6; 


P22  -  assume 

(Will  bo  altered  only  by  TEST  4) 


If  6  <  T22  <40 


LG  -  2 


If  40  <  T22 


aht  •  Ra 


ah2 


(46);:- 


144  W4  V22 
ah2 


(47)* 


A  P2 


NRe2  = 


V22  '  dh 


121/22-  fJ-  22 


_  370  x  10‘6 


V  2  22 


V  22 


"V- 


(tp/dh) 

NRe2 


A  P2  =  np  •  A  P*2 


P22  calculated  =  P21  -  A  P2 


(48) * 

(49) 

(50) * 

(51) 


TEST  4 

(51)  must  =  (42)  +  .  001 

If  (51)  >(42),  increase  (42)  and  iterate  from  (42). 
"  ''  <  "  ,  reduce  "  11  "  "  , 


TEST  5 

If  (51)  <  10,  stop  k  readout  message . "INCREASE  Ra" 

’*  "  5  1 ,  continue. 


/imZ 


8. 55497  x  10 


- 

1.643  1.643 

\ 

J  (122)  -  (T2 1 ) 

1,643 

s. 

(T22  -  T21) 


(52) 


Km2 


57.79  x  IQ"3 
[.00355  (T22  -T21)] 


(.00355  T22)1,642  -  (.00355  T21)1'642 
1.642 


Tcm 


Km2  \  1/.  642 

57.79  x  10~3  / 

.00355 


Pcm 


(P21  +  P22) 
2 


Call  cp  subroutine,  and,  with 
Tcm  and  Pcm . get 


cpcm 


(53) 


(54) 

(55) 


(56) 


Cc2  =  cpcm  •  W4 


(57)* 


NPrz  . 

Km  2 


(58)* 


NNU2 


l 

.104  \ 

r  (tp/dh)  i  | 

J  ' 

NRe2  •  NPr2j/ 

's 

ToT5 

r 

‘  (tp/dh)  I  b 

Nre2  •  NPr2j 

) 

>> 

i 

y 

h  12  NNu2  *  Km2 
5*6  "  dR 


(59)* 


(60)* 
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f 

If  Chi  <  C>  Z 
t'1  (>■  1 } 


i 


<!  ■' .  -  $> 

V 

! !  Chi  >  ft  2 


n" »  ■  r 


1 

If  Chl  sCf  2  •*  10  "b 
go  to  (67) 


T22  calculated 


+ 


(T8  -  T9) 


(61) 


TEST  6-A 
(6 1)  mus  t  = 

If  (61)>  (41). 

ii  M  <  n 


(41)  +_  .001 

increase  (41)  and  iterate  from  ( 4 1 ) . 
reduce 


Ntut 


1  1 

11  1  * 

* 

log 

r (T8  -  T22)~ 

e 

L(T9  -  T2 1 )  _ 

(62)* 


e,  . 


i 


(63)* 


Then  go  to  (70) 


If  Chl  >  Cc Z 
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T22  calculated 


T2 !  + 


ILL  -  .151 


(64) 


TEST  6  -  B 

(64)  must  -  (41)  _+_  .00) 

If  (64)  >(41),  increase  (41)  and  iterate  from  (41) 
"  "  <  ",  reduce  "  "  "  "  "  . 


Ntui 


log 

(T8  -  T22) 

e 

L  (T9  -  T21)  J 

/Cc2/  \ 

1  -(  /Cb  ij 


(*>*>) 


(66) 


Then  go  to  (70) 


If  Ch  i  s  Cc  2  +  10-6 


T22  calculated  ^  T8  -  (1 9  -  T21) 


(6  7) 


TEST  6-C 

(67)  must  =  (41)  +  .  001 

If  (67)  >  (4 1 ) ,  increase  (41)  and  iterate  from  (41). 
"  "  <  "  ,  reduce  "  . . . 


191 


8  -  T9 


r  (T8  -  T9)  -I 
'  [_  (T8  -  T21 ) 


Ntui 
+  Ntui 


Then  go  to  (70) 


T ' 2  -  Thm 


■ 

Thm  -  (yv  *yV1  Y'i  )  (Tlw"  • Tcm) 

Thm  •  (v^Vf^w)  <Thm  ' Tcm) 


LG  -  4 


If  (IP-122)  CODE 
go  to  LG  -  5 


If  (1  P-1 22)  CODE 
go  to  LG  -  6 


If  0  <  Thm  <  113 
go  to  72 1 


If  !13<  Thm  <  540 
go  to  (73) 


If  540  <  Thm  <  OO 
go  to  (74) 


__  1 _ J  _49.  [/ 

|”  .  I  (Thm  -  T’i  )  J  St  2  (_'* 


1  Thm)2  -  (.  IT'  i)2 


.  i_  r  ( 

3.47  '* 


iTu  \3»47  , 

IThm)  -  (. 


Then  go  to  (78) 


(.  1  Thm)2  -  (.  1  T'l  )2 


Kp  = 


[' 


1 

/ 

r(w-70,-(.iT'„2’708  i 

+  9.551 

T7TU-TTV 

_ 

2.708- 

- 

Then  gc-  to  (78) 


(73) 


Kp  =  111.74  =  constant 


(74) 


Then  no  to  (78) 


Kp  =  86.0  +  { 


6.25 


Then  go  to  (78) 


(76) 


Kp  =  92.25  =  constant 


(77) 


Then  go  to  (78) 
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.  m  ,l.iih  n  ■! » .oi  4P, ^ 


K  =  Kp  .93061 

2-3 


i/8r 


1 


If  (IP  -  122)  CODE  -  1.0 
go  to  LG-8 


* 


f 


If  (IP  -  122)  COD':  =  2.0 
go  to  LG-9 


* 


i 


IF  0  <T’i  <  113  If  113<T'i  <  540  lf540<T'i<oo 

go  to  (79)  go  to  (80)  go  to  (81) 

I  r  J  I  . 

*  “  C.lft'/.T'a)]  {t-  [(■'  T''  >2‘  <•'  r’>2]-  T3T  t-'  r'  >3‘47-<-1  t'*>3‘47]  } 


(79) 


Then  g<*.  to  (85) 


Then  go  to  (85) 


Kp  =  111,74  =  constant 
Then  go  to  (85) 


(81) 
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(85)* 


If  (IP  -  122)  CODE  *=  2.0 
go  to  LG  -  12 

I 

95 


If  0  <  r*2  <  113  If  1 13  <  r'2  <  540  if540  <T'2<oo 

go  to  (86)  go  to  (87)  go  to  (88) 


(87) 


Then  go  to  (92) 


Kp.-  111.74  s  constant  (88) 

Then  go  to  (92) 


lf0<T2<'.10  If  210  <  T'2  <  540  If540<r2<oo 

go  to  (89)  go  to  (90)  go  to  (91) 


i 
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W9E 


Kp  = 


1 

(TV  -  Tern) 


i 


6.  I  b  3 

2 


-  (Tem) 


‘]- 


(T'z)2,16- (Tem) 

- m — 


(89) 


Then  go  to  (92) 


Kp  =  86.0  +/6.2s(H‘fe.-t.  If™?/2]  ~210  \  > 

1  \  330  / 

(90) 

Then  go  to  (92) 

Kp  =  92.25  =  constant 

(91) 

Then  go  to  (92) 

K  =  Kp  .93061 
4*5 


l  -dh 
t  -  (dh/2) 


(92)* 


*  l_  +  1 1  +  1 2  +  1 3  +  1  \ 

h  K  K  K  h  (93)* 

t  1-2  2-3  3-4  '  4-S  3*6/ 


ns  =  np  +  1 


(94)* 


It 


ns  •  ts 


12 
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(95) 


r 

If  Chl  <  Cct 
go  to  (101) 


TEST  7  - 
If  (102)  >: 

n  **  < 


Til  = 


(96) 


K  l 


7.27  x  10 


-3 


r 


(T12 


-  Til  ) 


m»  l1-585-  (Tig 

- T58T 


1.585 


(99) 


If  Chl 


Ql 

^  =  3600  cKT 


(101)* 


€=  e*« 


(102)* 


1 .0,  stop  and  reodout  message 


"REDUCE -BORDER  DIMS. 
OR  INCREASE  ts" 


",  continue. 


198 


-<“L) 


*‘L  1-6 

Ntu  =  _ _ 


■H„ ) 


( 103)* 


Then  go  to  LG  -  14 


I  If  Chi  >Cc2 


x  ^  Q1 

3600  Cc2 


(104)* 


€  - €i 


1  -  * 


(105)* 


TEST  7-B 

If  (105)>  1.0,  stop  &  readout  message. 
11  "  <  11  ,  continue. 


"REDUCE  BORDER  DIMS. 
OR  INCREASE  ts" 


log  1 


Ntu  = 


1  -£ 


m 


(106)* 


Then  go  to  LG  -  14 


If  Chi  E  Cc2  +  10  '6 


QI 

3600  Chi 


•199 


(107)* 


€  =  6  i 


(108)* 


TEST  7 -C  ''REDUCE  BORDER  DIMS, 

If  (108)  i  1.0,  stop  &  readout  message . OR  INCREASE  t.g" 

"  11  c  "  ,  continue. 


Ntu  =  _ ^ _ 

1-6  -  (109)* 


Then  go  to  LG  -  14 


Ax  tot.  hot  side  =  3600  Ntu  Chi _ 

U  (110)* 


Note . If  a  fraction  results,  go  to  next  higher 

whole  number. 
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TESI  h-A 

(1  12)  must  -  (4)  y 

If  (11L'.}>(4),  increase  (4)  and  iterate  from  (4). 
*'  "  <  reduce  "  Sl  " 


Note . If  a  fraction  results,  go  to  next 

higher  whole  number. 

TEST  8  -  B 

(115)  must  =  (4)  +q 

If  (M5)  >  (4),  increase  (4)  and  iterate  from  (4). 
1.  U  <  I>(  reduce  >'  11  "  "  " 


Then  go  to  (119) 
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Ax  tot.  hot  side 


3600  Ntu  Chi 

U 


016) 


Axp  Ax  i  •  Af  i 


(1 17)* 


np  calculated 


14 4  Ax  tat.  hot  side 
Axp 


018) 


Note . If  a  fraction  results,  go  to  next 

higher  wnole  number. 


TEST  8  -  C 
(118)  must  -  (4)  q 

If  (118)  >  (4),  increase  (4)  ond  iterate  from  (*>). 
"  •<  <  reduce  "  "  "  ”  " 


Then  go  to  (1 19) 


width  X 

-  equation  (13)  =  inches 

after  closure  . 

019)* 

height  Y 

X  •  u 

~  — -  =  inches 

Fs 

(120)* 
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3  np  K 
2 


<  <X  *  tp  .93061  — — —  } 

'31  L«-(dh/Z)|  I 


TEST  9 


If  (125)  <.40,  readout  message . "INCREASE  (IP  -  112) 

"  "  >.60  ,  "  "  . "REDUCE  (IP  -  112) 


Do  not  stop  machine  on  TEST  9. 


Av  = 


Ax  tot.  hot  side 


=  ft2 /ft3 


X-Y-L 


1728 


Cp  T8  rft 


/p2 1  /  \  (S*-  1)/?- 


P8 1 


BTU/lb 


Turbine  output  -  1054.54  Wt3  AHt3  r  Wat 


3 


FINAL  TEST 


if  (41)  ?.  540,  stop  machine 

If  (41)  <  540.  call  h/.  (J)  and  continue. 

j*7.0 


The  FLOW  DIAGRAM  for  HX  -  6  is  exactly  similar  to  that  of  HX-2. 


\ 
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APPENDIX  VIII 


HX  7  PROGRAM  ....  CALL  HX  (J)  j  =  7.c  : 


Inputs: 

Call  numerical  values  from  APPENDIX  I,  Section  VIII. 
Also  call  last  result  far  T8,  P3,  T22  and  P22  from 
output  of  HX-6. 


Initial  numerical  assumptions: 


Equa. 

Initial 

No. 

value: 

(4) 

300.0 

(5) 

T8  +  150 

(6) 

1.02  P8 

(13) 

2.0 

(35) 

T22  + 150 

(36) 

.98  P22 

Notes: 

Readout  last  result  af  all  equations  morked  with  a  stor 
Do  not  stop  machine  at  TEST  9. 
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HOT  SIDE: 


906894 


dh2 

S 


n  = 


4  cr 

7fdh  2 


Axl  =  (N7|"dh  fp)  +  2  (1  -O') 


np  -  ossume 

(Will  be  altered  only  by  TEST  8) 


T7  =  assume 

(Will  be  altered  only  by  TEST  6) 


?  7  =  ossume 

(Will  be  ottered  only  by  TEST  2) 


ju7=  2.37888  x  10’7  (T7)  ,643 


V7 


12  NRei  47-7  fi7 
dh 


(10)* 


ah) 


144  W4  4?  7 
V7 


(11)* 


Afi  = 


oh  i 

<r 


(12)* 


X  =  assume 

(Will  be  altered  only  by  TEST  1) 


(13)* 


X'  =  X  -  (2  Bx) 


(H)* 


(15)* 


Y'l 


-  L(C-1) 
2  (C-l) 


'L. 


Nh 


(16)* 
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Afi  calculated  “  Nh  (X1  Y|l) 


(17) 


TEST  1 

(17)  must  =  (12)  *.001 

If  (17)  >  (12),  reduce  (13)  and  iterate  fram  (13). 

ii  ii  S  li  11  11  11  11 

"  "  <  1  ,  increase 


Y‘2  -  Y'l  Ra„ 

2 

(18)* 

Y'3  =  YV  *  Ra 

(19)* 

>- 

n 

(20) 

lt  =  “TT 

(21) 

»*■  A 

(22) 

lo”6  V72  /  (tp/ dll) 

#7  V  NRel 

(23) 

API  ■  np  •  AP'l 


(24)* 


P7  calculated  *»  P8  +  APi 


(25) 


TEST  2 

(25)  must  =  (6)  ± .  001 

If  (25)  >  (6),  increase  (6)  and  iterate  from  (6). 
"  "  <  ",  reduce  "  "  "  "  "  . 


Chi  =  cphm  *  W4  (31)* 
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NPr  I 


NNu!  3 


COLD  SIDE: 


cphm  ■  ^mi 
Km  I 


'32)* 


(33)* 


h  = 
1-2 


12  NNu  I  ‘  Km  i 
dFt 


(34)* 


T23  =  assume 

(Will  be  altered  only  by  TEST  5) 


(35)* 


P23  =  assume 

(Will  be  altered  only  by  TEST  3) 


(36)* 
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P23  calculated  -  P22  -  AP2 


(45) 


TEST  3 

(45)  must  =  (36)  *.001 

It  (45)  >  (36),  ir.creose  (36)  and  iterote  from  (36). 
tt  "  <  '  ,  reduce  . , 


TEST  4  "DECREASE  NRei  or 

If  (45}  <  10,  stop  &  readout  messoge  .....  INCREASE  Ra  " 

"  "  >  10,  continue. 


_  8.55497  x  10“^ 

*mZ  "  (123 -W" 


(T23)1*643-  (T22) 
1.643 


1.643 


(46) 


Km2 


57.79  x  10"3 

(.00355  T23)1  642  -(.00355  T22)1'642 

f  .00355  (T23  -  T22) 

L 

1.642 

(47) 


* 


Tcm 


t  1/.642 

Km2 

57.79  x  10~3) 

.00355 


(48) 


Pcm  =  JWJ-m 


m 


Cali  cp  subroutine  and,  with 
Tcm  and  Pcm . get 


epem  - 


(50) 
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TEST  5-A 
(58)  must 
If  (58)  > 

II  11  < 


TEST  6-A 
(59)  must 
If  (59)  > 

n.  H  < 


AX 


=  V]  ’  (78  1T22) 

Ccz)\] 

Chi  1 


'4*.  wm  saw!  uin  i*  nup. 


(V 


■nvminifH 


AY 


(Ccz  \ 

Chi  I 


•  AX 


(57) 


T23  calculated  =  T22  +  AX  (58) 


:  (35)  ±.001 

(35),  increase  (35)  and  iterate  from  (35). 
"  ,  reduce  "  . . . 


T7  calculated  =  T8  +  AY  (59) 


=  (5)  ±.001 

(5),  increase  (5)  and  iterate  from  (5). 
"  ,  reduce  .  "  " . 


-  Ntui  fl~(Chi/Cc2)l 

1  -e _ 1 _ J  (60)* 


{ 

1  [  Chi  I 

-  Ntui  p  -  (Chi/Cc2)j 

p 

Cc2 

V' 

Then  90  to  (71) 
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I - 1 

If  Chi  >  Cc2 


VI 


Ntui 

(78  -  T22)  e 


(60 


AX 


VI  -  (T8  -  T22) 


Chj_  _  i 
CTc2 


(62) 


AY 


Chi 


Cc2 


AX 


(63) 


T23  calculated  =  T22  +  AX  (64) 

TEST  5-B 

(64)  must  =  (35)  ±.001 

If  (64)  >  (35),  increase  (35)  and  iterate  from  (35). 

"  11  <  ",  reduce  . . 


T7  calculated  -  T8  +  AY  (65) 


TEST  6-B 

(65)  must  “  (5)  ±.001 

If  (65)  >  (5),  Increase  (5)  and  iterate  from  (5). 
"  "  <  ",  reduce  "  "  . . . 


•  _ 
a  - 


1“ 


-  Ntui 

1  -  e 

1-  (Cc2/Chl)] 

Cc2 

Chi 

*  Ntui 

I  6 

r1~(Cc2/Chi)j 

(66)* 


Then  go  to  (71) 


215 


|  If  Chi  s  Cc2  ±  10 


(T8-T22) 


1  +  Ntui 


T23  calculated  =  (T22  +  2)  -  (T8  -  T22) 


TEST  5  -  C 

(68)  must  =  (35)  ±  . 00 1 

If  (68)  >  (35),  increase  (35)  and  iterate  from  (35). 


■I  n  <  » 


,  reduce  "  " 


tt  (I  II 


T7  calculated  =  T22  +  Z 


TEST  6  -  C 


(69)  must  =  (5)  ±  .001 


It’  (69)  >  (5)  ,  increase  (5)  and  iterate  from  (5). 


"  "  <  ",  reduce 


ii  ir  ii  II  n 


1  +  Ntui 


Then  ga  ta  (71) 


T'i  =  Thm  -  |  y*'|  +  Y'~a)  (Thm  "  Tcm) 


T‘2  --  Thm  - 


Y’i  +  Y'2 


If  (IP  -  139)  CODE  *1.0 
go  to  LG  -  6 


(Thm  -  Tcm) 


If  (IP  -  139)  CODE  *2.0 
go  to  LG  -  6 


? 
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Kp 


223 


Ntu  - 


(1 


i 


log 

e 

1-* 

(Cc2/Chi) 

1  -« 

'-I 

f  Cc  2  / 

i  /ChiJ 

1 

Then  go  to  LG-14 


If  Chi  =  Cc2  ±10 


A  = 


Ql 

3600  Ch  i 


€  -  4 


TEST  7-C 

If  (109)  >1.0,  stop  &  reodout  messoge 
"  "  <  "  ,  continue. 


....  “REDUCE  Ntui 

l 


Ntu 


€ 

1  -  £ 


Then  go  to  LG-14 


(1 


(1 


( 


”1 
If  ChlS  Cc2  ±  10“ 
go  to  (1 17) 


F 


If  Chi  <  Cc2 


Ax  tot.  hot  side  =  ,.?6.PP.  Ntu ,  Ch_l (111)* 


Axp  =  Ax!  '  Af  I 


np  colculoted  =  Wft  t*-..h°t  »ld* 
r  Axp 


012)* 


013) 


Note. ...  If  a  fraction  results,  go  to  next 
higher  whole  number. 


+1 


TEST  8-A 

013)  mutt  =  (4)  ^ 

If  (113)  >  (4),  increase  (4)  ond  iterote  from  (4). 
"  "  <  "  ,  reduce  "*  "  "  “  "  . 

Then  go  to  (120) 


If  Chi  >  Cc2 


Ax  tot.  hot  side 


3600  Ntu  Cc2 
U 


(114V 


Axp  ~  Ax  I  •  Afi 


(1  IS)' 
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* 


np  calculated 


144  Ax  tot,  hot  side 
Axp 


(116) 


Note .  If  a  froction  results,  go  to  next 

higher  whole  number. 


TEST  8-B 
(116)  must  =  (4) 

If  (116)  >  (4),  increose  (4)  and  iterate  from  (4). 
"  "  <  "  ,  reduce  "  *'  "  "  " 


Then  go  to  (120) 


!fCh|  =  Cc2  ±10"6 


Ax  tot.  hot  side  = 


(1 17)* 


Axp  =  Ax  I  *  Aft 


(118)* 


np  calculated  =  144  hot  side  (119) 

Axp 

Note . If  o  fraction  results,  go  to  next 

higher  whole  number. 


TEST  8-C 
(1 19)  must 


(4) 


+1 

-0 


If  (119)  >  (4),  increase  (4)  ond  iterote  from  (4). 
"  "  <  "  ,  reduce  11  "  "  "  " 


Then  go  to  (120) 
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width  X  =  equation  (13)  = inches 
after  closure 


(120)* 


V 

height  Y  =  — —  “  inches 

(121)* 

care  length  L  =  £(np  *  tp)  +  (ns  *  ts)j  =  inches 

(122)* 

core 

weight  = 

.098  np  tp  j  (XY)  -  (X  'Y')]  +  [x'B1  (C-l)]  +  [Afi  (Ra+1)  (l-o')]j  + 

.  .078  ns  ts  j[(XY)  -  (X'Y’)]  +  [x'B1  (C-l)]j  =  lbs 

(123)* 

header  weight  =  .  196^(XY)  -  Afl  (Ra+1)  +  J~C  j  =  lbs 

(124)* 

/ 

total  weight  **  (123)  +  (124)  =  lbs 

(125)* 

TEST  9 

If  (126)  <  .40,  readout  message . "INCREASE  (IP  -  132)" 

"  "  >  .60,  "  "  . "REDUCE  "  " 


Do  not  stop  machine  on  TEST  9 


Ax  tot,  hot  side _ _ »  f*2  /  ft3 

7~x  •  y  •  i\ 

\  1728  I 


FINAL  TEST 

Final  test  does  not  opply 


STOP  ofter  (127). 


The  FLOW  DIAGRAM  for  HX-7  is  exoctly  similor  to  thot  of  HX-1 . 


(127)* 
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APPENDIX  IX 


BASIC  FORTRAN  MACHINE  LISTING: 


The  following  comprises  the  Fortran  machine  listing  for  the 
complete  program. 
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G 

C 

C 

c 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 
1  6 

17 

18 

19 

20 


i  iM  IN 


fom  t if; c  vnu  a::.',s_y  :  i cai_ 


FC1I  CRYOGENIC  TURBOMACMIi 


SYSTEMS 

VE>.'£  i  1  7JL_':  19  70 

DIMENSION  I  <t»)/AHl  (8  ) >  A F  1  <83 

DiM k: :  s i ok  xpr  cg  > /  ypr  < > /  anp  cs  )  *  dp  i  t s.)  /  ampri  c r  > 

DIMENSION  <:-,•>>  VMUUl  C  S.'3  /  H 1  2  <  3. 3  /  AMR  E2  <  8  ) 

Di  MEN  SI  OK  LP  2  < »  3  /-  AftP  r:  2  <  3  •>  /  ,*• ;  «•:  1 72  (C  3  /  H  5  6t  6  > 

DIKEN5I3H  A K 2 3  <  3  3 / AH  3 0 ( b  3 , AU  4  5  <  3 i  / U <  8  3 
DIMENSION  AXTKS(S)/ AXP<3)/ LHB)/AN$<8) 

DIMENSION  01,  CS)/AE<3)/  VX  C 3 3 / HY < 8 3 > ALL  Co  3 
Diii ENSIGN  C ! .’<b  3  /  K '/( R  3  /  TV.'( 8  3  /  AN F( 8  3  z  A V<  8  3 
DIMENSION  CHI  (8  )  /  CC2<8  )>  S2<  8  3/  AN2<8 3  /  AH2<8  ) 

Dll! FUSION  K  I  C  8  3  >  K2I  <8)/K3(8)zK4(3) 

DIMENSION  AN  TUI  3 3  /DHT<3)j AKLT(8)/ANTUI  <8 3 /  AKP <  5) /  TLX <  5) 
DIMENSION  XPK1 (33/ YPRl <S 3/ YPR2<3 3/ YPK3<8 ) 

DIM FUSION  ALANC»)/TX<3)  .  . 

COMMON  St  AN/AX  1  /  AH  I  /  AFI  /  XPR/  YP  Hz  ANP  t  DP  I  /  ANP  R I  t  ANPR 

COMMON  ANNU1 .Ml 2/ ANRE2/ DP2/ANPR2/ANNU2/H56/  AK23/AK24 

COMMON  AK45/ Uz AXTHSz  AXPz  El /  ANS/QL/AE/  LX/HY/ALL/ CW/HL’z  T  V 

COMMON  AHF/AV/ CHI t  CC2/52zAN2*AH2/Ki/K21/KCl/K4 

COMMON  T I  4/ Pi  4/  VBI  4/H!  4/ SI  At  fc.'l 

COMMON  Tl  5/PI  5/  VBl 5/ HI  Si  SI  5 

'COMMON  Tl  5/Pl  6/  VBI  6/Hl  6/SI  6  ■  . 

COMMON  DHz  TP»  51 G/ ANREl / C/  FS/RF/  ANTU/T5 
COMMON  ANTI  *  DH2z  S I  G1  /  SI 02*  V, 2/  WT1 


COMMON  DHT  t AK  WT *  ANT  U It M A  Tz  TLX 
COMMON  BX/BYzANK/AMC 
COMMON  XPRI t YPR1/YPK2/YPRS 
COMMON  ALAM/TX/BPR/RA 
FORMAT < lXz I 8HPR0GRAM  INPUT  DATA 3 
FORMATt IX# 1 1  HI  LOAD  DATA) 

FORMAT  Cl  X#  22H  AT  JT  VALVE  INLET*  3 
FORMAT < lXz  5H  Tl  4= FI 5.6/2X/ 5HDEG«R> 

FORM A T  < I X  z  SH  F I  4=  FI S . 8  z  2X  z  4HP S I A  3 
FORMAT < lXz  6H  VBI 4= Ft  5.8/ IX/9HCU.  FT* /LB 3 
FQRMATUX/5H  HI  4*F!  5.B/2X/ 6H81U/LB) 

F0RMAT<  I  Xz  5H  Sl4*Fl  5.8/ 2Xz  1  2HBTU/LB-  DF.G»R> 
FORMAT < 1 X#  4H  Vl«F15.8> 

FORMATC/lXz  70< 1H$)  3 
FORMAT< 1X/2IH  AT  JT  VALVE  EXIT! 3 
FORMA T< tXz  5H  Tl 5= FI  5 .8/ 2Xz 5HDEG. R 3 
FORMAT ( IX# 5H  Pi  S*F1  5.B/2X/4HPS1A3 
F0RMAT< 1 X* 6H  VB1S*F1 5.8/ 1X/9HCU. FT-/LB) 
FORMATt IX/ 5H  HI 5=F! 5-8/2X/ 6HBTU/LB) 

FORMAT < I Xz  5H  SI  5*Fl  5.8/2X/  I  2HBTU/LB-DEG'R) 
F0RMAT< IX/ 1  TH  AT  LOAD  EXIT*) 

FORMAT < 1 X/ 5H  Tl 6=F1 5.8/2X/ 5HDEG.R) 
F0RMAT(1X/5H  PI  e*FI  5.8/2X/  4HPSIA) 

FORMAT < IX/ 6H  VBI 6=  FI  5.8/ IX/9HCU. FT./LB3 
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21  FORMAT? IX,  5K  HI  6®  FI  5 .8,  2X,  6H3TU/LB) 

22  F0 RMAT <  1 X  >  5H  SI  6=  FI  5.  3,  2X,  1  2HBTU/LB-DEG.R) 

23  FORMAT? IX#  9H  HX1  DATA) 

2 4  FORMAT? 1 X, 5H  T1  7=F1  5.8, 2X,  5HDEG.  R) 

25  •  FORMAT? IX,  5H  DH=  FI  5. 3, 2X,  6HINCHES) 

26  FORMAT ? IX, 4H  TP  =  F15.8) 

27  FORMAT ?1  X,  4H  TS-F15.8) 

28  FORMAT? IX,  7H  SIGMA=F15.B) 

29  FORHATt IX, 6H  NRE1=F15.8) 

30  FORMAT <  1 X,  3H  C=F15.8) 

31  FORMAT? IX, 4H  FS*FlS.8) 

32  FORMAT (IX, 4H  RF-F15.8) 

34  FORM AT (IX, 1 6H  0UTPUT  FROM  HX1) 

55  F0RMAT(1X, 25H0UTPUT  FROM  HX24TURB2NE  I) 

56  F0RMAT< IX,  5H  P12=F15.8,5H  V12=F15.S,5H  T12=F15.8) 

57  FORMAT? IX,  5H  PI  3=  FI  5.  8,  5H  VI  3s FI  5«8,  SH  T13=F15.8) 

58  FORMAT? IX,  5H  P17=F1S.S,5W  "1  7=F1  5.8,  5H  T1  7®  FI  5.8) 

59  F0RMAT(1X,5H  P18=F15.8,5H  V18=Fl5.8,5H  T18-F15.8) 

60  FORMAT? IX,  5H  V14-F15.8) 

61  FORMAT? IX,  5H  NTt«ri5.8) 

62  FBRMATOX,  5H  DH2-F1  5.8,  2X,  6H INCHES)  • 

63  F0RMAT(tX,8H  SlGrtAl -Fl 5.8) 

64  F0RMAT(tX,8H  SIGMA2-F1 5.8) 

65  FORM  T(1X, 3H  WT, II , 1H= FI  5.8) 

66  FT  *T?1X,  2H  W,  1 1 ,  1H-F1  5.8) 

67  AT(1X,9H  HX2  DATA) 

68  loRMAT?lX,5H  VI 3*  FI  5. 8) 

69  F0RMATC12I2) 

70  FORMAT?  1X,9H  HX3  DATA) 

71  FORMAT?  IX,  16H  OUTPUT  FROM  HX3) 

72  FORMAT? IX,  5H  NTU-F15.8) 

C  LOAD  DATA* AT  JT  VALVE  INLET 
OPEN?  3,  INPUT,/AP DAT/3 

100  FORMAT?  5? IX, FIS. 8)) 

101  F0RMAT?2I 5) 

102  FORMAT?  IX,  1  5’riTl  4,P1  4,  Tt  6,P1  6) 

103  F0RMAT?1X,  9HT1  7,ANRE1  ) 

104  FORMAT? IX,  12HG0M  ,N0G0«-1  ) 

105  FORMAT? IX, 10HANT1 ,ANREl ) 

106  FORMAT?  IX,  1 1HANRE1  ,ANTUX> 

11-0 

12-1 

13-3 

TLX?1)«113» 

TLX?2)-2tO. 

TLX?3)»0.  ** 

TLX?4)-0. 

TLX?5)*0V 

READ?  13, TOO) Tt  4, PI  4, HI  4, SI  4 
READ? It, 101)  KR,MAT 


C  AT  JT  VALVE  EXIT 

KFAPCI3.  tOOTl  5*  Pi  5*  Ml  5#  SI  5 

C  AT  L0AD  EXIT 

REAi  13.  IG0JT1  6#P1  6#Hl  6#  SI  6 
Re  "DC  l  3.  100JDUL/ALW 
WRITEC 12#  1) 

WRiIE«I8#  2) 

WRI TEC  12  >  3) 

WRITEC 12#  4) T1 4 
WRITEC 12#  S)P1 4 
WRITEC 12#  COVBI4 
WRITECI2#  7)H14 
WR I TEC  12#  fci ) SI  4 
WRITEC 12#  9  )  m 
WRITECI2# 10) 

WRITECI2# 1 1 ) 

WRITEU2#  12)T1  5 
WRITEC 12# I 3)P1 5 
WRITEC 12# 1 4) VBI 5 
WRITEU2#  1  5)H1  5 
WRITECI2# 1 6) SI  5 
WRITr  ’■  12#  1  0 ) 

WRITECI2# 1 7) 

WRITEC 12# 1 3) T1  6 
WRITECI2#  1  9  '>  P 1  6 
WRIT£CI2#W)V816 
WRITEC 12#  21  )H1 6 
WRITECI2#22)S1 6 
WRITEC 12# 10) 

C  HX1  DATA 

READ! 13# 69)  J 

READCI3#  1 00 )  DH#  TP#  TS#  SI  G#  ANREI 
READC13# I00)ANTUICJ)#C#ANH#ANC# FS 
READCI3# i 00 ) RA#  RF#  BX#  BY#  BPR 
READCI3#  i 00) Vil 
WRI TEC  12*23) 

WRITEC 12#  PS) DH 
WRIT  El  iii#  26#  TP 
WRIT EC  12# 27) TS 
WRITEC i 2# 28) SI G 
WRITEC IS# 39) ANREI 
wKITEC 12# 30 )C 
WRITEC 12# 31  )FS 
WRITE* I2#32)RF 
W5i  1  TEC  1 2#  24)  T1  7 
WRITEC IS# 10> 

READC  U#  »00)XI#ANPI 
SEl  4=1. 6«'0 
SFI  6=2-925>S 

1 000  7  CONTINUE 


* 


- 


r 
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C*********+**INITIAL  CONDITIONS*********** 

T13I=T14+10. 

PI 31=1 »02*P1 4 
T1 7I*T1 6+10 • 

°1 71*0.  98*P1 6 

CALL  HX1  *  1  #  P]  3#  VI  3#  T1  3#P  1  4#  VI  4# T1  4#  PI  6#  VI  6#  T1  6# PI  7/  VI  7*  T1  7# 
1  XI#ANPI#P1  3I#T13I#P17I#T1  7I#CPCM#CPHM) 

SE1  3*1  •  66+<  <CPHM*AL0G*T1  3/Tl  4>)-«l  .98  718*AL0G<P1  3/PI  4))) 
SE1  7=2.9225+* <CPCM*AL0G<T1  7/T1  6> )- <  1 .98  71 8*AL0G*P1  6/P1  7))) 
WRITE*I2#10> 

WRITE*I2#34) 

IF*KR>900# 1900*901 
901  .  G0  T0  903 

900  CALL  OUTPUT*  1) 

90  3  CONTINUE 

WRITE*  12#  57>P  1  3#  VI  3#  T1  3 
WRITE*  IS#  58)P1  7#  VI  7#T1  7 
WRITE*  12#  60>  VI  4 
WRITE* 12# 10  > 

1900  CONTINUE 

WRITE* 12# 104) 

READ*  11 # 101 )  NGO 
IF*NGO)1000  7#1 000  7#  10008 
10008  CONTINUE 
C  HX2  *TURB1NE  DATA 
READ*  13#  69)  J 

READ*  13#  100)DH#  TP#TS#SI  G#  ANRE1  . 

READ* 13# 100 )C#ANH# ANC# FS 
READ*  13#  1 00  )RA#RF#BX#BY#  BPR 
READ* 13# 1 00 ) WT1 # W2# ANTI 
WRITE*  12#  73)  J 
WRITE* 12#  61 )ANT1 
WRITE*  I2#25)DH 
WRITE* 12# 62) DH2 
WRITE*  I2#26)TP 
WRITE*  I2#27)TS 
WRITE* 12#  63) SI G1 
WRITE* 12# 64)SIG2 
WRITE*  12#  29)ANRE1 
WRITE* 12# 30 )C 
WRITE* 12# 31  )FS 
WRITE* 12# 32)RF 
WRITE* IL# 66) J# K2 
UXA*J-1 

WRITE* 12#  65) JXA#  WT1 
WRITE* 12# 10 ) 

READ* 13# 100 )X1# ANP1 
10005  CONTINUE 

C*********INITIAL  conditions******* 

Ti2I=T13+5. 
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P 1  21=  1  «02:l-P  1  3 
Tl  8i=Tl  7+5. 

PI G 1=0 • 98  '1-P  1  7 

CALL  HX2C2,P12,V12>T12,P13,  V13,T1  3,  PI  7,  VI  7,T1  7, PI  8,  VI  8,  Tl  8, 
1  X I , ANP I , P 1 2 1 , T 1 2 1 , P  1  G I ,  T18I>CPCw>CPK:-  ) 

SE12=SE1  3  +( (CPHii :  ALO^Cl  12/T1  35 5-  (1 .58  718*ALG G(Pl  2/P1  3>  )  5 
SE18=SE1  7  ♦  <  C  CP  CM*ALQGC  Tl  8/T1  7>)-*l .  98  71 8*AL0  G<  P 1  7/P18))) 

WRITE*  12,  1 0  ) 

WRITEU2,  74)  J 
IFCKR)904,1  904,905 

905  GO  TO  90  6 

904  CALL  OUTPUT < 2) 

906  CONTINUE 

1  WRITE*I2,56)P12,  VI2, T12 
WRITE(I2,59)P1S, Vl8,TlB 
WRITE*!2, 68 ) VI 3 
WRI TE< 12, 1 0 ) 

1904  CONTINUE 

WRITE*  12#  104) 

READ*  11,101)  NGO 
I  F(NGO)  1000  5,  1000  5,  1000  6 
10006  CONTINUE 

C  HX3  OATA 

REAO* 1 3, 69) J 

READ*  13,  100  )  OH,TP,  TS,  St  G,  ANRE1 

REAO* 1 3, 1 00  )ANTUI * J), C,ANH,ANC, FS 

REA0CI3,  100)  RA,  RF,  BX,  BY,  BPR 

REAOC 1 3,1 00 )W2 

WRITE*I2, 73) J 

WRITE* 1 2, 25) DH 

WRITE*I2,26)TP 

WRITE*I2,27)TS 

WRITE* 12, 28 )  SI  G 

WRI TE* 1 2, 29 )ANRE1 

WRI TE(  1 2, 30 ) C  -  'N 

WRITE*I2,31 )rS 
"  WftnEfrg,32)RF 

WKITECI2, 72)ANTUI*J) 

WRITE*  1 2,  10) 

READ*  13,  l 00) XI, ANP I 

C*>Mi******INITIAL  CONDITIONS************** 

TllI«T12+7S.  - . . 

Pll I«l .02*P12 
T19I*T11*75. 

P191»0.98*P18 
10003  CONTINUE 

CALL  HX 1*3, Pll,  VI 1,  Til, PIC,  V12,  T12,P1 8,  V18,  Tl  8,Pl  9,  V19,T19, 
1  XI, ANPI, Pill, Till, P19I,T19I,CPCH,CPHM) 

SE1 1  **SF1  2  +* *CPHM*AL0G*T1 1/T12) )-C 1 .98 718*AL0GCP1 1 /P12) ) ) 
SE1  9°SE18  +C  *  CPCM*AL0G<  Tl  9/T1 8)  )•<!  .98  71 8*AL0GCP18/P1 9)  ) ) 


234 


KRITECI2# 1 0  > 

73  FORMAT?  IX*  3HMX-#  12*  5H  DATA) 

74  FORMAT ( 1 X#  1 5H3UTPUT  FROM  HX-# 12) 

WRITE? 1 2# 74) J 

IF?KR)907#  1907#903 
90S  GS  TO  90  9 

90  7  CALL  OUTPUT? 3) 

909  CONTINUE 

WRITF.?I2#  75)P1  1#  V1UT1 1 
KRITECI2# 76)P19* V19#T19 
WRITECI2#  77) VI 2 

75  FORMAT? 1 X#  5H  P11  =  F15.8#5H  V11  =  F15.8*5H  T11  =  F15.8) 

76  FORMAT? IX#  5H  P19=F15.8#5H  V19=F15.8#5H  T19=F15.8) 

77  FORMAT? 1 X#  5H  VI 2= FI  5.8) 

1907  CONTINUE 

WRI TE? I 2#  10  4) 

READ?  1 1  #  1  0  1  )  NG0 
IF?NG0)10003#  10003#  10004 
10004  CONTINUE 

WRITE? 12# 10) 

IF?T19-  540*)50  4#  999#  999 
504  CONTINUE 

C  MX 4  DATA 

READ? 1 3#  69)  J 

READ?  1 3#  1 00')  DH#  TP#  TS#  SI  6#  ANRE1 
READ?  13#  100)C#ANH#ANC#  FS 
READ?  1 3# 100 )RA#RF#  BX#  BY#  BPR 
READ?  13#  100)  W2#  WT2#  W3#ANT1 
WRITE?  12#  73)  J 
WRITE? 12#  61 >ANT1 
WRITE? 12# 25) DH 
WRITE? 12# 62) DH2 
WRITE? 12#  2 6) TP 
WRITE? I 2# 2 7) TS 
WRITE? 12# 63) SI G1 
— WRITE? 12# 64>SIG2 
WRITE? 12# 29>ANRE1 
WRITE?  12#  30  >C 
WRITE? 12# 31 )FS 
WRITE? 12#  32)RF 
WRITE? 12#  66) J#  W2 
JXA* J- 1 

WRITE? 1 2#  65)JXA#WT1 
WRITE? 12# 10) 

READCJ.3#  100)XI#ANPI 
10504  CONTINUE 

C****s***#lNITIAL  CONDITIONS************ 

T10I*T1 1  +10  • 

P10IM  .02*PV1 
T20I*T1'9  +  10  • 
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r1  PO I  “  .  ■'■jvP:  9 

CALL  A;  -.'(4,F10j  VIC-  TIC.. FI  1  ,  VI 1  *  T 11  , PI  9,  VI 9* T1 9>P20,  V20,T20> 
1  XI,  ANPI, P 10  I  >  T10 I, P20 1 , T20 1  ,  CP CM, CPHM ) 

SE1 0=  SiT  1  1  +{ CCPHM*AL0GCT1O/T11  )  ) - C 1 . 93 71 8*AL0 GCP 10/P l 1))) 
SE20*SE19  +C  CCPCH*AL0GCT2O/T19) )-Cl . 98 71 S*AL0 GC?1 9/P20 ) ) ) 
WRT  TECT  2,  10) 

WRI TEC  12, 74) J 
IFCKR)910, 1910,91 1 

911  GO  T0  912 

910  CALL  OUTPUT (4) 

912  CONTINUE 

456  FORMAT C! X,  5H  P10=F15.8,5H  V10*F15.8,5H  T10-F15-8) 

459  FG RMA T  C i X ,  5H  P20=F15.8,5H  V20*F15.8,5H  T20*F15.8) 

468  .  F0RMATC 1 X,  5H  VU«F15.8> 

WRITEU2,4561P10,V10,T10 
WRITEC 12,  459  )¥20>  V20,  T20 
WRITECI2, 468) VI 1  . 

1910  CONTINUE 

WRITEC 12,104) 

READC 11,101)  N6C  .  . 

IFCNG0MO5U4,  10504,  10002 
10002  CONTINUE 

WRITECI2, 10) 

I  F(  T20  -  540*  )50  5,  9  9  9,9  9  9 
505  CONTINUE 

C  HX5  DATA 

REA0CI3,  69)  J 

READC  1 3,  100 )  DH,  TP,  TS,  SI  G,ANRE1  . 

REAOCl 3,  100 )ANTUI ( J), C, ANH,ANC, FS 

REA0CI3,  100)RA,RF,BX,BY,  BPR 

READC  13,  100  )W3 

WRITEC 12, 73) J 

WRITEC 12, 25) DH 

WRITEC 1 2, 26) TP 

WRITEC 12, 27)TS 

WRITEC 12, 28) SI G 

WRI TE  C  1 2, 29) ANRE1 

WRI TEC  1 2, 30  >  C 

WRITF.Ci£,31  )FS 

t.:n  fTr  #  t  rt  r»p 

I  t*\  k  G.*  UL  /  im 

WRITECI2, 72>ANTUICJ) 

WRITECI2>10) 

READC  13,  100)XI, ANPI 
T9I  -T10+100. 

P9IM .02*P10 
T21I=T20+100. 

P21 I=0.98*P20 
10  50  5  CONTINUE 

CALL  HX 1 C  5, fj 9  ,  V9  ,T9  ,P10,  VI 0,  T10#P20,  V20#  T20,P21>  V2l  ,  Tei , 
1  XI  ,ANPI,P9I,T9I  , P21 I,T2l I,CPCH,CPHH> 
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5E9=SEiO  +  !  !CPHH*AL0G!T9/TIO))-!1.98  718*AL0G!P9/PIO>)> 
SE21=SE20  ♦  !  !  CPCM*AL0G!  T21 /T20)  )  -  (1  •  98  718*AL0  G<  P20/P2 1  )  >  ) 

WRITE! 12# 10  ) 

WRITEU2#  74)  J 
IF  KR)913#  191  3# 9 1  4 

914  G0  T0  915 

913  CALL  0UTPUT <  5) 

915  C0NTINUE 
WRH‘E<I2*575>P9*  V9*T9 
WRITE!I2#576)P21#V2l#T2t 
WRITE! 12#  577> V10 

575  FORMAT! IX#  5H  P9  =  F15.8#SH  V9  =  F1S.8#5H  T9  =  F15.8) 

576  FORMAT!  1 X*  5H  P21  =  F15.8#$H  V2l»F15.8#SH  T21  =  F1S.8> 

577  FORMAT! IX#  5H  V10*F15.8) 

1913  C0NTINUE 

WRITE!  12#  10  4) 

READC11# 101 )  NG0 
I  F!NGG>  1 0  50  5#  10  50  5#  1 000 1 

10001  CONTINUE 

WRITE!  12#  10) 

1FCT21-540.  )50  6#  999#  999 

506  CONTINUE  " 

C  HX6  DATA 

READ!  I  3#  69)  J 

READ!  1 3#  100  )DH#  TP#  TS#  SI  G#  ANRE1 
READ!  I  3#  100)C#ANH#  ANC#  FS 
READ!  I  3#  100  )RA#  RF#  BX#  BY#  DPR 
READ!I3#  100)  W3#  WT3#  W4#ANT1 
WRITE! 1 2#  73) J 
WRITE!  12#  6DANT1 
WRITE!I2#25)DH 
WRITE! 1 2#  62 ) DH2 
WRITE!  1 2#  26) TP 
WRITE!  12#  27)  TS 
WRITE! 12#  63) SI G1 
WRITE! 1 2# 64) SI G2 
WR1TE!I2#29)ANRE1 
WRITE! 12#  30 )C 
WRITE!I2#  31  )FS 
WRITE! 12#  32)RF 
WRITE! 12# 66) J# W2 
JXA* J- 1 

WRITE! 12#  65)JXA#  WTl 
WRITE! 12# 10) 

READ!  1 3#  100)Xi#ANPI 

C***+****-1'*9iNi  T1AL  CONDI TI0NS************* 

T8I=T9+10. 

PBI»l«0e*P9 

T22T=T21+10. 

P221*0.9B*P21 


CALL  HXGC  6, PS  * V8  ,  TG  ,P9  ,  V9  ,T9  ,P21 * VP1 ,T21 ,P22, V22, T22> 

1  XI,  ANPI,P3I,  T8I  ,P22I, T28I, CPCM, CPHM ) 

SE8=SE9  +  C  CCPHH*ALGGCT8/T9>  )-C  1  .96 71 8*ALGGCP8 /P9> > ) 

SE22=SE21  +  C  C  CP CM+AL0CC  T2R/T2I  ) )  -  ( l‘ . 98  71 3*AL0  GCP21  /P22>  >> 

WRITE< IP, 10) 

WRITEC 12,  74) J 
I  F  C  K  R  )  9 1  6,  191  6,91  7 

917  G0  T0  918 

916  CALL  0UTPUTC  6) 

918  CONTINUE 

WRITEC 12, 656>P8,V8,T8 
WRITEC 12, 659>P22,V22,T22 
WRITEC  12,  668)  V9 

656  FORMAT C  IX,  5H  PB  =  F15.8,SH  V8  *F1S.8,SH  T8  *=F1S.8> 

659  F0RMATC  1 X,  5H  P22=?15.8,5H  V22=F15.8,5H  T22*F15.8> 

668  F0RMATC 1 X,  5H  V9  «=F1S.8) 

1916  CONTINUE 

WRI TEC  12, 10 ) 

IFCT22-540. IFVFFVLRM7  CONTINUE 
C  HX7  DATA  " 

REA0CI3, 69) J 

READC 1 3, 1 00 ) DH, TP, TS, SI G, ANRE1 
REA0CI3,  100)ANTUIC  J),C»ANrl,ANC,  FS 
REA DC  13, 100 )RA, RF, BX, BY, BPR 
READC 1 3, 100 ) W4 
WR1TEC12, 73) J 
WRITEC I2,25)DH 
WR1TEC12,P6)TP 
WR1TF.C12,27)TS 
WRITEC 12, 28) SI G 
WRITEC 12, 29)ANRE1 
WRI TEC  12,  30  )C 
WTEC 12,  31  )FS 
WRITEC 12, 32)RF 
•  WRITEC  12,  727ANTUICJ) 

WRI TEC  12, 10) 

READC 13, 100)X1,ANPI 

C* t ! t  J ! 91NITIAL  CONDITIONS************ « 

T7I  *T8+1  50  • 

P  71* 1 »02*P8* 

T23J»T22+150. 

P23I“0  »98*P22 
I FCKR)  200  5,  2006,2005 
200  6  GO  TO  1999 
200  5  CONTINUE 

CALL  HX 1  C  7, P 7  ,V7  ,T7  ,P8  ,V8,T8  ,P22,  V22,T22,P23,  V23,T23, 

1  XI,  ANPI,P7I, T71  ,P23I, T23I, CPCM, CPHM) 

SE7*SE8  +CCCPHM*ALOGCT7/T8))-C1.98718*AL0GCP7/P8))) 
SE23*SE22  *C  C  CP  CM*AL0  GCT23/T22) )- C 1 .98  718+ALO  6CP22/P23)) ) 


Wil  (12*10) 

VJR.il  j.<  12*  74)  J 
IM  KH)  919*  1  919*920 

920  G0  TO  921 

919  CALL  OUTPUT  (7) 

921  CONTINUE 

WRITEC 12*  775)P  7*  V7>  T7 
.  WRITEC12* 776)P23* V23*T23 
WR1TEC12*  777) V8 

775  FORMA  TC IX*  5H  P7  =  F15.8*5H  V7  =  F15.8*5H  T7  s=F15.3) 

776  FORMATC IX*  5H  P23=F15.8*5H  V23“F15.8*5H  T23=F15.8> 

777  FORMA T(  IX*  5H  V8  =  F15.8) 

1919  CONTINUE 

WRITE* 12*10) 

IFCKR)999* 1999*1001 

1999  SEM«=AMAX1  (SE7*  SE8*SE9*  SE10*  SE1 1*  SE12*SE1  3*  SE1  4*SEI  5*SE1  6* 

1  SE1  7*  SEt  8*  SEl  9*  SE20*  SE21  *  SE22*  SE23> 

TM»AMAX1  (  T7*  T8*  T9*  T1 0*  T1 1  *  T1 2*  T1  3*  T1  4*  T1  5*  T1  6*  Tl  7*  T1 8*  T1 9* 
1  T20*  T2l*T22*T23) 

DX«SEM/10. 

DY»540./10. 

1001  READC1V* 101)  J 
CALL  OUTPUT* J> 

GO  TO  1001 
999  STOP 
END 
END 

EOF  READ($MAIN$)  100 
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g  »  yi  fc  «  M  —  00)00 


SUB  KOn't  i..r  HM  <J*PAO,  VAG*  TAO»  PA  I *  VAI*1  A 1 *  Pbl  *  VBI *  TBI*  PC0*  VC0* 
1  TC«*Xi»ANPI*PASI* 7AGI*PC0I*TC0I*  CPCK*CPHN> 

TURBINE  DESIGN  PRO  GRAM-HEL  1  UM 

CRYOGENIC  TURBOMACHINE  SYSTEMS  WITH  JT  VALVE - GAS  AND 

/OR  VAPOR  PHASE  LOAD 

VERSION  15JULY  19  70  .  . 

DIMENSION  S*8)*AN*b>/  iX 1 C8 ) * AH1 <8 > * A Fl (8 ) 

DIMENSION  XPR<8  )* YPR ( 8  )  * ANP* 8 ) *  DP  1 C8)*ANPRI (8) 

DIMENSION  ANPR (8 ) * ANNU1 (8  >  *H I  2 (8  >» ANRE2C8 ) 

DIMENSION  DP2(S ) >  ANPk2  <  8 ) *  ANNU2C8) *  HS 6(  8 ) 

DIMENSION  AK23*B)*AK34(6)*AK45(8)*U<6) 

DIMENSION  AXTHSC6 )* AXP <8  )*  El < 8 > * ANS( 8 > 

DIMENSION  QL(8)*AE*8)*VX(8)*HY(8)*ALL(8) 

DIMENSION  CV.(8)*Hl.(8)*TV(6)*ANF(6)*AV(8) 

DIMENSION  CHI (8)*  CC2(8)*S2*8)*AN2(8)*AH2(8) 

DIMENSION  K1 <B)*K21 <8>*K3*3)*K«(8) 

DIMENSION  ANTU(8)*DHT(8)*AKVT(8)*ANTUI (8 ) * AKP*  5)*  TLX (5) 
DIMENSION  XPR1 <6)*YPR1 <8)» YPR2(3)*YPR3C8) 

DIMENSION  ALAM*8>*  TX(8) 

COMMON  S*  AN*  AX1 *  AH1 *  A  FI *  XPR*  YPR* ANP*  DP  1  *  ANPR l *  ANPR 
COMMON  ANNUI*H12*ANRE2*  DP2* ANPR2  *  ANNU2  *  H56*  AK23*AK3* 

COMMON  AK45*U*AXTHS*AXP*EI*ANS.QL*AE*WJS*HY*ALL*  CK*HW*  TW 
COMMON  ANF*AV*CH1  *CC2*S2*AN2*AH2*K1  *K2I*K3*K4 
COMMON  Tl  A,  PI  At  VB1  A*HI  A,  SI  4*  LI 
COMMON  Tl  5*PI 5*  VBI 5* HI  5*  SI  5 
COMMON  Tl  6# PI  6*  VBI  6# HI  6*  SI  6 
COMMON  DH*TP*SIG*ANREI*C*FS*RF*ANTU*TS 
COMMON  ANTI *DH2*  SI  Gl *  SI G2* W2* VT 1 
COMMON  DHT*AKWT*ANTUI*MAT*  TLX 
COMMON  BX*  BY* ANH* ANC 
COMMON  XPRI*YPR1*  YPR2*  YPR3 
COMMON  ALAM*  TX*  BPR*  RA 
FORMAT (IX* 10HERR0R  ST0P*2X*F1 5.8) 

FORMAT* IX* 4* FI  5.8* 2X)) 

FORMAT < I  X*  27HREDUCE  ANREI'  OR  INCREASE  RA) 

FORMAT* IX* 1 2HREDUCE  ANTUI) 

FORMAT* 1 x , » IHINCREASE  FS) 

FORMAT* 1X*9HREDUCE  FS) 

FORMAT* 1X*15HX I  IS  TOO  LARGE) 

IS*  1 

C  INITIAL  VALUES 
KTi  =0 

ANP*J)-ANPI 
TAO=TA0I 
PC0*PCG  I 

C*****vINI TIAL  X  VALUE  IS  CALCULATED  NEXT*** 

X»FS*< ( <C* 1 • )*BPR)+ANH+*2.*BY) ) 

XI  =  X 

PA0=PA0I 

TCO*TC0l 
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Ki <J)=0 
K2-0 
K3< J>-0 

ka:j>=o 

K5=0 

K21 <J)=0 
NN3=0 
Nl  *0 
N2=0 
N3=0 

PI-3.14159 
•  E«2. 718282 
C  H0T  SIDE 

C  EG. 14 

S<  J>«  SQRT<0 .90 6894*<  <DH*+2>/SI G) ) 

AN<  J>  =  <  4SI G>/<P I*<  DH**2 > ) 

AX  1 ( J>*<AN< J)*PI*DH*TP)+<2» *< 1 »-SI G>  > 

C  EG.  4 

200  ANP  < J  >-ANP  <  J> 

300  TA0=TA0 

42  PA0-PA0 

Kl <J)=0 
C  LG1 

I  F<TA0-4O •  > 50#  50#  60 
C  EG.  7 

50  VBA0-<2* 765/E**<0  »03*< <PA0**O  •333333>/<TA0/4O»>>  >  >*<TA0/PA0) 

60  T0  61 
C  EG. 8 

60  VBA0-2.863O92*<TA0/PA0> 

61  CONTINUE  „ 

C  EQ.  9 

UA0=  <2  <•  37888E-0  7)  *<  TA044O  •  643  ) 

VA0«< 1 2. 4ANRE1 * VBA04UA0 )/pH 
AH1  <J)*U  44.*W1*V8A0)/VA0 
AF1<J)*AH1<J)/SIG  . 

C  EQ. 13# ASSUME  X 

100  X«X 
XPR<J)-X-<2.*BX) 

YPR< J>B<  (X/FS>"<2»*BY>> 

C  EG.  16 

YPR1 < J>3  <YPR< J>- <  <  C- 1  * )*BPR)-ANH>/<2. *C  C- 1  • > ) 

A FI C=ANH*<XPR< J)*YPR1 <  J>  > 

Kl<J>*Kl«J>+3 
C-  TEST1 

IFtAFl C- < AF1 <J> +0.001 >  >101 #101 #10101 

101  I  F<AF1  O  <AF1  <J)-0. 001  >>10100#  150#  150 
10101  IFCK1 <J>-1 >2000*101 11*10102 

10111  X-X/l *5 
Kl <J>=0 
60  T0  100 
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10100  IF(KTl-4>10105/ 1010  A, 1010 4 

10105  IKK  l  <  J)-l  )?000*  10  1  1  5>  1011  6 
10115  11XA=X 

T1  YA  =  AF  C 
T1  ZA  =  AF l < J) 

101  It  , 

GO  TO  100 

10102  IF<KTI-4>I0106/ 10104/  10104 

10106  1 1  XB=X 
TlYb-AFI  C 
TlZB*AFl <J) 

K  T 1  =  3 

1010*  CALu  nERA<TlXA,TIXB/TlYA/TlYB/71ZA/TlZ8/TlDYAl/TlDYA2/ 

1  T1DYB3/T1DYB4/T1DYI/X/AF1C/AF1  <J>/KTl > 

KT1 =KT  1+1 
G0  TO  100 
C  E0.18 

150  YPR2<  J>c  YPR1 <  J>  *<  RA/2  •  > 

YPR3<J>=YPR1 < J>*RA 
AL1=YPR1 < JJ/24. 

KT1  *0 

K1 < J)*0 

AL2*BPR/12. 

AL3=YPR3<J>/<12.*RA> 

C  HOT  SIDE  PRESSURE  DROP 

DPI  A*  <  VA0**2>/VBAG 

DPI  P*370  •0E-06*  DPI  A*SQRT(  <  TP/DH>/ANRE1  > 

C  E0  +  24 

DPI <J)»ANP<J)*DP1P 
PACOPA1+DP1  <J) 

NN3»NN3+1 
C  TEST2 

474  IF<PA0C-<PA0-1 .0E-03>>476/  477/477 

477  1F<PA0C-<PA0+1 .0E-03) >478/ 478/476 

476  lF<NN3-$>  4761/  4762/ 4763  *  v 

4762  PA0A*PA0 
PA0CA»PA0C 
PA0ZAcPA0 
60  T0  4761 

4763  I F<NN3-6> 4765/ 4765/ 4766 

4765  PA0B*PA0 
PA0CB«PA0C 
PA0ZB»PA0 

4766  K31>NN3-3 

CALL  1TERA<PA0A/PA0B/?A0CA/ PA0CB/PA0ZA/ PA0ZB/ T3DYA1 / T3DYA2/ 
1  T3DYB3/ T3DYB4/T3DY1/PA0/PA0C/ PA0/K31 > 

60  T0  42 

4761  PA0  -PA0C 

GO  T0  42 
C  E0.26 


I 


I 
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4  . 


c 


478  CONTINUE 

UM I  =  ( 8 • S549  7E-04/ < TAB- TA  i  >  )  *  ( (<TA0**| .  643) - < TA I ** t . )/ j  .  64  3) 
AKm  =  <57.79E-03/(0.00  3SS*C  TAS-TAI  >))*<<<  <0.00 3S5*TAG>**1  .  t  42) 

*  “((0.00355*TAI)**1.642))/J .642) 

THM=<CAKl*il/S7. 79E-03)**C  1  ./.  642)  )/0  .00355 
PHM=  <  PAB+PA I )/2» 

NN3-0 

CALL  CPSUB(THM^PHMiCP) 

350  CPHM=CP 

C  EQ.29 

CHI <J)*W1*CPHM 

ANPR1C J)=(CPHH*UM1 ) /AKM1 

ANN U1 C J)=3. 66+<<<0.l04)/CCTP/DH)/<ANREl*ANPRl ( J) ) ) )/C 1 •♦( <0 .01 6) 
1  /( <  CTP/DH  >/< ANREl *ANPR1 (  J)  )  ) **0 « 8  )  )  )  ) 

H12(J)«(12.*ANNU1 < J)*AKM1 )/DH 
C  '  COLO  SIDE 
4000  TC0«TC0 

400  PC0=PCO 

C  LG2 

I FC TCO- 40 . ) 44 W 441 , 442 

441  VBC0* (2* 765/<E**<0  »03*C CPC0**O  *333333) /( TC0/ 40 • ) )  )  )  )*<TC0/PC0) 

GO  TO  443 

442  VBC0e2.863O92*<TC0/PC0) 

443  AH2CJ)*RA*AH1 (J) 

VC0*< I 44.*W1*VBC0)/AH2< J) 

C  EQ.41 

UC0* < 2* 37868E-0 7) *<TCe**0. 643) 

ANRE2(J)«(VC0*DH)/(12.*VBC0*UC0) 

DPM«<VC0**2)/VBC0 

DP2P«370.0E-06*DPM*SQRT<<TP/DH)/ANRE2<J)) 

DP2<J)»ANP(J)*DP2P 
C  EG.  45 

>*C0C»PBI-DP2CJ) 

K2«K2*1 
C  TEST3 

4862  l F C P C0 C“ C P C0  +0  • 00 1  )  )  4863>  4863*  4861 

4863  XF<PC0C>CPC0~O.OO1))4861*'483»483 

4861  XF<K2-5) 4881  *481*  4882 

481  PC0A-PC0 

T2YA-PC0C 
T2ZA«PC0 
G0  T0  4881 

4882  IF«K2-6)487*  48  7#  4883 

487  PC08»PC0 

T2YB»PC0C 

T2ZB»PC0 

4883  KH*K2*3 

CALL  I TERA<  PC0A,PC0B* T2YA# T2YB* T2ZA* T2ZB* T2DYA1 * T2DYA2# 

1  T2DYB3*  T2DYB4*  T2DYI*PC0*PC0C*PC0*KM) 

I F(PC0  )2000*  400  *  400 
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438 1 
C 

48  3 
C 

48  30 
4831 


550 

C 


C 


5570 

C 


557 

C 

556 


C__ 

560 


C 

558 


5750 

C 

575 

577 


* 

% 


PCVJ-Pf'.C 
GS  TO  *.0» 
tO  ••‘'6 
CONTINUE 
1  E3T  4 

I  F<PCGC“  10- >4830,4831  *48 31 
LRITE< 12, 3) 

G0  TO  9799 
CON  t  i  Nl'E 

UM2s<B.5549  7E-0  4/<TC0-TBI>)*<<<TC0**l  .643)-<TBI**l .643) >/l .643) 
AKM2«<  57.79E-O3/<O.OO355*<TC0-TBI)>>* 
i  <<(<0.0O355*TC0>**l.642)-<<0.00352*TBI)**1.642))  /J.642) 
TCM=(<AKM2/57. 79E-03)**< 1 ./• 642) ) /O .00355 
PCMS  <PB 1 +PCO )/2» 

K21 <J)*X21 < J) +K2 
K2=0 

CALL  CPSUB<  T CM*  PCM* CP ) 

CP  CM  =  CP 
EQ.49 

CC2  <  J  )  =  W1  *  CP  CM 
ANPR2<J)=<CPCM*UM2)/AKM2 

ANNU2< J)e3« 66+<<<0.l04)/(<TP/0H)/<ANRE2< J)*ANPR2<J) ))) 

/<1 .♦<<0.01 6)  /<  <  < TP/DH) /<ANRE2< J)*ANPR2< J) ))**0«8)))) 

EG  •  70 

H56<J)M12.*ANNU2<  J)*AKM2)/DH 
I F<X3< J)- 5000) 5570/ 5570.9999 
CONTINUE 

LG3  '  .  « 

K3< J)«K3C J)+l 

ir<CHUJ)-<CC2<J)-l  »0E-0  6) )  556#  557*557 
I  IF <  CHI  <  J)  - < CC2<  J)  +1  .D E-0  6) )  560  *  560  #  558 
EG.  50 

DELI=<TAI-TBI )*<  E**< ANTUI <J)*<1«*<CHI < J)/CC2< J) )) )) 
DELX=<0EL1-<TAI-TBI))/<<CC2<J)/CH1  <J) )-!.) 

DELY«  <  CC2<  J)/CH1 <J))*DELX 

TC0OTBI  +  DELX  -  ■ 

GO  TO  5750 

EQ.62  '  . 

Z*  <  TA I  -T  BI >/(l.-< ANTUI <J)/<I .f ANTUI < J) ) ) ) 

TC0C=<TBI+Z)-<TAI-TB1) 

GO  TO  5750 
EG. 56 

OEL 1  =  CTA I -TBI )*<  E**< ANTUI < J) *< 1 •-<  CC2<  J)/CHI <J)>) )* 

DELX=<DEH-<TAI-TBI))/<<CH1  <  J)/CC2<  J)  )*I  • )  ■  ' 

DELY=CCH1  (J)/CC2CJ))*DELX 

TCO C=TBI+OELX 

N  3=  N  3  ♦  1 

TESTS 

IFCTCfiC  <7CC~1  .0E-03))576»S77#577 
IF<TC0C-(TC0-*1  .0E-03))573*578*576 


r*' 


1 

:  i 

t  '  *’  1 

)  244  ^  ’  ■<  . 

* 


576  IF(N3-5>5761/5762/5763 

5762  TC0A-TC0 
TC0CA*TC0C 
TC0  ZA=TC0 
60  T0  5761 

5763  7  F<N3“6>5765/  5765/  5766 

5765  TC0B- TC0 
TC0CB--TC0C 
TC0ZB=TC0 

5766  K31*N3-3 

CALL  1TERACTC0A/ TC0B/ TC0CA/ TC0  CB/ TC0ZA/ TC0ZB/  T3DYA1 / T3DYA2 
1  T3DYB3/  T3DYB4/  T3DYI/  TC0/TC0C/TC0/K31  > 

G0  T0  4000 
5761  TC0  *TC0C 

'  60  T0  4000  , 

578  IF<CH1<J>-CCC2CJ>-1  .0  E-0  6>  >  5781  /  5782/  5782 

5782  IF<  CHI  CJ>-CCC2<J>  +  l.OE-0  6>  >5783/ 5783/ 5784 

5783  TA0C=TBI+Z 
GC  T0  5785 

5781  TA0C*TAI+DELY 
60  T0  5785 

5784  TA0C*TAI+DELY 

5785  N4=N4+1 
N3*0 

C  TEST6 

IF<TA0C-<TA0-1.OE-O3>>  57600  /  57700  /  5  7700 
57  700  I  F(TA0C-<TA0+1 .0E-03)  >  57800/  57800/57600 
57600  IF(N4- 5)57610/ 57620/ 57630 
57620  TA0A*TA0 
TA0CA»TA0C 
TA0ZA*TA0 
60  T0  57610 

57630  1  F<N4-6>  57650/  57650/  57660 

57650  TA0B=TA0 
TA0CB«TA0C 
TA0ZB«TA0 
57660  K4l  -  N4"3 

CALL  1TERACTA0A/TA0B/ TA0CA/ TA0C8/TA0ZA/TA0ZB/ T4DYA1 /T4DYA2 
1  T4DYB3/T4DYB3/T4DYI/TA0/ TA0C/TA0/K41 > 

60  T0  300 
57610  .  TA0  =TA0C 
G0  T0  300 

57800  I  Ft CHI  ( J>- <CC2< J>- 1 »0E“06>  >5791 / 5790/ 5790 

5790  1FCCH1 CJ>-CCC2CJ>+1.0E-06>> 5792/ 5792/ 5793 

5791  EIA*< CHI  (  J>/CC2(J) 

El  C  J>“  <  1  •  -  <  E**C-ANTU1 { J>*(  1  •- EIA> >  >  > 

1  /( 1  •- <EIA*<  E**< -AN7UI (J)*(l»- E1A> ) ) > ) 

GO  T0  5795 

5792  E1(J>~ANTU1(J>/(1*  +ANTU1 ( J>  > 

G0  T0  5795 


♦ 


5793  EIA  =  CC2( J) /CHI ( J> 

El <  J>s ( 1 .-<E>*(-ANTUI<J)*(I • *  EIA ) )  )  ) 

I  /< I • *  C  EIA*  <  2**  ( -ANT  Ui <  J) * ( 1 .-EIA)) >) ) 

C  EQ  .66 

5795  CONTINUE 
Nl  *0 
N3=Q 
N3=0 
N4a0 

C  LG4 

TPR I  *  THM- ( ( Y  PR  1 < J>/( YPR1 ( J)*YPR2< J> «-YPR3( j) > >♦< THM-T CM ) > 
TPR2*THM-(((YPR1(J>*YPR2(J)>/(YPR1( J)+YPR2( J)+YPR3( J) > >* 

S  (THM-TCM) ) 

I F( THM  )  601 »  602*  602 
602  IF(THM-TLXCMAT)  >666#  666>  604 

604  IF( THM- 5*0. >667. 667* 6671 

60 1  CO  TO  201  1 

605  C0  T0  2012 

C  EO.  71 

666  AKP(1>«(1 ./(.l*(THM-TPRI>>>*(((49./2. >*(((. I*THM)**S>>- 

S  ( ( • |*TPRl >**2>  >-(  « I ./3.47)*(((.l*THM)*4>3.47>-((  .1*TPRI )**3.47 
S  >)>> 

C  EG* 74 

AKP(2>*( I ./CTHM-TPR1 > >*( ( (2. 765/2. > *(<THM**2>-( TPRt**2> > >- 
S  <<<THM**2.1 6)-tTPRI**2.16))/2.t 6)) 

GO  T0  668 
C  EO. 73 

6671  AKP< l  )■ 1 1 1 • 74 

C  EO. 76 

AKP(2>*92.25 
GO  TO  668 

C  EO.  72  <. 

667  AKP( l  >*  ( l •/( . I *(THM-TPRt >))*((*•((( . l*THM> ••2*  708)-( ( .1 *TPR1 > 

S  **2«  70S) >/2. 708) ♦( *9S5l  *(  ( ( . I *THM)**2>- (( . I 4TPR1 >442)))} 

C  EQ. 75 

AKP(2>«86.M 6.2S*(((THM+TPRt>/2>-2IO.>/330.) 

C  EO. 77 

668  AK23( J>aAKP(MAT)*0.9DQ6l*( (  S( J)-DH)/(S( J>- (DH/2*  >) > 

C  LG8 

IF(TPR1  > 660  j  66l »  661 

661  SFCIPRI  -TLX (MAT) >670? 670* 662 

662  1F<  TPRl  -  540 »  >  671 / 671  *  671 1  -  L 

660  GO  TO  2013 

663  GO  TO  2014 

C  EO. 79 

670  AKPC1  >■(! ./< •! *(TPR1 -TPR2)>>*( (( 49*/2» )*<<<• 1 *1 PR 1 >**2>>- 

S  ((. l*TPR2)**2)>-((l ./3.47>*((( .t*TPRl >**3. 47) -( < . 1 *TPR2>**3. 47 
1  >>>> 

C  E082 

AKP<2)«(1 ./(TPRl -TPR2>>*< ((2*765/2. >*< CTPRI **2> M TPR2**2> > > - 

*  ' 
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4  . 


mm 


C 

6711 

C 


C 

671 


C 

672 
C 

674 
676 

673 
678 
C 

675 


C 

6771 

C 


C 

677 


C 

679 

C 


C 

708 


*  ^CI!R1**2,,6>-(TPR2**2.1  6))/2.l  6>) 

G0  T0  672 
E0.81 

AKPC1  >=in  .74 
£0.84 

AKP<2>*92 .25 
G5  T6  672 
EQ.80 

1  Sir*  *  1;*:  "« •  « J  ?  •*  ,* ;  r**"** 

AKP{2>  =  86.  +  c6.25*((<TpRi+TpR2>/2,-2lo.>/330o 

.  AK34(s)>=AKP  (MAT) 

LG10 

IFCTPR2  >673*674,674 
IFCTPR2-TLX(MAT) >  675*  675*  676 
1  FCTPR2- 540.  >  677*  677*  6771 
G0  T0  201  5 
G0  T0  20 1  6 
EQ  .  8  6 

.  SIMS'  JMjMliSaMSISffiKMIMSSS'lS,.., 

EQ.89 

£0.88 

AKPCJ  )*| 1 1 . 74 

E0.9I 

AKP<2 >*92.25 
60  T0  679 
E8.87 

.  KBIISaiStoMilKSKSKMRnSRai;!,"” ' 
s3!3t:iSKiasrR.sss“"/“-’ 

EQ.  92 

AK45<J>*AKPCMAT>*0.93061*<CsCJ>-DH>/<S(J>-<DH/2.>>> 

■  >*<AL2,AK34< J>  >+<AL3/AK45( J)  , 

EQ.94 

ANSCJ>*ANPCJ>+1 . 

N2*0 

ALE*CANSCJ>*IS>/12. 

TLl*<TA|+TBl>/2. 

TL2*  C  TA0 ♦TC0  > /2 • 
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Y=X/FS 

AKL=  (  ( t  X  rY  )•  <XPR<  J>*YPK(J>  >>**■<<  C- 1  .  >*XPR<  J>*BPR>  >/l  44. 

BA=1 • / C  TL2-TLI ) 

AKBL=(  7.2?£-03)+BA*<  ( (TL2/1 . 58  5>**l  *58  5>-<(TLI/!  .  585>**l . 5B5> > 
C  EG • 1 00 

QL< J) - (AK8L*AKL)/ALE 
C  LGK 

1 F<  CHI ( J> - <CC2<  J)- 1 .0 E-06>  >90 1  *  902# 902 

902  I F( QH 1  <  J)  •  ( CC2 ( J ) +1  .0E-06)  >903/ 903/904 

C  E0.101 

901  ALAM( J)«OL< J>/< 3600 . *CH1 ( J>> 

AE< J>*  El < J  >/< I • -ALAM< J> ) 

C  TEST7A 

1F<AE<J>-1  0  9040/9041/9041 

9041  VR1TE< 12/ 4) 

G0  TO  9999 

9040  CONTINUE 

ANTU<J)»<AL0G((1 .-<AE<J>*<CHt < J>/CC2< J)))>/(1  .-AE<J>>>/ 

1  (1.-<CH1<J)/CC2(J)>>) 

GO  TO  7020 
C  E0.107 

903  ALAM< J>*OL< J)/<3600.*CHl < J>) 

AE<  J)*=El  <  *J>/(  1  »-ALAM< J)  > 

C  TEST7C  -  . 

IF(AE<J>-1  .>9042/9043/90  43 
90  43  WR1TE<  12/4) 

60  TO  9999 

9042  CONTINUE 

ANTUC J)CAE< J>/< l.~AE(J)> 

GO  TO  7020 

904  ALAM< JlfcQLC J)/< 3600«*CC2< J> ) 

AEt J)b£1C J>/< 1 . -ALAM<d> > 

C  TEST  78 

1F(AE<J>“1*  >90*>5z  9044/  9044 

9044  WR1TEU2/4) 

GO  TO  9999 

9045  CONTINUE 

ANTU(J>«<Al.OG<<  1  .-<AE<J>*<CC2<J>/CH1  <J>>>>/(  1  .-AE<J>>>/ 

1  <) .-<CC2<J>/CH1 <J>>>> 

C  LG1 4 

7020  K4<J>=K4< J>*1 

1FCCH1  < J )• <CC2< J>- 1  .0E-06)  >701/  702/  702 
702  1F(CN1  <  J)MCC2<  J>  +  1  .0E-0  6)  >703/ 703/  704 

C  E0.I10 

701  AXTHSC J>* (3600.*ANTU( J)*CH1 <J>  >/U<J> 

AXP  < J)*AX 1 <J>*AF1 <J> 

ANPC«<144.+AXTHS<J>)/AXP< J> 

anpci*aint<anpc> 

1F\  (ANP C-ANPCl  >>0.5  )  701 1/70  1  2  /  70  1  2 

7012  ANPOANPCm  • 


G0  T0  701  3 
7011  hNPC=ANPC1+1 . 

7013  CONTINUE 

G0  70  70  5 
C  EQ.116 

703  AXTHSt J)«f3600.*ANTUt J)*CH1 tJ))/UtJ> 

A.'Pt  J)*AXlt  J)*AF1  tJ> 

ANPC*tl  44.*AXTHStJ))/AXPtJ) 

ANPC1  eAlNTtANPC) 

I  Ft  t ANPC-ANPC1  >-0.5)  7031  #7032*  7032 

7032  ANPOANPCl  +  1  * 

60  *T0  7033 

7031  ANPC*ANPC1 +1  * 

7033  CONTINUE 
GO  TO  70  5 

C  EQ.113 

704  AXTHSt  J)«t3600.*ANTUtJ>*CC2tjn/UtJ> 

AXPtJ>*AXl  t.J)*AFl  tj> 

ANPCS  1 1  44**AXTHSt J) ) /AXPt J) 

ANPC1 -AlNTtANPC) 

I  F  t  t  ANP  C- ANP  Cl  >  -0  .  5)  70  41  #  70  42#  70  42 

7042  ANPC“ANPC1+1  • 

G0  TO  7043 

7041  ANPOANPCIO  • 

7043  CONTINUE 

C  TESTS 

705  K5*K5+1 

JFt ANPC“ANPt  J)  )  70  6*  707#  70  7 

707  IFtANPC-tANPtJ)+l  0)1000#  1000#  706 

706  ANPtJ)*ANPC 
GO  TO  200 

C  SIZE  £  WEIGHT 

'  C  EG  •  1  1  9 

1000  WXt J)=X 

HY<  J) *X/FS 

ALLf J>  =  1 t ANPt J)*TP)+tANSt J)*TS) ) 
XYZ*tX*Y)-tXPRtJ)*YPRt J>)  ** 
jm*XBRU)tBPK*tC-lO.  .  .. 

XYA»AF1 t J)*tRA+l 0*tl*-SIG) 

XYT1*XYZ+XYB+XYA 

XYT2=XYZ+XYB 

CWtJ)=t0.098*ANPt J)*TP*XYT1 )+t0.0 78*ANSt J)*TS*XYT2> 
HWtJ)=0O96*ttX*Y>-tAFl  tJ)*tRA+l  O  )  +  tX*Y/80  ) 

J  5 - 

ANFl’YPRl t J>**2 

ANF2*  AXPt J)/t ANH+YPR1 t J) ) 

ANF3*H12t  J)*ANF2*ANF1 
ANF4=tStJ)-DH)/tStJ)-tDH/20> 

ANF5*XPRt  J)*TP*0. 93061  *ANF4 
~ANF6="3 . *ANP  < J)*AK23t J)*ANF5 


249 


ANF7=ANF3/AhF  6 
C  EQ.12S 

ANF< J)«l »/t I .+ASF7) 
IFCANF(J>-.4)1970#1971#1971 

1971  IFCANF<J)-.fe)J972#i972#1973 
1970  WRITEU2#5> 

GO  TO  1972 
1973  WRITEC 12# 6) 

1972  CONTINUE 

A VC J)sAXTHSC J>/ <  <  V.X  <  J)*HY ( J)*ALL ( J> >  / 1  728.1 
GO  TO  9999 

2000  WRITEC 12# 1>THM 

2001  WRITEC 12. 1 )PHM 

2002  WRI TEC  12* 1 )PHM 
•2003  WRITEC 12# 1>THM 
6007  WRITECI2#1>TCH 

2008  WRI TEC  12.1) TCM 

2009  WRI TEC  *2# 1  >PCM 

2010  WRITEC 12# 1 >PCM 

2011  WRITEC 12# 1 )THM 

2012  WRITEC 12# 1 >THM 

2013  WRITEC 12# 1 )TXC J> 

2014  WRITEC 12# 1 >TXC J) 

2013  WRITEC 12# 1 )TCM 
2016  WR I TEC  1 2#  1  ITCH 
9999  CONTINUE 

RETURN 

STOP 

END 

END 


SUBROUTINE  HX2C  J*P,,  .}*  v/AC*  7AU*PAI *  VAI * TAI*P l?i *  VB I*  »ul*.  *,** 

1  VB0*  T30*XI * ANP I  * PAD I*  TAOl  *PBf)  I*  TB0I  >  CPCM*CPHM  ) 

C  TURBINE  DESIGN  PRGGRAW-HE.L  IUM 

C  CRY0GENIC  TURBOMACHINE  SYSTEMS  WITH  JT  VALVE - GAS  AND 

C  /0R  VAPOR  PHASE  LOAD 

C  VERSION  16JULY  19  70 

•DIMENSION  S  <  8  )  *  AN  (  8  )  *  AX  1  C8)*AHl  (8)*AF1  (6) 

DIMENSION  XPR<8)*YPR(8)*ANP<8)*  DPI <6)*ANPRl (8 ) 

DIMENSION  ANPR<8 )/ ANNU1 <8  )*H1 2<6)* ANRE2<8) 

DIMENSION  DP2 C6 )*ANPR2C8)* ANNU2 <8)*H56<6) 

DIMENSION  AK23<8)*AK34(8 )  *  AK45( S ) *  U(6  ) 

DIMENSION  AX i'HS  <8  )  *  AXP  (6  )*  El  (8  )*  ANSC8  ) 

DIMENSION  QL<8 )* AE<8  > * WX<  S)  *HY<8  )*  ALL  (8  ) 

DIMENSION  CW<8  )*HV.'<8 )  *T  W(  8  )*  ANF(8 )  *  A  V(8 ) 

■  DIMENSION  CHI <6 )* CC2<6 )* S2<3 )* AN2C8) *AH2<8 ) 

DIMENSION  K1 (8)*K21 <8 ) *K 3 < 8 ) * K4< 8 ) 

DIMENSION  ANTU(8)*DHT<8)*AKWT(6)*ANTUI<8)*AKP<5)*TLX<  5) 
DIMENSION  XPR1 (8)*YPR1 (8)* YPR2<8)*YPR3(8) 

DIMENSION  ALAM (8)# IX (8) 

COMMON  S* AN* AX  1 *  AH  1  *  A  Fl  *XPR*  YPR*  ANP*  DP  1  * ANPR1 *ANPR 

COMMON  ANNU1 *H1 2* ANRE2*  DP2* ANPR2* ANNU2*  H56*  AK23*AK34 

COMMON  AK  A 5*  U* AX TH S* AXP *  EI*ANS*OL*  A E*  WX * HY* ALL#  C  W*  H  W* T  W 

COMMON  AN F*  A  V*  CHI  *  CC2*  S2* AN2*AH2*K1 *K21 *K3*K4 

COMMON  T 1  A# Pi  4#  VBl  4* HI  4*  SI  A*  Wl 

COMMON  T1  5#  PIS*  VBl  5,  HI  5/ SIS 

COMMON  T 1  6*  P 1  6*  VBl  6*  H 1  6*  S 1  6 

COMMON  DH*TP*5I  G*ANRE1*C*  FS*RF*ANTU*TS 

COMMON  ANT1*.DH2*  SI  Gl*  SI  G2*  W2*  WT1 

COMMON  DHT*  AKWT*  ANTUI*MAT*  TLX 

COMMON  BX*  BY* ANH*ANC 

COMMON  XPR1 * YPR1 *  YPR2*  YPR3 

COMMON  ALAM*  TX*  BPR*  RA 

1  FORMAT (1X*1OHERR0R  STOP* 2X* FI  5.8 ) 

2  F0  RMA  T(1X*4(F15*8*2X)) 

3  FORMATS  1  X*  1 1  HIM  CREASE  RA) 

A  FORMAT* IX*  39HREDUCE  BORDER  DIMENSIONS  OR  INCREASE  TS) 

5  FORMATUX*  1 1HINCREASE  FS) 

6  FORMATUX*  9HREDUCE  FS) 

7  FORMATUX*!  5HX I  IS  TOO  LARGE) 

12*1 

C  INITIAL  VALUES 

KT1  *0 

C******INITIAL  VALUE  OF  X  IS  CALCULATED  NEXT*** 

X«FS*( <  <  C- 1 . ) *BPR) +ANH* ( 2*BY) ) 

XI«X 

PBOsPBOI 

TBB«TBOI 

TA0=TA0I 

ANP( J)*ANP 1 

K1<J)*0 
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K2-0 
K3CJ>=0 
K4< J?=0 
K5a0 

K21 ( o)«0 

K31  =0  '  • 

N1=0 
N2C0 
N3-0 
N  4=0 
NN3=0 

Pl*3. 14159 
E*2»  718282 
C  H0T  SIDE 

C  EG  .  1  * 

S< J>=  SQRT<0*906&94+<  < DH*+2 J/SI G> > 

AN  <  J  >  *  <  4 • *  5 1 G> / <  P 1 *  <  DH++2  >  > 

AX1<  J>  =  <AN<  J>*PI*DH*TP>*<2.*<1  .-SIG>> 

C  EG.  4 

600  ANP  <  J  >'*ANP  (  J  ) 

3000  TA0*TA3 

42  PA0«PA0 

K1 < J)=0 

C  LGl 

IF<TA0-4O.  >50  ,  50#  60 
C  £0.7 

SO  VBA0=<2*  765/E*+C0*03*<  <PA0*1>O«  333333)/<TA0/4O*  >  )> >*<TA0/PA0) 
60  10  61 
C  EQ.B 

60  VBA0=»2.8  63O92*<TA0/PAO> 

61  CONTINUE 

C  E0.9 

UA0=<2.37SBB£-O7>*<TAG**O.643> 

VA0  =  <  1 2  .  *  A  N  R  E 1  *  VBA0  *  UA0  ) /W 
AH1 < J)»<1 34.*W1*VBA0>/VA0 
AFl<J)*AKl CJ1/S1G 
C  EQ.  1  3^ ASSUME  X 

100  X=X 

XPRCJ)=X -<2.+BX> 

YPF:<J>=<X/FSX2.*BY> 

C  EQ .  1  6 

YPR1<J)=<YPR<J>-<<C-1.>*BPR>-ANH>/(2.*CC-1 •>> 

AF1 C=ANH*<XPR< J)*YPR1 <J)> 

K1 <  J )  =  K 1 <J)+1 
C  TEST1 

1FCAF1 C“<AF1 <J)40.001 ) ) 10 1 # 10 1  *  1010 1 

101  I FCAF1C~<AF1 < J)-0. 001 >>10100/1 SO* 1  50 

10101  1F(K1 (J) -1 )2000/ 10111*10102 

10111  X=X/1.S 

K1  C J)=0 
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C 0  T0  100 

10100  IF<KT1-4>10105/10104/10104 

10105  1F<K1(J>-1 >2000/10115/10116 

10115  T1XA=X 
tlYA*AFlC 
T1ZA=AF1 <J>  . 

10116  X*1  •  5*X 
G0  T0  100 

10102  IFCKT1-4)1010  6/  10104/  10104 

10106  TIXBSX 
T1YB-AF1C 
T1ZB=AF1(J> 

.  KT1-3 

10104  CALL  1TERA(T1XA/T1XB/ Tl YA/Tl YB/T1 ZA/T1ZB/T1 DYA1 / T1 DYA2/ 

1  Tl  DYB3/T1  DYB4/T1  DYI/X/AFl  C/AF1  (  J>/KT1  > 

KTl*KTl +1 
G0  T0  100 
C  EQ«  IS 

150  YPR2<J>*YPR1  <J>*<RA/2.> 

.  YPR3C J)«YPR1 < J)*RA 
AL1*YPR1 <J>/24. 

KTI*0 

K1CJ>»1 

AL2*BPR/12. 

AL3*YPR3<J>/C12»*RA> 

C  HOT  SIDE  PRESSURE  DROP 

DPI A*<  VA0**2)/VBA0 

DP1P»370.0E-0  6*DP1A*S0RT<<TP/DH)/ANRE1  > 

DPI t J>°ANP  < J)*DP1P 
C  EQ.25 

PA0C«PAI+DP1 <J> 

C  GET  CHI 

C  TEST2 

474  I F<PA0C- CPA0- 1  •OE-03)  >476/  477/  477 
477  IF(PA0C"  CPAO  +  1  .0E-03 >  > 478f/ 478/ 476 
476  IFCNN3- 5) 4761/ 4762/ 4763 

4762  PA0A*PA0 
PA0CA*PA0C 
PA0ZA=PA0 
GO  10  4761 

4763  I  FCNN3~6>4765/  4765/  4766 

4765  PA0BsPA0 
PA0CB-PA0C 
PA0ZB*PA0 

4766  K31*NN3-3 

CALL  ITERACPA0A/PA0B/PA0CA/PA0CB/PA0ZA/PA0ZB/T3DYA1/ T3DYA2 
1  T3DYB3/ T3DYB4/ T3DYI/PA0/PA0C/PA0/K31 > 

GO  T0  42 

4761  PA0  =PA0C 

GO  TO  42 


253 


t, 


t 


t 


m 

478  CALL  CPLUIU  T.V.  ,PA8#CP> 

C  FL.Pb 

3  500  Cr'  1  P  -  CP 

CALL  CPSUp<  I.Mi#Phi#CP> 

7800  CPI  7= CP 

NN3=0 

CP :i-  <  Cl  1  Pi  CP  l  7)/L. 

GAMIi-1  ./(  I  •  -  (  •  *'.9  i  44 7 4rf  7/ CP b )  ) 

C  EG*  30 

T l  7 T  l  2=  1  *-AN1  1  ■*  <  i  .  -  (PB1/PA0>**(CCAMB-1  .  l/GAMB)) 

TA0C=THI/<  Tl 7T1  2  > 

C  TEST3 

K2=K2+ 1 

28  62  IFCTAOC- (TA0+O .001  )  )26  63#  28  63#  28  61 

2863  IF<TAGC-<TA0-0.001 >)  26  61  #  28  3#  283 

2861  IMK2-2>2661#2Bl»2882 

281  7A0A=TAG  *  , 

T2YA=TA0C 
T2ZA=  TAO 
GO  TO  2581 

2882  lF<K2-6>287,267#2863 

28  7  TA0B=  TAO 

T2YB=TAuC 

T2ZB-TA0 

2883  KH=K2*  3 

CALL  ITEKACTA0A#  TAGB#T2YA#  T2YB#  T2ZA#  T2ZB#  T20YA1 #  T20YA2# 

1  T20YB3#T20YB4#T2DYI#  TAO#  TAO C#  TAG# KM) 

I F< TAO  >2000# 3000# 3000 
2881  TAO=  TAOC 

GO  TO  3000 
C  EQ.32 

283  CONTINUE 

UM1"<8. 5549 7E“0  4/ (TAG-TA1 ))+<(< 1AS+*1 . 643)-<TAl**l  .643))/!  *643) 

AKM1*<<  57.  79E-0  3)/< *00 355+< TAO - TAI >))*<<<< .OO355*TA0)**l . 642) 

1  •  << .0035S+TAI )**1 . 642) )/l  • 642) 

TH.V,  'AKM1  /S7.  79E-03)**1  ./•  642)2.00355 
PHM=CPA0»PAI)/2. 

K2=0 

CALL  CP  SUB f THM#PHM#CP) 

350  CPHM=CP 

C  EG. 35 

CHI  <J)=U1  4CPHM 

ANPR!  <J)=<CPHM*U?I1  )/AKMl 

ANNU1 (J)=3*  66*C<C0. 10 4)/< CTP/Oh )/ CANREl *ANPR1 < J) ) ) )/< 1 .♦< <0.01  6) 

1  / <  <  <  TP/DH ) /<ANRE1 *ANPR1 <J)))**0«8)))) 

H12< J)=<12.*ANMU1 < J)*AKM1 )/BH  ’ 

C  COLO  SIOE 

C  EG «  41 

4000  TB0-TB0 

C  EQ .  42 
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400  PBO=PB0 

C  LG2 

IFCTB0-4O  »  )441 , 441  ,  442 

441  VBB0=  (  2  •  765/(E**(0  .0 3*  C  <PB0**O  •  333333)  /<  TB0/4O  .)>)))* C  TEG/P L<fi  > 

G0  T0  443 

442  VBB0=2 • 8  63092*C  TB3/FPG ) 

443  UB0=(2  -8 6  E- 0  7 )  *  (  Tb0 *  '=0  •  6 4  3  > 

AH?.£  J>-aHl  ( J)*RA 

VB0=C1 44.*W2*VBB0>/AH2£J> 

ANRE2C J)= £ VB0*DH  )/£ 1 2. * VBB0*UB0> 

DPM=£ VB0**2>/VBB0 

DP2P=3  70  *  0  E-0  6*DPM*SQRT  (  <  TP/Dri  >/ANRE2CJ)  ) 

DP2<J>=ANP£ J)*DP2P 
PB0C*PBI-DP2(J> 

K32=K32  +  1 
C  TEST4 

48  62  I  F<PB0C-  <  PB0+O • 00 1  >  )  4863*  48  63,  48  61 

4863  IFCPB0C-<PB0-O.OO1  )>  4861 ,  483,  483 
4861  I F£K32-2>  488 1 , 4C 1  ,  4882 
481  PB0A=PB0 

PB0YA=PB0C 
PB0ZAkPB0 
60  T0  4881 

4882  I  F(K32~  6)  48  7,  48 7,  4883 

487  PB0B=PB0 

PB0YB=PB0C 

PB0ZB-PB0 

4883  KMl=K32-3 

CALL  ITERA (PB0A,PB0B, PB0YA,  PB0YB, PB0ZA, PB0ZB, PB0DYA1 ,PB0DYA2, 

1  PB0DYB3,PB0DYB«,PB0DYI,P80,PB0C,PB0,KM1  ) 

IF<PB0  >2000,400,400 
4881  PB0=PB0C 

G0  T0  400 
C  E0.52 

C  TESTS 

483  IFCPB0C-1O.  >4831, 4832,4832 

4831  WRITEf 12, 3) 

G0  T0  9999 

4832  CONTINUE 

UM2=<<8. 5549  7E-0  4>/C  TB0-TBI >>*((<  TB0**1 . 643>-CTBI**l  .  643)) /I .64? 
AKM2=C <57. 79E-O3>/£.OO3S5*CTB0-TBl >>)*(<  <  £ .OG335*TB0>**1 . 642 j- 
1  <(.0035S*TBI>**1 . 642>>/l  . 642> 

TCM=<CAKM2/S7.79E-03>**1 ./. 642>/. 00355 
PCM*£PBI+PB0>/2» 

K31 CK3I +K32 
K32=0 

CALL  CPSUB(TCM,PCM,CP> 

550  CPCMsCF 

C  £0.54 

CC2  <  J  )  a  W2*  CP  CM 
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C 


557 

C 

556 

C 

560 

C 

556 

5750 

C 

575 
577 

576 
5762 


5763 

5765 


5766 


5761 

578 

5702 

C 

578  3 


C 

5731 


1  Kc..'  <.  >-50'."/‘>5570/  5570/99. >9 
uNF  K2tJ)*(  Cl  ‘  CM  *  U-l  21/AKM2 

ANNl)2tJ)a  3.  66+ <  <  (0  •  1 0  4)/<  C  TP/DH > /<  ANRE2 C  J)  *ANPR2(  J)  )  )  > 

/( 1 .+< <0.01 6)  /<<<TP/DH)/<ANRE2<J)*ANPR2(J)))**0.8)))) 

H  6<  J)  = i 12.*ANNU2< J1+AKM21/DH 

l.  0  J  •  .  . 

K  4  <  J  )  -  X  4  <  J  )  + 1 

1F<CH1<J)-<CC2< J)-l .0E-0  6))5S6/ 557/  557 
1F(CWI  <J)-CCCS<  I>  +  1  .OE-O  0)560  /  560  /  558 
EQ.  56 

TB0C=1.:I+C<  TA0-TAI)/<CC2<J)/CH1  CJ>>) 

GO  TO  5750 
E0. 62 

TB0C-TA0- < TAJ-TBI ) 

GO  TO  5750 
EQ .  59 

TB0C-  TBI  +  <<TA0-TAI  )/<CHl  (J)/CC2<  J)  )  ) 

N4=N4+5 
TEST  6 

I  FCTBOC-<TB0-i  .0E-03))  576/57  7.  577 

I FCTB0C-(T80+1.OE-O 31)578/ 578/ 576 

I F(N 4- 5)5761/ 5762/ 5763 

T80A=TBC 

TB0CA=TBOC 

T80ZA=TB6 

« 

GO  TO  5761 

lF(N4-6)5765/S765z  5766 

TBOB-TB0 

TB0CBMB0C 

TB0ZB=  TB0 

K41 =N  4- 3 

CALL  I TERA < TB0A z TB0B/TB0CA/ TB0C8/ TB0ZA/  TB0ZB/  T3DYA 1 / T3DYA2/ 

1  T3DYB3/T3DY84/T3DYI/T80/TB0C/T80/K41) 

GO  TO  4000  *  . 

TOO  - TBOC 
GO  T0  4000 

I  F <  Cr‘ !  < J)-(CC2<  J)-l  .0E-06)  )  5781/5782/5782 
1  F  C  Cl ! I  C.»)-<CC2<  J)  +  l  .0E-06)  1578  3/  578  3/  5784 

EC  •  &  3 

ANTUI CJ ) ~  <  < TAB*  TA I ) /<  TA0-TBI ))/<! .“<<TA0“TAI ) /<  TA0-TBI ) ) ) 

El  C  J)=ANT'JI  <J)/<1  .+ANTUKJ)) 

C;':  T0  578  5 
E0.57 

AN  TUI < J)~ AL0G<  <  TA0-T80 )/<TAI-T8l) )/<l  • •<  CHI  <  J) /CC2< J) ) ) 
EIA=C>'1  <  J1/CC2CJ) 

E I C .  J  ) - ( i  .-<t**<-ANTUI<J)*<l.-EIA>))) 
l  /<l .-<EIA*<E**<-ANTUI<J)*<1  .-£IA>)  >>) 

GO  TC  5  76  5 
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C  FO .  60 

5784  ANTUI C J)=ALOGCCTA0-TB0)/CTA1-TB1 >>/Cl . - C CC2C J) /CHI CJ>>> 
E1A=CC2C J)/CH1 C J> 

EU  J)=C1  .-CE**C-ANTUCJ>*C1  .-E1A>>>> 

1  /Cl .-CE1A*CE**C-ANTUCJ>*C1.-E1A> ) >>> 

C  E0.65  .  • 

5785  CONTINUE 
N 1  =0 
N2=0 
N3-0 

N  4=0 

T1PR=THM-  (CYPR1 CJ1/CYPR1 C J5+YPR2C J>+YPR3C J) > )*C THM-TCM) > 

-  T2PR=THM-  C  <  CYPR1 C J)+YPR2C J) )/CYPRl CJ)+YPR2CJ)+YPR3CJ>) )* 

1  < THM-TCM >  ) 

C  GET  PLATE  THERMAL  CONDUCTIVITIES 

C  LG4 

IFCTHM  >601/  602j  602 
602  I FC  THM-TLXCMAT)  >666/  666/  604 

604  IFCTHM- 540.  >667,  667/  6671 

601  G0  T0  2011 

605  G0  T0  2012 

C  EQ  •  70 

666  AKPC1 >«C1 ./Cl .*CTHM-T1PR)>>*CCC49./2.)*CCC . 1*THM>**2>- 

1  CC.1*T1PR>**2)>>-CC1./3.47>*CCC.1*THM)**3.47>-CC.1*T1PR>**3. 47 

2  >>>> 

C  EQ.73 

AKPC2)=C1 ./CTHM-T1PR)>*CCC2.765/2.)*CCTHM**2)-C  CT1PR>**2 
1  )))-<<< THM**2« 1 6)-CTlPR**2»l 6) ) /2» 1  6)) 

G0  T0  668 
C  EQ .  72 

6671  AKPC1 )=1 1 1 .74 
C  EQ. 75 

AKP  C2) =92*25 
G0  T0  668 

C  F.Q.71 

667.  AKPC1 >=Cl ./Cl .*CTHM-T1PR>))*C-  <  C  C  C  .  1  *THM>**2  .  708  >- 

1  CC.l*TlPR>**2.706>>/2.708>+C9»551*CCC.l*THM>**2>-C  <• 1* 

2  T1PR>**2>)>> 

C  EQ  •  74 

AKP  C2)=86.+C6«25*CCCC  THM+T1 PRJ/2*  J-210* )/330. >  ) 

C  EQ  *  7  0 

660  AKS3CJ*- AKP CHAT  )*Q.930  61*CCSCJ)-DH>/CSCJ>- CDH/2.))) 

TX  C  J  > - C  T  HM+TCM ) /2 . 

C  LG7 

IFCTXC.j)  >660  /661*  661 

£63  1  F  C  TX  C  J ) -  TLX  C  MAT)  >  670  /  670  /  662 

662  IFCTXCJ)-540.:  671/  671/  67J1 

660  G9  T0  2013 

6  63  ■  GO  70  20 1  4 

C  EQ .  78 
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6  VO  AKP  C  1  )-  <  1  . /(  I  .  •;  (  }!  r  k-T;v  ■:>  >  )  *  C  ( (  4V  .  /2  . )  *  C  C  C  .  1  *T1PR>*$2>  - 

1  C  C  . !  *  T  2P  ;•<>*  + 2  >  )  )  -  C  C 1  ./3. 4  7)  *(  (  C.  1+T1PR>**3.  47>-  <  <  .  1  *Tf>PR> *->2 . 47 

2  )  )  )  ) 

c  EC- .  a  i 

AKP  C  2 )  -  ( 1  ./CT1PR-T2PR>  >*C  C  C2. 765/2.  >  *  C  C  T1  PR**2>  -  C  C  T2P  R )  *8 
1  )>>-((  (TIPt^  +  2. 1  6>-tTP'»K**2.  !  6)  )/2. 1  6)  ) 

G3  TC  672 

C  EC. 80  • 

6711  AKPC 1 >=1 1 1 . 74 
C  E0.S3 

AKP  C2)=  92. 25 
GO  T0  672 
C  EO .  79 

671  AKPC1 >  =  C1 ./<1  .*CT1PR-T?PR>)>*<-  C { C C . 1 *71  PR ) **2.  708 ) - 

1  CC  .  1 *T2PR>  +  *2. 708) >/2. 706)  +  <S. 551  * C  C C .  1  *T1  PR ) **2> - C C  .  1  * 

2  T2PR>**2>>>> 

C  E0.S2 

AKPC2>  =  86.  +  C  6.25*<  C<  CllPR+T2PR>/2. >-2lO.>/330.>  > 

C  EG. 8 A 

672  AK34C J>=AKPCMAT> 

C  LG10 

I  FC TCM  >  67 3*  674#  674 

674  I FCTCM-7LXCHAT) ) 675#  675>  676 

676  I FCTCH- 540 •  )  677#  67  7*  6771 

673  G3  TO  2015 

678  GO  T0  20  1  6 

C  EO  .  8  5 

675  AKP Cl >=C1 ./Cl .*CT2PR-TCM>> >*C CC49./2. ) * < C < . 1 *T2PR >**2> - 

1  CC.l*TCM>**2>>>-CU  ./3.47>*CCC.1*T2PR>**3.47>-<C.1*TCM>**3.47 

2  >>>) 

C  E0.88 

AKP <2 >= Cl ./CT2PR-TCM))*CCC2.76S/2. >*C  CT2pR**2>“ CCTCM>**2 
1  >>  >-CCCT2PR**2. 1 6)- CTCM**2.1  6>)/2.  1  6)> 

G0  T0  679 

C  EQ.6  7 

6771  AKPC1 >=111 .74 

C  EQ  .  90 

AKP  C2>  =  92*  £5 
G0  T0  679 

C  EO .  8  6 

677  AKPC1 >=C1 ./Cl .*CT2PR-TCN>>>*C-  C  C  C  C  .  1  *T2PR  >  **2.  708  >  - 

1  CC.  1*T«.)**2.  708)  >/2. 706>  +  C9.55l*CCC.  1 *T2PR)**2>-C  C.  1* 

2  TCM>**2>>>> 

C  EO  *  89 

AKPC2>  =  36.  +  C  6.25*CC  C  CT2PR+TCM  >/2. >-210. >/330. >  > 

C  Efl  •  91 

679  AK45CJ>=AKPCMAT>*0.93061*CCSCJ>-OH  >/CSCJ>-CDH  /2*>>> 

C  EO • 92 

UCJ>=1./CC1 ./HI 2C J>  >  +  CALl /AK23C J>  >  +  CAL2/AK34C J>  >  +  CAL3/AK4  5C J>  > 

1  +C1 ./H56CJ>>> 
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C  GF.T  ENDWISE  HEAT  LEAKAGE 

C  ECU  93 

ANSCJJ--ANPC J) H . 
iJ  2-0 

AL  E=  (  AN  S  C  J  )  *  T  S  >  / 1  2  » 

TL1=(TAI+TBI )/2. 

TL2=  (  TAG  +  Triw  )  /2  • 

AKL=C  C  CX+CX/FS)  )  -  (XPR  C  J)  *YPRC  J)  )  )  +  C  <  C“  1  .  )*XPRCJ)*BPRJ)/1  4  4  . 
BA=1 . /CTL2-TL! ) 

AKBL=  < 7.27E-03)*BA*C (C  CTL2>**1 . 535>-C C TL 1 >**1  • 585) >/l  • 535) 

C  EQ.99 

GLCJ>=CAKBL*AKL)/ALE 
C  LG1  3 

.  I  Ft  CHI  (  JJ-CCC2C  J)-l  .OE-O  6)  >901,902#  902 

902  I  Ft CH! C JJ-CCC2C J>+1 .OE-O 6) ) 903# 903# 90  4 

C  E0 . 1 00 

901  ALAMC J)=0L< J ) /C  3600  . *CH  1 C J)) 

AEC J)=E1C JJ/Cl .-ALAMC J) ) 

C  TEST  7A 

IFCAECJ>-1 .) 9001 #9002# 9002 
9002  WRI TEC  12*  A) 

G0  T0  9999 
9001  CONTINUE 

ANTUC J)  =  CALGGC  C 1  - - CAEC J) *CCH1 C JJ/CC2C J) > ) )/C 1 .-AECJ>))/ 

1  Cl  .-CCH1  CJ)/CC2CJ) ))  ) 

G0  T0  70200 
C  EQ.10  6 

903  ALAMC J>=0LCJ)/C  3 600* * CHI  C J) ) 

AEC J)=EI C J>/C 1 . -ALAMC J>) 

C  TEST7C 

IFCAECJJ-l .0)9003# 9002# 9002 
900  3  CONTINUE 

ANTUC J)=A EC J)/C 1 .-AEC J)) 

G0  T0  70200 
C  EO. 103 

904  ALAMC J)=GLC J>/C 3600. +CC2C J) ) 

AECJ)=EIC J) /Cl .-ALAMC J) ) 

C  TEST7B 

IFCAECJJ-l. >900 4# 9002# 9002 
9004  CONTINUE 

ANTUC J>  =  CALGGC Cl.- CAEC J)*CCC2C JJ/CH1 C J))J )/Cl  .-AEC J>>>/ 
1C1.-CCC2CJJ/CH1 CJ))J) 

C  LGl  4 

70200  K5CJ)  =  K5C  JJ  +  l 

IFCCH1  CJJ-CCC2CJ>-1  .OE-O  6)  )  701  #702#  702 
702  IFCCH1  C  J)  -CCC2C  JJ  +  l  .OF.-O  6)  )  70  3#  703#  704 

C  EQ.109 

701  AXTHSCJ)=C  CANTUC J)*CH1  C  J)  J/UC  J)  )*3600. 

AXPCJ)=AX1 CJ)*AFl C J) 

ANPC=  C 1 44.*AXTHSC J) )/AXP  C  J> 
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AjrrCi  -a I i  *  '  :  *r.> 

i  •  ■  ri  >-o.5>7on*  701  St  7012 

vois  a:>-c  .  >  1  :  . 

GO  Ti;  70  1  3 
70  11  AMi-'C»A:iP;i  +1 . 

70  l  0  COM  I  MU'- 
GO  IT)  7u5 

C  £0.113 

703  AXTHOf J)=C(ANTU(J)*CH1 < J> >/U< J) >*3600 . 

AXP( J)-AX1 C 0>*AF1 (J> 

ANPC- Cl 44.*AX7HSCJ>  >/AXP< J> 

ANPC1  =  A I N  T ( 3  MP  C  > 

I  FCCANPC-ANP  Cl  >-0.5)70  31  *70  32*  70  32 
70  32  ANPC=ANPC1  +  1  * 

G0  TO  70  33 

70  31  ANPC=ANPC1 +1  . 

7033  CONTINUE 
GO  TO  70  5 
C  EQ .  1  1 2 

704  AXTHSC J>=( CANTUC J>*CC2( J > >/U( J>>*3600. 

AXP( J>=AX1 ( J>*AF1 <J> 

ANPC“ ( 1 4  4*  *AX  THSCJ>>/AXP(J> 

AN  PCI  =AINT<ANPC> 

I  F(  (ANPC-ANPC1  >-0. 5>  7041 *  7042*  7042 

7042  ANPC=ANPC1  +  1  . 

GO  TO  70  43 

7041  ANPCC  ANP  Cl  + 1  • 

7043  CONTINUE 

C  TEST8 

70  5  I  FCAMPC-ANPC  J>  >  70  6*  70 7*  70  7 
70  7  I FCANPC- (ANP( J)  +  l . >  >1000# 1000 *  70  6 

706  ANP(J>  =  ANPC 

GO  T0  200 

C  SIZE  Z  WEIGHT 

C  EG  •  1  1  3 

1000  WX<J>=X 

HY  (  J  >  =  X  /  F5‘ 

ALL  v  J>  -  <  (  Af.P  CJ>*TP>  +  (ANSCJ>*TS>> 

XYZ-CX>CX/FS) >- CXPRC J>*YPR(J> > 

XYB-XPCC J>*BPR*CC-1 . > 

XYA-  A?  S  <J>*  (RA+1  .  >*n  .  -  Si  G> 

C.A  Jj  - ^0 .oy6*ANPC J>*TP*CXYZ+XYB+XYA>  >+<0 .0 73*ANS< J>*TS*<XYZ+XYB> 
M  V.  ( J  1  76*(  CX*HY  (J>)-((RA+1  .>  *AF1  (  J>  >  +  <  <X  :HY<  J>  >/8*  >  > 

TUC  J>-CV(  J>  J> 

AN-'l-YPKI  C  J)*  +  2 

ANF2=  AXPt J)/<ANH*YPR1 <J>  > 

ANF3=H1  2(.1)*ANF2+ANF1 
A N F 4= C  3C J>-OH>/(SC J>- (OH/2. >  > 

AN  I1  *.  -XPPC J)*TP*0 .93061  *  AN  F4 
A  N  r  '  ■  »  3 .  ••  .AN1 °  * .  j  >  *  A  K  2  3  <  J  >  *  A  N  F  5 


ANF7=ANF3/ANF6 
ANFfJ)=  1 */t 1 «  +ANF7) 

C  TZST9 

IFtANFtJ)*-»4>700!  *  7002*  7002 

7002  IFfANF  CJ)-.  6)7003*  7003*  7004 
7001  WRI TEt 12*  5) 

G0  T0  7003 
700  4  WRI TEt 12*  6) 

7003  CONTINUE 

AVtJ>=AXTHSt  J)/t  tWXt J)*HYt J)*ALLf J)>/1 728.) 

C  TURBINE  1  OUTPUT 

DHTf J)=CPB*TA0*ANTl*t l « - < (PBI/PA0 >**<  < GAMB-1 . )/GAMB)>  > 
G0  T0  9999 


2000 

WRITECI2* 

1  )  THM 

2001 

WRITE(  12* 

1  )PHM 

2002 

WRITEt 12* 

1  )PHM 

2003 

WRITEt 12* 

1  )THM 

2007 

WRITEt 12* 

1  )  TCM 

2008 

WRI TE( 12* 

1  >TCM 

2009. 

WRITEt 12* 

1  )PCM 

2010 

WRITEt  12* 

1  >P  CM 

2011 

WRITEt  12* 

1  )THM 

2012 

WRITEt 12* 

1  >THM 

201  3 

WRITEt 12* 

1 )TXtJ> 

2014 

WRITEt 12* 

1  )  TX  <  J  ) 

201  5 

WRITEt 12* 

1  )  TCM 

201  6 

WRI  TEC  12* 

1  )  TCM 

9999 

CONTINUE 

RETURN 

STOP 

END  ' 

FfID 

S  Ul>l  v  ^  L1 1  i  i  i  f1 

N  ) 


I !  Y  A  l  ,  DYA2,  DYB3,  DYBA,  DYi ,  M.iu 


i  1  i1. ! ,<  i  ( A  >  M  h  #  Yi t  ,  Y  H >  i  ' *  t  ; t 


C 

C  PRG  RAM  KLGUi  RI.5  INPUTSC FROM  MAIN)  --XA,  XB,  YA,  YB,  ZA,  ZB#  YI*  Zi 
C  AND  OUt^UTSCiO  MAlN)--.Yi 
IFCN-4)  10,20,  100 

10  CALL  I  t  r  ;t  l  (  XA,  XB,  YA ,  Y  B,  ZA,  Lb*  X  I  ,  DYA  1  *  DYA2,  DYB3,  DYB  A) 

N-N  +  l 
HEIURN 

2D  CALL  ITER2(XA,XB,YA,YB,ZA,  ZB,  X I  ,  DY  I,  DYA  1  *  DYA2,  DYB3,  DYB4,  YI  ,  ZI,N> 

N-N  +  l 
RETURN 

100  CONTINUE 

CALL  ITF.R3<XA,XB,YA,YB,  ZA,  ZB,  X 1  ,  DY  I ,  DYA  1  ,  DYA2,  DYB3,DYB4,  Yi>Zi,iO 
■  N=N+ 1 
RETURN 
END 
END 

SUBROUTINE  I  TER  1 CXA, XB, YA, YB, ZA, ZB, X I ,  DYA1 , DYA 2#  BYB3, DYB4) 

DYA  1  =  ZA- YA 
DYA2=YA-ZA 
DYB3=ZB- YR 
DYB4=YB-ZB 

C  WRITE(D,2TDYA1,DYA2,  PYB3,  DYB4 

2  FORMAT  (1X,4F15*7) 

I F( YA- ZA) 99, 99, 1 DD 

99  XI  =  XA+ ( (DYA 1 / ( DYA1 +  DYB  4) )*<XB-XA>  ) 

C  WR ITEC D,  2  )X  I 

RETURN  '  ■ 

1DD  XI=XA+<(DYA2/(DYA2+DYB3))*(XB-XA)> 

C  WR  1  TE(0, 2  >  X I 

RETURN 
END 
END 

SUBROUTINE  ITER2(XA,XB,  YA, YB, ZA, ZB, XI, DYI, DYA1 , DYA2, DYB3, DYB4, 
4YI,ZI,N) 

X  I A=X  I 

!FfYA-ZA)99,99,lD0 
99  DYl=ZI- YI 

C  WRITEtD, I )DYI 

I  FORMAT (IX, 2 FI  5.7) 

IF<DYI)200,20D,30C 
IDD  DYi=YI-ZI 

1F<DYI)4D0, ADD, 500 

20 D  XI=XB“C ( DYB4/CDYB4+DYA1 ))*(XB-XA)) 

C  .  WRITECD, 15X1 
RETURN 
3D0  RETURN 

40D  XI=XB- C(DYR3/CDYB3+DYA2))*CXB-XA)) 

RETURN 
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500  RETURN 
ENO 
END 

SUBROUTINE  ITER3CXA,  XB,  YA,  YB,  ZA,  ZB, X  I ,  DYI ,  OYA 1  ,  DYA2*  DYB3*  0Yt!4, 
*YI/ZI>N) 

I F<  YA- ZA >99/99/1 00 

99  '  DYI=ZI-YI 

C  VJRI  TEt  0  /  1  >  DY I 

1  F0RMATUX>2F1  5.7) 

I F<  DYI >200  *  200* 300 

100  DYI=YI -ZI 

C  WRI TE<0# 1 ) OYI 

I F<  DYI ) 400*  400*  500 
200  X I A=X I 

DYI=  ABS(DYI) 

XI=XlA-< (DYI/CDYI+DYAl >>*CXIA-XA>) 

C  WRITE<0* 1 )XIA*XI 

C  WRI TE<0* 1 ) DYI*  OYA 1 

RETURN 
300  XIA=XI 

DYI=  ABS(DYI) 

XI=XIA+((DYI/<0YI+QYB4> >*<XB-X I A ) ) 

RETURN 
400  XIA=XI 

DYI=  ABS<  DYI > 

XI*XIfi-<CDYI/(DYI  +  DYA2)>*CXIA-XA>> 

RETURN 
500  XI A=XI 

DYI-  ABSv'OYI) 

XI=XIA+<(DYI/(DYI+DYB3) >*CXB-XIA)> 

RETURN 

END 

END 


B0F  READ<  SNA  I  NS)  100 


s|:  cr i.'iL  rivu'Hj) 

C  Vt  b_  K  Lt  1  wtt  j:  LY  19  70 

DIMENSION  ^  C  a  >  >  AM  <8  >#  AX  1  (  8  ) #  AH1  ( 8  )#  AF1  (  8  ) 

DIMENSION  XPRC8)#YPRC8)#ANPC8)#DP1C8)#ANPR1C8) 

DIMENSION  ANPRC  8 ) #  ANNU1 C8 )#  HI  2  C  8 ) # ANR E2 ( 8 ) 

DIMENSION  DP2CC)#ANPR2C8)#ANNU2C8)#H56C8) 

DIMENSION  AK23C8)#AK34<S)#AK45C8)#UC8) 

DIMENSION  AXTHS C8  )  # AXP C 8  )  #  El ( 8  )  #  ANSC  8 ) 

DIMENSION  0LC8)#AEC3)# WXC8).HYC8)#ALLC8) 

DIMENSION  CUC8)#HWC8)#  TWC3>#  ANFC8)#AVC8) 

DIMENSION  CM  C8)#CC2C8>#  S2C 8 ) # AN2C8 ) # AH2C3 ) 

DIMENSION  K 1 C8)#K21 C8)#K3C3)#K4C8) 

DIMENSION  ANTUC8)#DHTC8)#AKWTC8)#ANTUIC8)#AKPCS)#TLXC5) 
DIMENSION  XPR1 CS># YPR1 C8)#YPR2C8)# YPR3CC) 

DIMENSION  ALAH C  8  )  #  TX  <  8  ) 

COMMON  S#  AN#  AX  1  #  AH1  #  AF1 #XPR# YPR# ANP# DPI#  ANPR1 #  ANPR 

COMMON  ANNU1 j  H 1 2/ ANRE2#  DP2# ANPR2# ANNU2#HS6#  AK23#AK34 

COMMON  AK4S#  U#  AXTHS/ AXP*  El# ANS#0L#AE#  IJX#  HY# ALL# C W#H Vi#  T W 

COMMON  AN E#  A  V#  CH 1 # CC2# S2# AN2# AH2#K 1 #K2l#K3#K4 

COMMON  T1  4#  P 1  4#  VB 1  4#  H 1  4#  SI  4#  VI 

COMMON  T1  5#P1  5#  VB1  5#H1  S#  SI  5 

COMMON  T1  6# PI  6#  VB I  6# HI  6#  SI  6 

COMMON  DH#  TP#  SI  G>  ANRE1  #  C#  FS#  RF#  ANTU#  TS 

COMMON  ANTI #  DH2#  SI G1 #  SI G2#  W2#  WT1 

COMMON  DHT# AKWT# ANTUI#MAT# TLX 

COMMON  BX#  BY# ANH# ANC 

COMMON  XPR1#YPRI# YPR2# YPR3 

COMMON  ALAM#  TX 

35  FORMATC 1 X#  8HWI DTH  X=Fl S»8#2X#  6H1NCHES) 

36  FORMAT C  IX  # 9HHEI  GHT  Y=Fl 5.8#2X# 6HINCHES) 

37  FORMAT  (IX#  9HLENGTH  L=F1  S«8#2X#  6HINCHES) 

38  FORMATC  IX#  12HC0RE  V.’El  GHT=  FI  5. 8# 2X#  3HL8S) 

39  FORMATC  IX#  14HHEADER  VEl  GHT=Fl  S*8#2X#  3HLBS)  • 

40  FORMATC IX# 1 3HTCTAL  WEI  GHT=  FI  5.8#2X#  3HLBS) 

41  FORMATC  IX#  NF-F1  5.8) 

42  FORMAT C1X#4H  AV=F1S«8) 

43  FORMATC  1 X#  3H  S=  FI  5.8#  3X#  4H  N=F1S.8#3X#  5H  AF1  =  F15.8) 

44  FORMAT  <  1  X#  SH  XPR=F1S«8#3X#  5H  YPR=  FI  S.  8#  3X#  5H  ANP*F1S.8) 

45  FORMAT C I  X #  5H  DP  1  =  FI  5.8) 

46  FORMATC  IX#  5H  CH 1  =  FI  5.8# 3X#  6H  Npftl  =  FI  S.S#  6H  NNU1  =  F15.8) 

47  FORMAT C 1 X#  5H  HI 2~ FI  5.3) 

48  FORMAT C I  X#  6H  NRE2= FI S. 8# 3X# 5H  DP2=F1S.8) 

49  FORMAT C 1  X#  5H  CC2=  FI  5*8#  3X#  6H  NPR2*  FI  S. 8#  3X #  6H  NNU2=F15.8) 

50  FORMATC  IX#  5H  H56- FI  5. 8#  3X,  5H  K23=  FI  S.  8#  3X#  5H  K34=F15.8> 

51  FORMATC  1  X#  SH  K45=  FI  S.  8#  3X#  3H  U=  FI  5.8#  3X#  7H  AXTHS«F15.8) 

52  FORMATC  1  X#  SH  AXP=  FI  5-8#  3X#4H  EI*F1  5.8#3X#  4H  NS*F1  5.8) 

53  FORMATC  1  X#  4H  0L=  FI  5.8#  3X#  3H  E=FlSc8) 

54  FORMAT C 1 X#  7*10) 

55  FORMATC  1 X #  5H  NTU=  FI  5. 8#  3X#  6H  NTU1  =  F15.8) 

56  FORMAT  C  1  X#  5H  AX  I  =  Ft  5. 5#  3X#  5H  AH  1  =  FI  5. 8#  3X.»  5M  AH2*Fl5.8) 


ST  F0RriATClX#4H  S2*F1 5. 8# 3X# 5H  AN2=F15*8) 

58  F0RMATC  1 X#  5H  DHT=  FI  5. 8#  3X#  5H  KWT=F15.8> 

59  F0RMATC 1 X#  6H  ALAM»=F1  5.8#  3X#  4H  TX=F15.8> 

60  FORMAT C 1 X *  6H  YPR1  =  F1  S.8#  3X#  6H  YPR2=F15.8> 

61  F0RMATC1X#  6H  YPR3=F1  5.8) 

WRI  TEC  12#  35) WX  C J) 

WRITECI2#36)HYCJ> 

WRITEC 12#  37)ALLC J) 

WRITEC  1 2/  38>CW<J> 

'  WRITEC 12#  39)HWC J) 

WRITEC  I2#40)TWCJ) 

WRITECI2#  55>ANTUCJ)#AijTUICJ) 

WRITEC  12#  41  )ANF(  J) 

WRITEC  12#  42>AVC  J> 

WRITEC I2#43)SCJ)#ANCJ)#AF1 CJ) 

WRITEC 12#  *4)XPRCJ)#YPRCJ)#ANPCJ) 

WRITEC  I2#45>DP1  CJ) 

WRITEC 12# 46JCH1 CJ>#ANPR1 CJ)#ANNU1 CJ) 
WRITF.CI2#47)H12CJ> 

WRITEC 12  #  48) ANRE2C J)#  DP2C J) 

WRITEC 12#  49 ) CC2C J) # ANPR2C  J) # ANNU2C  J) 

WRITEC  12#  50)H56C  J>#  AK23C J)#AK34C  J) 

WRITEC 12# 51 )AK45CJ)#UCJ)#AXTHSCJ) 

WRITEC 12# 52)AXPCJ>#  EICJ)#ANSCJ) 
WRITECI2#53)QLCJ)#AECJ) 

WRITEC 12# 56)AX1 CJ)# AH  1 CJ)#AH2CJ) 

WRITEC 12#  57)S2C J)# AN2C J) 

WRITECI2# 58)DHTC J)#AKWTC J) 

WRITEC  12#  59)A1.AMCJ)#TXC  J) 

WRITEC  12#  60 ) YP R 1  CJ)#YPR2C*J) 

WRITEC 12#  61 )YPR3CJ) 

WRITEC  1 2#  54)K1  CJ)#K21CJ)#K3CJ)#K4CJ) 

RETURN 

END 
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“Ij'i-.m.  «  iui*  CPi»UlU..  y;  .ZuJ, 

ip<24/>  n?i  .i)>qp<2«>  i  4> 
0PiM«5>  INPUT*  /CPDAT/  ) 

1  FORMA i (SI  7) 

N-  1 

NA-8 

13=5 

12=1 

UJ  10  KIP  = 1 * 3 

READ? I 3*  1  > C IPCM1P >>MIP=N> NA) 

N=NA+1 

NA=NA+3 

10  CONTINUE 

2  FORMAT?  717) 

N=  I 

NA=  7 

DO  I  1  KIP=1 *  2 

READ? I3>2> ?IT?MIP)>MIP=N>NA) 

N=NA+ 1 
NA=NA+  7 

11  CONTINUE 

3  FORMAT?  1  1F6.  3) 

D012  KCP=  1  >  23 

12  REA0?I3,  3)  ?CP  ?KCP>KCT)>KCT=  3>  1  3) 
CP  ?  1  >  1  >  » 1  .  90  6 

DO  1  3  KCP=  6*  23 

13  READ?I3>3)CP?KCP>2) 
IXT~IFIX?1000.*XT> 

JYP=IFIX? 1000. *YP> 
IF?XT-9.468)100>101>102 

101  D01O3  K=l>23 
IF?IP?K)-IYP )200>  201 >  201 

200  I F?K-2  3) 1 0  3*  202#  202 

103  CONTINUE 

201  KB=K 
KA=K-1 

IF?KA>210>211 >212 

211  ZCP  =  CF?1>3) 

G0  TO  900 

212  DPB=  FLOAT ? IP  ?KB>~ IP  ?KA> ) 

DPP=  FLOAT  ?IYP-IP?KA) ) 

DCP=  CP  ?KB> 3)-CP?KA,3) 

ZCP=  ?  DCP*?  DPP/DPB) >  <  CP  ?KA> 3) 

GO  TO  900 

102  D0300  LT=  3> 1 3 
IF?IT?LT)-IXT)301>302>302 

301  IF?LT-13)300>303>  303 

300  CONTINUE 

302  LTD=LT 
LTA=LT- 1 
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DTB=  FLOA't  (IT Ci.l 6)  -  II  :LTA)  ) 

DTT=  KLO(\T<  IXT-ITCL1A)) 

400  D04O1  LP= 1.23 
IFCIP<LP)“IYP)-h02>40  3>40  3 

402  IFCLP“23J  401  >  40  4>  404 

401  CONTINUE 

403  LPB=LP 
LPA=UP - 1 

IF(LPA)410>  41  1  >412 

411  DCPTB=CP(LPB>LTB)-CPCLPB>LTA) 

ZCP=C  DCPTB*C  DTT/DTB ) )+CP (LPB*LTA > 

GO  T0  900 

412  DPB=FL0ATC IP(LPB)-IPCLPA) > 

DPPr  FLOATCIYP-IPCLPA)) 

•  DCPTA=CPCLPA>LTB)-CPCLPA>LTA) 
DCPTB«CP(LPB>LTB)-CPCLPB>LTA) 
ZCPTA=CDCPTA*C DTT/DTB > )+CPCLPA>LTA) 
ZCPTB=C  DCPTB*(  DTT/DTB) )+CP(LPB>LTA) 
DCP=ZCPTB-ZCPTA 
ZCP=<  DCP*C  DPP/DPB ) J+ZCPTA 
G0  T0  900 

ICO  IFCYP-40.)500>  501  >501 

501  D0503  LP=6>23 

I  FC  IP  <LP  )  - 1 YP  )  504>  505>  505 

504  I FCLP-23  )  503>  50  6>  506 

503  CONTINUE 

505  LPB=LP 
LPA*LP-t 

I F<LP  -  6)510>  51 1>  51  7 
511  DTB*  FLOATC ITC3)-ITC2>) 

DTT*  FLOATC IXT-ITC2)) 

I  FC  DTT)  51  5>  51  6>  51  7 

516  ZCP=CPC6>2) 

GO  TO  900 

517  DPB*  FL0ATC  IP  CLPB)  ~IP  (LPA  )  ) 

DPP*  FLOATC IYP- IP  CLP A) ) 

DCPTA=CPCLPA>  3)-CP(LPA>2> 
DCPTB*CPCLPB>3)-CPCLPA>2) 
DTB=FL0AT<ITC3)-ITC2>> 

DTT=FLOATC  IXT-  ITC  2)  ) 

ZCPTA=CDCPTA*C  DTT/DTB)  )+CPCLPA>  2  > 
ZCPTB= C  DCPTB*C  DTT/DTB) ) +CPCLPB>  2 ) 
DCP  =  ZCPTB-ZCPTA 
ZCP=CDCP*CDPP/DPB))+ZCP7A 

GO  TO  900 

500  YP1*CC30./1 .106>*XT>-1  77. 

IYP1=I FIXCYPl > 

DCPl'O.  616-1  .906 
CP1=CCCYPI-10.)/30.)*DCP1 >+1.906 
DO  700  LPS 1  >  6 
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?«:  i 
70  c 
7D2 


1  F f  T;  (Lp)-IYP)701»  702*  702 
I  Ft  P- 6)700,  704,  704 
C0to7iMUF 
Lf'8=LP 
LpA=LP- 3 

- FcLPA )  710,  71  1*  712 

711  DTT=X  T- 6*89  4 
DT6=9. 468-6.894 
DCP=1 . 781-1 .906 
ZCP=DCP*(DTT/DTB>+i .906 
GO  T0  900 

712  XTB='YP  +  177.)*<1 .106/30.) 
DTT-XT-XTB 

.  DTB=9  *  468- XTB 

DPP=FLGATCIYP-IPCLPA)> 

DPB=FL0ATC IP  CLP  B)- IP (LPA) ) 

YP1 =YP 

CP1=CCCYP1-10.)/30.)*DCP1 )+l .906 
DCP=CP  CLPB*  3 >- CP (LPA/ 3) 
ZCPA=DCP*CDPP/DPB)+CPCLPA,  3) 
DZCP=ZCPA-CP1 
ZCP=DZCP*C  DTT/DTB)  +  CP 1 
60  T0  900 

C  ERROR  CODES 


202 

WR1 

TEC  12, 1001 ) 

G0 

TO  9999 

210 

WRi 

TEC  12*1101) 

G0 

T0  9999 

303 

WRi 

TEC 12,1 201  ) 

* 

G0 

T0  9999 

404 

WRI 

TEC  12* 1301) 

G0 

T0  9999 

410 

WRI 

TEC I2> 1 4) 

GO 

T0  9999 

50  6 

WRI 

TEC  1 2, 1 5) 

GO 

T0  9999 

510 

WRi 

TEC  12* 1 6) 

GO 

Tfl  9999 

51  5 

•aRI 

TEC  12, 1 7) 

GO 

TO  9999 

70  4 

WRI 

TEC  12#  18) 

G0 

TO  9999 

710 

WRI 

TEC  12*  19) 

C-0 

TO  9999 

1001 

F0RMATC IX,  7HK 

GT 

23) 

1  101 

FORMAT C 1 X#  6HK  1 

LT 

1  ) 

1  201 

FOR 

liiATC  1X,8HLT 

GT 

1  3) 

1301 

FORMAT  C 1 X,8HLP 

GT 

23) 

14 

FORMAT  C  1  X, 7HLP 

LT 

1  ) 

15 

F0RMATC IXjBHLP 

GT 

23) 

\  t 

F0RMATC  IX,  7HLP 

uT 

6) 

i  7 

FQF'-IATC  iX»SHI>T  i 

LT 

0) 

18 

FOR'iAT <  1  X*  7HLP 

G7 

6) 

19 

FORMA T<  1Xj>  7HLP 

LT 

!  ) 

900 

CONTINUE 

CLO  iE<  5) 

RET  JRN 

9999 

STOP 

END 
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